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ABSTRACT 



Passive solar heating systems can supply a major portion of a house 
heating load if properly designed. The four basic concepts used are dire 
thermal storage in wall or roof, solar greenhouse and convective loop. I 
applications some of these concepts will be combined to give better overa 
performance. Technical advances will make passive solar system function 
better. How far the technology advances to, will depend on how well pass 
solar systems become accepted. To further the use of solar energy an int 
approach will need to be taken. One which will combine the best of activ 
passive systems to produce a system that supplies ICO" of the house's hea 
load. 
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PREFACE 

Passive solar heating systems can supply a major portion of 
a house’s heating load if properly designed. The four basic concepts 
used are direct gain, thermal storage in wall or roof, solar 
greenhouse and convective loop. In most applications some of these 
concepts will be combined to give better overall performance. Technical 
advances will make passive solar system function even better. How 
far the technology advances to, will depend on how well passive solar 
systems become accepted. To further the use of solar energy an 
integrated approach will need to be taken. One which will combine the 
best of active and passive systems to produce a system that supplies 
100% of the house’s heating load. 
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INTRODUCTION 



Passive solar heating is defined as using solar energy to heat a 
building, without the use of non renewable energy. The heating is 
basically done by letting the sun's radiation into a building's interior, 
to be stored in some kind of thermal mass or to heat up the air space. 

Using the sun for heating and cooling goes back to "the ancient Greeks, 
who enunicated the basic principle of passive heating: face south, 
where the sun spends the winter; keep the winter winds away by embankments 
or vegatating; shade against the summer sun; let cooling be done by 
evaporation water; work with, rather than against, nature. 5A This 
project will only cover passive heating, but most passive systems 
can do some cooling. The four basic concepts are direct gain, thermal 
storage in a wall or roof, solar greenhouse and convective loop. 

The direct gain concept allows the solar radiation through south 
facing windows or clerestroies to be stored in the walls or floor of 
the house. Thermal storage in a wall or roof allows the solar radiation 
to come through south facing glazing and then directly onto thermal 
storage, with no direct solar radiation entering the living space. 

Masonry material or water are normally the thermal storage mediums. 

The solar greenhouse concept uses the solar radiation captured by a 
south facing greenhouse for heating. Convective loop is similar to 
flat plate air solar collector, except the air is circulated by the 
thermosiphon process instead of by a fan. Many houses use a combination 
of the above concepts to make up their passive solar heating system. 

These concepts will be covered in detail later. 

A house that uses a passive solar system for heating will have to 
be sealed as well as possible to reduce the^ inf iltration heating load. 

In this time of dwindling energy resources any house should be sealed 
to insure the minimum heating load. Generally the economics of passive 
solar systems will not be addressed in the project, because it is 
difficult to assess the cost of the system. The economics analysis 
that has been done has sho\>m that passive solar heating is economically 
attractive. 

How the sun's path across the sky changes with the seasons, aides 
the operation of passive solar heating systems. During the winter 
months the sun is low in the sky, which is good for the collection of 
solar radiation on a vertical surface. During the summer months the 
sun is high in the sky, which will help prevent the collection of solar 
radiation, on a vertical surface. Thus the sun's position in the sky 
directly enhances passive solar systems operation, by allowing good 
collection during the winter months and poor collection during the summer 
months. The farther north a house is located the lower in the sky the 
sun will be, which will improve the collection of available solar 
radiation on a vertical surface. For a house to have an effective 
passive solar system it will have to have a wall that is within 20 
degrees east or west of south, which has no obstacles in the path of 
the solar radiation to the wall area. The house shape that will provide 
the optimum collection of solar radiation is a form elongated along the 
east -west direction. 

Exhibits 1 , and 2 , will give some idea of how well passive solar 
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(Based on national averages; not to be 
used for design purposes) 
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ANNUAL MEAN DAILY SOLAR RADIATION 
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bright enough to cast a shadow across a 
surface, divided by the number of hours the 
sun is above the horizon, 

ANNUAL MEAN PERCENTAGE OF POSSIBLE SUNSHINE 




heating during heating season : Excellent sunshine is combined with high 
heat demand in area 1, moderate heat demand in area 2, and low heat 
demand in area 3. Good sunshine with high heat demand in area 4, 
moderate heat demand in area 5, and low heat demand in area 6. Fair 
sunshine is taken from TRW, Phase 0 Study, “Solar Heating and Cooling 
of Buildings," Nation! Science Foundation, 1974 

EXHIBIT 2 
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systems will perform in different parts of the country. 3 These 
maps show that passive solar systems could be used throughout the 
United States, with the worst performance in area seven of 
exhibit 2- 

An important factor to consider when designing a passive solar 
system is the microclimate of the area the house is located. Often 
the weather patterns surrounding one house will be different from a 
neighbors house. These factors may be natural wind breaks, large 
bodies of water, mountains, valleys, open plains, landscaping and or 
proximity to surrounding buildings. l^Jhatever the microclimate is it 
should be considered when designing a system. Another thing to consider 
is the benefit gained from broad leaf trees. These trees will be bare 
during the winter months, thus allowing solar radiation through, and 
full of leaves during the summer months providing shade to keep the 
house from overheating. The biggest problem with passive solar systems 
is the temperature fluctuations. For a passive system to work there 
has to be changes in temperature, but if the system is improperly 
designed the temperature changes can be too great. Initially, the lack 
of adequate design parameters resulted in some poorly designed systems, 
which had large temperature fluctuations, and has given passive solar 
systems a bad name. The expericene gained from the initial design 
attempts has produced better design procedures, but the procedures are 
still too complicated for general use to be adopted. 

The problem of overheating has resulted in the use of fans to 
transport the excess heat to a separate storage medium. The use of a 
fan in passive systems has created a conflict among designers. The 
purest position is that the use of an electric fan makes the system an 
active one. Others try to put a coefficient of performance (50) limit 
on when a system is passive or active. As far as this project is 
concerned, as long as the basic system falls within one of the passive 
concepts, then it is a passive system. If a fan can increase the 
system's performance then it is the logical thing to do. The systems 
with a fan will be referred to as hybrid passive solar systems. 

People who have lived in houses with passive solar systems, that 
operate correctly, have liked the interior environmental conditions. 
Passive systems do take some time to learn how to operate, to get the 
most comfort, but once the system is learned the interior areas can 
be kept comfortable. A factor to consider with systems that use 
radiant heat, is the interior air temperature can be lower than 
conventional systems and still be comfortable. The reason for this is 
the radiant heat from the wall or floor can be felt as warmth. Also 
to be considered with a passive system is that there is no noise from 
mechanical equipment . 
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DIRECT GAIN 



The direct gain passive solar system is the simplest approach when 
using solar energy. The system simply allows the solar energy to enter 
through south facing windows and be absorbed in the mass of the house. 

The reason this system works is solor energy (short wave raidation) will 
pass through the glass, with very little loss in energy, and be absorbed 
by the mass of the house, which will then emit long wave radiation that 
will not pass easily through the windows. The windows will be a large 
source of heat loss due to the inside-outside temperature difference and 
the window heat transfer coefficient (U) . A method used to overcome this 
heat loss is to use movable insulation, to be covered later, to cover 
the windows when there is no solar energy available, mostly at night. 

A direct gain system will need an overhang to prevent the summer sun, 
which comes from large angles from the horizontal, from coming into 
the house. 5B 

The problems with direct gain systems are large temperature swings, 
expensive thermal storage mass, strong directional daylighting, glare, 
and ultraviolet degradation of fabrics. There will also be temperature 
swings of 15 to 22 degree fahrenheit during the winter months, which 
can be unacceptable to some families. These large temperature swings 
are also an indication that the system is not working very well. The 
movable insulation can help the situation, but there is still the 
overheating which is energy that could be stored at a lower temperature, 
and thus last longer and be available when there is no solar energy. 5B 

The best and most economical place to put the thermal mass is the 
floor, but a masonry floor that can't be covered by rugs, furniture 
and other things normally used in a living space would be unacceptable 
to most families. The area generally used is the walls, which is more 
expensive to install. There is also the problem of getting sufficient 
area that is directly in the sun's rays. Thus the walls will require 
extensive and expensive mass. The roof can be used for thermal storage, 
which would be very effective because heat raises. The obvious draw 
back is the expense of using heavy thermal mass on the roof. With the 
use of a phase change material the excessive weight could be somewhat 
eliminated. The following exhibit '3, shows the basics of direct gain 
passive solar heating. 

HOUSE AT CROOKED CREEK (SEVERE CLIMATE) 

A house constructed in Dubois, Wyoming, which has a 10 month heating 
season and is very cold, was designed to use passive solar energy. The 
house was designed as a long structure with three large clerestories 
facing south. The clerestories allowed the direct radiation to enter the 
rooms above eye level, thus reducing glare and filling the space with 
reflected light. The house is very well insulated and includes an air- 
lock entry. Reinforced concrete block forms most of the house, and along 
with an 8" concrete floor slab make up the thermal storage of the house. 
The insulation in the wood walls and roof is R-22, in the block walls 
(urea-f ormaldahyde foam sprayed in wall cavity) is R-30, and a 4" styro- 
foam board at the slab edge (R-18) is used. All the glass is thermopane. 
Roofing is corrugated metal. 5C 
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EXHIBIT 3 



The collection of solar energy is by the clerestories which had 
to be sized to insure sufficient energy could be collected. The glass 
area calculations were based on a rainimun exposure to a 15 degree 
angle of incidence (December 21 average, 45 degree north latitude). 

After heat loss calculations, the window calculations resulted in 100 
square feet for each living area and 85 square feet for each bedroom. 

Total insolation was calculated to be 2,300 BTU per square foot per 
day, of which 2,000 BTU per square foot per day (December 21) is direct 
radiation, plus 15% for clear day reflected radiation, and also some 
diffuse radiation. Transmission oD radiation through thermopane glass 
is 75%, which gives a total radiation entering the house of 1,725 
BTU/ sq. ft. /day. The windows were inset 12*’ to try and maintain 
a cushion of still air in front of the windows and thus minimize 
daytime heat loss. 5C 

The storage of energy is in the walls and floor. Direct radiation 
strikes the plastered block walls (sand finish to help insure diffuse 
reflection) and reflects the radiation to the floors. 25 to 30% of the 
direct radiation will be absorbed by the walls as determined by the 
thermal admittance and acceptance of the masonry masses. The acceptance 
and admittance formulas are by Mazria; Solar Energy Workshop Workbook , 
December 5, 1975. The masonry wall can admit, or store, 45-50 BTU/sq.ft./ 
hr. of the approximate total heat exposure of 200 BTU/sq. ft. /hr. 
before overheating at the surface. Which is why the light color for the 
wall surface was selected to allow only 25 to 30% of the radiation 
striking to be absorbed. The black slate finished floor accepts 90 to 
95% of the radiation striking it. This process of calculating and 
balancing color allows the maximun storage of energy, by minimizing 
surface overheating. This helps to insure energy for nighttime heating 
The mass in the building is calculated to store one days heat within 
a 10 degree fahrenheit temperature swing (65 - 75^F) . The amount of 
mass was determined from the thermal capacity of the walls and floor. 

All calculation assumed R-8 Technifoam shutters over the south clerestories 
at night. 5C 

The actual performance data was collected without the shutters. 

During clear days the temperature fluctuated about 15 degrees (60 - 75) 
with outside temperature ranging from 10 to 25^F daytime and -10, -O^F 
at night. There was a need for approximately 2 hours of nighttime 
heating with a ’*Franklin" woodstove. The house used one half cord of 
white pine from August to the end of November. There was only a 
maximun 5 degree temperature difference between the floor and the peak 
of interior spaces, because the mass radiates very slowly and gradually, 
causing the heat to migrate slowly. The house was monitored for a four 
day period of continuous snow. During this time outside temperatures 
ranges from -lO^F to -35^F, while the house never fell below 50'^F 
without any auxiliary heating. The reason this temperature can be 
maintained is there is enough diffuse radiation reflected into the house 
to maintain this temperature. This also shows what can be accomplished 
with a well insulated house. 5C 

When the walls and floor are warmer than the air space, as it is 
with a direct gain house, the house is comfortable at 65®F as compared 
to 75^F in a conventional house. Also around sundown when the house 
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changes from direct solar radiation heat to thermal energy (from 

the walls) the house will experience a chilly period of about 30 minutes. 

Once the walls begin to radiate, the house feels warm again. 5C 

HOUSE IN CENTRAL CALIFORNIA (MILD CLIMATE) 

A passive solar house that is located in a mild climate (central 
California) has been designed to give essentially 100% passive solar 
heating and cooling. Only the heating will be addressed in this project. 
The location has 589 degree F days of heating with a relatively 
dependable sunshine and large dirunal swings (28 degrees in winter and 
32 degrees in summer) . The house was built with materials and construction 
techniques standard to the area, and has 1,500 square feet of living 
area. This requires only 400 square feet of south facing glass and 
2,000 square feet of thermal mass (mostly 2’* thick) , which will allow 
the system to operate 24 hours a day. The temperature swing predicted 
for this house is 22 degrees F (64®F to 84‘^F) . This does not take 
into account movable insulation on all south facing windows, which 
would improve the performance. 5G 

The thermal storage is based on diurnal storage, which is optimal 
at about 2”. The floor is a 4” slab on grade, on concrete block wall 
capping, laid on a wood decking, filled with grout and sealed with wax, 
to give a light colored semi-reflected finished, so that the insolation 
will be spread over all the high mass surfaces. The thermal mass on 
the walls are 2" plaster where the sun strikes the walls from one side 
and 4” plaster where the sun strikes the surface from both sides. 5G 

Collection of solar energy is accomplished as shown in exhibit 4. 

This design shows that with the proper climate conditions a ^house can 
be heated and cooled almost 100% by passive solar energy, and be built 
using materials and construction techiques standard to the area, 
keeping costs down. An example of how costs compared to a conventional 
house, is the cost of the high mass wall which was 70c per square foot, 
as compared to 63c per square foot for a typical 1/2’* gypsum wall 
with 2 X 4 studs. There is no actual performance data available. 5G 

RETROFIT NEW MEXICO ADOBE 

The retrofitting of a direct gain system to a house will usually 
take considerable remodeling if the system is going to be effective 
in using solar energy. The house will also have to have a large amount 
of south facing wall. To install the needed south facing windows and 
storage mass in the house will probably require structural changes, 
which would be very expensive. I^There retrofitting can be effective is 
if the house already needs extensive remodeling, and is situated such 
that there will be sufficient south facing wall area. An example of 
this type of remodeling is an adobe house in Alcalde, New Mexico. 

This house needed extensive repairs. The owner decided to convert to 
solar energy (direct gain) with the remodeling. The house was ideally 
suited in that it was made of adobe, giving it mass, and had sufficient 
south facing wall. The solar retrofitting with remodeling cost 
$15,000, as compared to $14,000 that a contractor estimated for just 
the needed repairs. The estimated total heat loss, with a 15 degree F 
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average 24 hour temperature, is 185,000 BTU/day. The estimated solar 
gain through the 180 square feet of window is 180,000 BTU/day. Tlie 
180,000/day is based on 1,000 BTU/sq. ft. /day, which is a clear day 
value. The south windows have reflectors which will increase the 
energy gained (30% by Steve Baer’s rule of thumb), and should make 
the house very close to 100% passive solar heated, considering the 
area of New Mexico. 5K 
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THER>IAL STORAGE IN A WALL OR ROOF 



TROMBE WALL 

The Trombe Wall, named for Dr. Felix Trombe one of it’s developers, 
is a passive system that uses thermal storage on the south wall of 
the house. The exterior of the wall is painted a dark heat absorbing 
color, to insure as much energy is absorbed as possible. There will 
be either one or two layers of glazing mounted several inches from the 
wall. If it is desired to have both daytime and nighttime heating the 
wall will have to have vents at both the top and the bottom of the wall. 
These vents will have to have dampers to prevent the reverse flow of 
heat through the vent and out the windows at night. Exhibit 5 shows 
basically what a Trombe Wall is and how it works. 

As exhibit 5 indicates solar radiation (short wave) passes 
through the glass and is absorbed in the wall. The thermal energy, 
in the wall being at a low temperature, producing long wave radiation, 
will not pass easily through the glass, thus capturing the solar energy 
in the house. Exhibits 5 and 6 also show that the daytime heating 
is done by a natural convective loop, and nighttime heating is done 
by radiation from the thermal storage wall. 5B 

In the Trombe Wall, the selection of how much and what kind, of 
material to use for the wall is critical. Normally masonry is used for 
the material, but with the development of phase change materials 
(to be covered later) there may be a switch to this less space consuming 
material. Temperature swings, as experienced with direct gain systems, 
can be reduced, due to the smoothing achieved as the temperature wave 
diffuses through the wall. There is also a time delay between when 
the solar energy is absorbed and when the energy enters the house. 
Usually this time delay is from 6 to 12 hours, which gives the maximum 
heating when most needed, at nighttime. Exhibits 7 and 8 show the 
smoothing effect and the time delay of th Trombe Wall. 5B 

The results in exhibit 7 are more extreme than would be observed 
in a house, since test rooms have a large ratio of collector area to 
load (4.3 sq. ft. of collector per BTU/hr. ^F), and consequently the 
inside temperature averages about 50 degrees F above the outside 
temperature on sunny midwinter days. 5B 

The thickness of the wall can be some what flexible because of 
the time delay. The thickness that gives the maximum annual energy 
yield to the building is 12 inches. With a 12" wall there will be too 
much* of a temperature swing, which results in uncomfortable periods 
when there is an extended cloudy period. The usual design uses a 
16" wall, which will give more storage and less of a temperature swing. 
The exhibit 7 lists the characteristics of a solid concrete wall 
during sunny days with double glazing. 5B 

The Trombe Wall can be used a number of ways as a passive solar 
system. The Trombe Wall can be used by itself, but will need the vents 
for daytime heating. The Trombe Wall can be combined with an attached 
greenhouse, still using the vents. The Trombe Wall can also be used 
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THERMAL 
STORAGE 
MASS 
BTU/ “F 


THERMAL 
STORAGE 
SURFACE AREA/ 
GLAZING AREA 


INSIDE 

DAILY 

TEMP.^’F 

SWING 


TIME OF 
INSIDE 
TEMP. 
PEAK 


DIRECT GAIiN 


37 


2.8 


38 


3:00 P.M. 


16" TROMBE WALL 
(WITH vt:nts) 


32 


0.84 


26 


4:00 P.M. 


16" TROMBE WALL 
(NO VENTS) 


32 


0.84 


9 


10:00 P.M. 


WATER WALL 


35 


1.01 


25 


4:00 P.M. 



THICKNESS 


INSIDE SURFACE 


TIME DELAY 


INCHES 


TEMPERATURE 


OF PEAK ON 




SWING 


THE INSIDE 


8 


40 ”F 


6.8 hrs. 


12 


20°F 


9.3 hrs. 


16 


lO^’F 


11.9 hrs. 


20 


5<»F 


14.5 hrs. 


24 


2^F 


17. 1 hrs. 



EXHIBIT 7 
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EXHIBIT 8 



TEMPERATURES MEASURED IN A TWO-STORY THERMU STORAGE WALL. THE WALL 
IS MADE OF 12" HOLLOW CONCRETE BLOCK FILLED ^-TETH MORTAR. THE WALL IS 
DOUBLE-GLAZED AND HAS NO VENTS 
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just as a nighttime heating device, no vents, and some other form of 
heating during the day. How the Trombe Wall is used will be up to the 
designer. 

The problems x^;ith the Trombe Wall is the high cost of construction, 
it takes up valuable space within the house, and it may not have a 
pleasing appearance to some people. The problem of loss of heat at 
night, through the windows, can be reduced by using movable insulation 
between the wall and the glazing. Any overheating that takes place, 
can be resolved, but with proper design shouldn’t happen. 

DOUG KELBAUGH HOUSE 

The Doug Kelbaugh house uses the Trombe Wall with vents as its 
passive solar heating system. The Trombe Wall is constructed of 15 
inches of concrete^ which has been painted with 3M Nexel (special 
black coating) over masonry conditioner. There is double glazing in 
front of the wall, giving a total of 600 square feet of collection 
area, plus the attached greenhouse. The vents are located at the top 
and bottom of each floor level. After the first winters occupancy it 
was discovered that dampers were needed to stop the reverse thermo- 
siphoning. The dampers were made so that they would operate passively, 
and are made from a light cloth over a screen that would only allow 
flow in one direction. The house also has an attached greenhouse, 
which will be discussed at a later time. Since the house is two stories 
there was a problem with the upper story always being hotter than the 
lower story. This was resolved by closing off the lower from the 
upper level. Exhibit 9 is a schematic of the house’s heating system. 

The house has 2,100 square feet of floor area and a volume of 19,000 
cubic feet. The heat loss for the house is about 65,000 BTU/hour at 
O^F. The backup heating is a gas fired hot air furnace that can 
deliver 58,000 BTU/hour. There are also three 250 watt heaters in the 
bathroom to provide instantaneous heat. The additional cost estimated 
to be attributable to the solar heating and cooling system is $8,000.00 
to $10,000.00. This cost for a solar system may seem high, but the 
following exhibit 10, shows why it is worth the extra money. 5DD 

The house had a temperature swing of 5 to lO^F over a 24 hour period. 
Daytime temperature setting was between 60 and 64^F and the nighttime 
setting was 58^F. These settings were acceptable because the radiant 
heat from the wall allows lower room temperatures to be comfortable. 

The estimated downstairs average was about 63^F while for upstairs it 
was 67‘^F. 5DD 

Since this system obtains a high percentage of total heating from 
solar it makes the high cost acceptable. When maintenance costs, 
reliability and convenience are also considered the system is even more 
and more appealing.. 

The Trombe Wall concept can be used with no vents directly into 
the living space, but with a rock bed storage system in which the air 
between the Trombe Wall and glazing is passed through. The Bruce Hunn 
house is an example of this kind of system, but since it is considered 
a hybrid system it will be covered later. 
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1976 - 1977 Winter Perf prmance 
Heat Loss; 5,556 Degree Day; Design Temperature 
(65 Inside, 0®F Outside) 
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RETROFITS 



WICKEL HOUSE 

Retrofitting a Trombe Wall on an existing house is not too costly 
if the house already has a south facing masonry wall. The system will 
probably not be as effective as a new construction system, but it 
should reduce the heating bill. Two houses that were retrofitted with 
Trombe Walls are the Wickel House in New Jersey, and the Upper Black 
Eddy House in Pennsylvania. The Wickel House is a converted vacation 
lodge with over 7,000 square feet of living space, with three stories. 

One gable end is exposed to 5 degrees west of true south. Since the 
house was originally designed as a lodge, it was possible to limit 
the heating load for the solar system to 2,200 square feet of living 
space in five separate zones. The total heating load for this portion 
of the house was estimated to be 205 million BTU per heating season. 

The south end of the basement and first floor (432 sq, ft, ) was 
constructed of hollow core cement blocks, which was ideal for retrofitting 
a Trombe Wall. 5D 

The first step in constructing a Trombe Wall was the preparation 
of the surface of the wall. The surface already had a rough coat 
plaster job, which was nice for solar absorbance and creation of air 
turbulance for efficient heat transfer. The wall was blackened with 
a thin slurry of mortar mix, with black masonry pigment added. The 
air space for the Trombe Wall was formed with 2 X 4 ^s^ giving a three 
inch air space. The perimeter was then caulked inside and out with 
caulking^and lined with one inch fiberglass drop-ceiling insulation. 

Lead anchor shields were placed in the masonry at 30 inch centers 
across the wall. The glazing was attached ds indicated in exhibit 11. 
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distribution. It can also compensate for irregularities in the wall 
by varying the pipe spaces length. All the forming members were painted 
white to reduce the amount of heat they absorbed. Three 6" X 12** 
vent holes were cut at the top of the collector to allow air flow 
into the third and fourth floor areas. 5D 

The system operates by convective air being fed through the 
basement windows, then by opening or closing selective windows and 
vents, the heat can be directed to any of five living areas as desired. 
Dampers are used to prevent reverse thermosiphoning. The 8’* hollow 
block doesn’t provide very much thermal mass, which is why the 
insulation and sheet rock were left on the south wall. By leaving 
the south wall the way it was helped reduce the overheating problems. 5D 

The air entering the Trombe Wall via basement windows was 38 
degrees, with the exhausted air ranging from lOO^F to 122^F. Since no 
air flow measurements were taken, these temperatures have very little 
meaning. During winter days, use of the wood stove has been virtually 
eliminated. The limited thermal mass gives very little storage, which 
means when the sun goes do\^m the system is of very little use. The 
fact that the area has nearly perpetual snow cover and clear skies 
improves the systems performance. The Wickel family has been pleased 
with its performance. The best feature of the retrofit is the low 
cost. The total expenditure for materials was $378.55, resulting 
in 88c per square foot of collector area. 5D 

UPPER BLACK EDDY HOUSE 

The Upper Black Eddy House is a two story stone farmhouse with 
1,000 square feet of living space. The heating load is 110 million 
BTU per heating season. The primary source of heat is a small log- 
burner wood stove. The house has one side that is 5 degrees east to 
true south. After deletion of obstructions there was 28’ by 5 ’ of 
stone wall exposed to the sun, of which one third is window space. 

To make full use of the limited space available for Trombe Wall it 
was decided to use double glazing. The construction of the Trombe 
Wall was similar to the Wickel House but had to be modified to accept 
double glazing^in lieu of single Kal Wall. The builder, taking 
advantage of good deals^got the materials, for the project, for $113.00, 
or 81c per square foot of collector. The materials used caused the 
labor time to increase considerable over the Wickel House. The house 
has a 15 degree pent roof sloping away, for 10 feet, from the base of 
the collector, which will improve its efficiency. 5D 

The system operates by air being fed into the collector via the 
lower sash of the three windows and exhausted through vents cut into 
the eves, where it ties into an existing unused insulated hot air duct 
system. The hot air is then fed into the second floor near the ceiling. 

The second floor will quickly reach 60 to 65^F with moderate insolation, 
but as inlet air temperature increases the heat transfer becomes 
inefficient, thus slowing down the thermosiphoning. There are plans to 
install a circulating fan to help direct the air through out the house 
for a more even heat distribution, and to help maintain the thermosiphoning 
loop. This system will not meet a major percentage of the overall 
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heating needs, but it will be significant enough to make the investment 
worthwhile. No actual performance data was available. 5D 

Retrofitting a Trombe Wall to an existing house will not give 
the same kind of performance as a new construction system, but as 
shown with the Wickel House and the Upper Black Eddy House satis- 
factory performance can be obtained at a minimun of cost. This is 
important because it allows homeov-mers in the lower income strata 
to take advantage of solar energy. It should be pointed out that 
the above only applies to houses that have a masonry wall that faces 
very close to the south. 



WATER WALL 



The water wall system works on the same principals as the Trombe 
Wall, except water is used as the storage medium instead of masonry. 

The main advantage of water is its efficiency in picking up heat, 
storing it, and releasing it, resulting in less storage space. The 
surface of the water will also remain cooler than the masonry. The 
cooler surface will result in less heat losses and will have less 
tendency ot overheat during sunshine hours. The lower surface 
temperature will also be more responsive to solar gain. 5Y The water 
wall can be directly exposed to the living space or it can be enclosed. 
If the water wall is enclosed the air can be circulated naturally or 
forced. Tests conducted by the Energy Center of the University of 
California, San Diego indicate that the water wall bin with forced 
convection gives the best solar load fraction, which is dependent on 
a person’s subjective assessment of the comfort enhacement by direct 
radiation. The water wall bin with forced convection does provide 
the best room temperature control. During the daytime temperatures 
rarely exceed the thermostat setting by more than 2^F. The placement 
of an insulating wall (even with vents) between the living space and 
the water, with natural convection, will give poorer performance than 
the exposed wall. There was a significant loss of energy through 
the windows by conduction, even with R-8 night insulation. This heat 
loss can be reduced from 50% to 100% by the use of night insulation 
on the exterior of the glazing. Exhibit 12 is a schematic of a water 
wall bin. 5Z 

The primary problem with the water wall is getting the water in 
the vertical configuration. Initial attempts had continous water from 
the top to the bottom of the wall, such as vinyl bladders inside of 
walls and various types of tubes made of steel and fiberglass. These 
methods put too much head pressure on the system, making leaks a 
major problem. Corrosion is also a problem when metal containers are 
used. Drums are usually used, in that the water pressure can be kept 
down and still form a wall by stacking the drums. If metal drums are 
used to contain the water the corrosion can be minimized by using 
anticorrosive stabilizers in the water, or by lining the containers 
with something that is not corrosive in water. If the system is an 
exposed water wall there will have to be something to shade the wall 
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during the sumiter when no heating is desired. With the water wall 
bin the system can be designed so that the solar heat is vented 
outside, it would still be a good idea to provide some kind of 
shading. 5Z 

FOCUSING ROOF APPERATURE /HOUSE IN NEW JERSEY 

A technique used to minimize the thermal storage area is the 
focusing foof apperture. The roof aperture will allow the storage 
on the north wall. The designing of the aperture to give optimun 
distribution of solar energy is very complex, and will not be covered 
by this project. There is on going study, by others, in this area. 

The aperture will not focus the solar radiation on to the water storage 
at all times, therfore, a concrete floor should be used to insure all 
the focused energy is absorbed. The technique of focusing can 
minimize the storage area by a factor of up to two to one. This will 
save money as storage costs about $2.00 per square foot. The focused 
radiation will require a high capacity material like contained water 
or phase change material to avoid overheating. A masonry wall however 
thick, (8^* to 16”) will not provide enough storage to prevent over- 
heating. The contained water should be situated on the north wall 
such that it will receive the most radiation. An advantage of the 
focused aperture is that during the summer months when the sun is 
higher, solar radiation will not be collected. 51 

A house located in New Jersey has a passive system that uses both 
a water wall and a focused roof aperture. The location has 4,980^ 
days per year with 57% yearly sunshine. The house has 1,400 sq. ft. 
of living space and a heat loss of 8.3 BTU/DDF/sq. ft. which is a 
well insulated house. Exhibit 13 is a schematic of the house’s 
heating system. The water wall operates by absorbing the solar energy 
through the glazing, and then radiating it from the v/ater storage 
through the inner wall. The system also uses natural convection, 
from the water tube space, through the vents into the living space. 

The vents have dampers to prevent reverse thermosiphoning. There 
are exterior reflectors to increase the solar energy collected, which 
can be pulled up at night to give nighttime insulation. 

The focusing roof aperture operates by passing the solar radiation 
through the outer glazing material and reflecting it, by one or more of 
three reflectors, to the water storage or concrete floor. The three 
reflectors are indicated in exhibit 13 as (T) (^, and(^. There is 
internal glazing at the ceiling to prevent convective coupling of 
the skylight with the occupied space. Reflector two is movable so 
that it can also act as insulation. The water storage and concrete^ 
transfers the heat to the house by natural convection and radiation. 
This system has a good heat distribution because there is a heat 
source at both the south and north walls. The above system has a ratio 
of collector area to building load of about 1.7 sq. f t . /BTU/hr . ^F. 51 

The simulation used for this house had many simplifications. 

The validity and effect of these assputions are as of yet unknown. 
Future simulations are going to be done by computer, which should 
give more accurate results. The results of the simplified simulation 
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are as follows: 



YEARS TOTALS 

HEATING HEATING LOAD PRECENT SOLAR LOAD 

DEGREE DAYS (BTUXIO^) SraSHINE FRACTION 

4,980 58.42 57 90% 

Total cost of combined water wall and roof aperture system was 
$7,858.00 which includes materials and labor. Considering costs, the 
solar load fraction, possible tax incentive and the low maintenance 
requirements makes this system very attractive. The above solar load 
fraction does not include the energy needed for hot water. The hot 
water required 3.04X10^BTU per year, which is 25% of the total annual 
requirements. The system could be modified to use the water storage 
as a pre heat tank for the hot water, which would decrease the 
conventional energy required for hot water. How this would effect 
the overall performance would have to be simulated. 51 

RETROFITS 

Retrofitting a water wall or water storage system should be 
effective. The construction of the water wall or water storage should 
cost about the same as a new construction system. With a water wall 
it may take some structural work to get a wall with a lot of glass 
area, but it shouldn’t be too much of a problem. The house will have 
to have a suitable wall facing south. If the system is properly 
designed it should be as effective as a new construction system, 
unlike a retrofitted Trombe Wall. The house will have to be well 
insulated, but a house should be well insulated no matter what kind of 
energy is used for heating, 

WATER ROOF 

The water roof system is a concept used in areas where flat roofs 
are acceptable. Most of the work on this system has been in the south- 
western part of the United States. The system uses a water storage 
pond on the roof of the house. The solar energy is absorbed during 
sunlight hours, and a movable insulation is used to keep from losing 
energy to the air at night. There are movable insulation systems 
that can be moved with nonrenewable energy, and others that can be moved 
with manpower. Exhibit 14 is a schematic of the water roof system. 

As can be seen from exhibit 14 , the heating is done by radiation 
only. The system radiates heat into each room, making the heat 
distribution no problem. With the heat source in the ceiling it may 
cause a problem of the floor being uncomfortable, in the systems 
built to date this has not been a problem. The water roof system 
should be able to handle 100% of the heating load in areas where it 
is practicable, and also 100% of the cooling load, if the weather 
conditions are right. The construction of a house capable of carrying 
the load of a water pond on the roof has not been a problem, but 
have had trouble keeping the system from leaking. There is ongoing 
work to try and improve the systems water tight integrity. A family 
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that lived in a completed house felt that the system was superior to 
conventional heating and cooling. Of particular note was the lack of 
blowing air and noise. 4 

The water roof wystem doesn*t lend itself to retrofits, since 
there would have to be extensive structural work done to be able to 
take the added load of the water. The cost of the structural changes 
would make the retrofit prohibitive. 
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SOLAR GREENHOUSE 



The solar greenhouse is a room that has extensive glass area that 
faces south. It is really a direct gain space that is thermally 
separated from the living space. This allows the greenhouse to have 
large temperature swings, while the living space, being buffered, will 
have smaller temperature swings. There will be a theraal storage 
wall separating the greenhouse from the living space. The wall can 
either be insulated on the living space side or left uninsulated. If 
the wall is left uninsulated it will operate on the same pricnipals 
as a Trombe Wall. If the wall is insulated, to prevent heat loss when 
the greenhouse is cooler than the living space, the heating is done 
by convention through vents, windows or doors adjoining the greenhouse. 
The periraiter of the greenouse will be insulated along with the floor 
if the climate is severe enough. Additional thermal storage, such 
as water filled drums, may be used in the greenhouse to improve its 
performance. Exhibit 15 shows the basics of greenhouse heating. 5B 

The greenhouse with its extensive glass area is a good collector 
of solar energy, but the glass is not a very good insulator to thermal 
conduction. Shutters can be used to increase the greenhouse resistance 
to conduction. There has been recent development of transparent 
insulation, to be covered latter, that could make the greenhouse a 
better thermal storage medium. The greenhouse also can be used as a 
plant growing area, which could be used to grow vegetables year around. 
If the greenhouse is open to the living space there could be a problem 
with insects, if proper care is not taken. 

DOUBLE BOX 

A unique greenhouse system that is being developed is the 
double box. The double box is a greenhouse that completely surrounds 
the living space. Exhibits 16 and 17 are schematics of the double 
box. The structure will have to be well insulated. The air is heated 
on the south side through the glass raising to the top. The air is 
then forced down the north side through an air plenum space into a mass 
storage area at the base of the greenhouse. The air will then proceed 
into the greenhouse and up through a vent to be heated again. The air 
will circulated at about twice that of a normal central conditioning 
system, which will give air exchanges in the greenhouse of about 
one per minute. The reason for these high air speeds are; 5AA 

1. Lower delta temperatures of operation, reducing heat loss 
back through the glass. 

2. Rapid Air FLow 

a. Rapidly looping air prevents heat from building up in 
the greenhouse 

3. Increased storage mass , 

a. Storage at low temperatures increase operation time of 

heat supply. 

The double box gives two modes of operation; the air loop and the 
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air buffer, which provide 3 main temperature zones; the inner box, 
the outer box and the outside air. The outer box^s temperature will 
have big temperature swings, while the inner box will be maintained 
at a fairly constant temperature. During the heating season when the 
sun is shining the outer box will be warmer than the inner box, thus 
the inner box will retain its heat. When the sun is not shining the 
outer box will be warmer than the outside air, due to thermal storage, 
and will funtion as a warm air blanket thus reducing the heat loss 
of the inner box. The energy used to move the air in the outer box 
can be achieved by either an electric fan, 10 to 15 KWH per month, or 
by setting up a natural convective loop. The natural convective loop 
can be obtained by leaving a portion of the north air plenum 
uninsulated. The warm air will escape out the north wall and be 
replaced by the warmer air moving down from the attic, thus forming 
the convective loop. 5AA 

A demonstration project is being built, that uses the double box 
concept. The project, when fully operational, is expected to have a 
potential gain that is more than twice the heat losses. There will 
be an excess heat storage medium (low delta T) that should have the 
capacity (2 million BTU) to last through the cloudy days. The double 
box configuation can also provide summer cooling. Since this project 
is not completed there is no actual data available. There is also no 
figures for the cost of construction. This system appears to have 
potential, but how much will have to be assessed when actual data 
becomes available . 



DOUG KELBOUGH* HOUSE 

The usual application of greenhouses is in combination with other 
systems. The Doug Kelbaugh House is primarily a Trombe Wall System, 
but has an attached greenhouse for additional heating and growing 
ornamental and edible plants. Initially, the greenhouse was single 
glazed and opened directly into the living space. After the first 
winter’s operation it was found that the greenhouse was losing an 
excessive amount of heat, about 40% of the total heat loss of the 
house. To help remedy the heat loss problem a second layer of glass 
was added, which was estimated to cut the heat loss from the green- 
house in half. The greenhouse also had an excessive temperature 
Swing, up to 25^^?. Eight 55 gallon drums painted black where added 
to try and reduce the temperature swing and to store energy. The 
Doug Kelbaugh House showed that if a greenhouse is going to be open 
to the living spaces it should have some means of uncoupling the 
greenhouse at night. Shutters can also be used to reduce heat losses 
and give privacy. Therefore, if a greenhouse is going to be used 
in combination with another system and is going to be open to the 
living space it should only be considered as a daytime heating 
source. 4 



HYBRID SYSTEM 

Another interesting combination is to use the greenhouse attached 
to a Trombe Wall and use the heated air in the greenhouse to charge 
a rock bed storage system. This system will have the wall and the 
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floor giving off radiant heat, which will allow the room air temperature 
to be lower and still be comfortable. Exhibit 18 is a schematic of the 
above system. The concrete floor will provide a time delay for the 
heat flow. The warm floor will also offset the tendency of the room 
to stratify in temperature. 

RETROFITS 

Perhaps the best use of the greenhouse is in retrofitting it 
to existing houses. The greenhouse can be an effective solar system 
if the house it is retrofitted to is well insulated and has sufficient 
south facing wall area. It will not take care of the major heating 
needs, but will reduce it. The cost of construction will be similar 
to a new construction system. 
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CONVECTIVE LOOP (TliERIlOSIPHON LOOP) 



The natural convective loop operates similar to active air 
systems, except no conventional energy is needed. The collector will 
be similar to the active system. The collector has one or two layers 
of glazing, a heat absorbing plate and a method for the heat transfer 
fluid to flow across the plate. As the fluid passes across the plate 
it will absorb the thermal energy in the plate and transfer that 
energy to the storage medium. As the fluid is heated up it will 
raise and go into the storage bin at the top of the collector. The 
storage medium will absorb the thermal energy, cooling the fluid ^ 
causing it to flow do^^ and into the bottom of the collector. This 
system can be used with either water or air as the heat transfer 
fluid. 5B The usual use for the above system is for passive hot 
water heaters. 

The natural convective loop air system can be used as shown in 
exhibit 19 . The house will have to be well insulated as usual. 

The collector will heat up the air causing it to raise into the 
storage bin. The rock bed storage will give up its heat by allowing 
air to raise through vents into the living space. Once the air is in 
the living space it will raise, hearing the second floor spaces. As 
the air cools, it will return through the cold air return ducts. 

The hot water can also be heated in this system by installing water 
pipes in the collector. The system has dampers on the vents to control 
the heating and prevent any reverse thermosiphoning throughout the 
system. A house has been built using the above concept, but no 
performance or cost data was available. 3 

For the natural convective loop system to be justified it will 
have to have better performance than an active system, and be com- 
parable to other passive systems. The reason for this is that the 
system requires a collector that has no other function than to collect 
energy, as in the active system. The other passive systems have good 
performance and serve more than one function. There is a use for 
the convective loop system in retrofitting existing houses. The 
retrofits can be small convective air heaters for use during sunlight 
hours only, or a simple window box heater if it can be designed not 
to lose heat at night. Work is on going in this area of small 
convective air heaters, which can play an important role in using 
solar energy, espically in existing houses. 5B 



page 32 



EXHIBIT 19 



Loop 

( /) 



J. N^u fA~tf c>aJ 





Coc>leJ Aih- ifJp 

bee/ 



ouU^ ql/AZ/FO^ 

L/Ay^reJ Mesk 

Sov-bcf- 
piAi^ 




'^'^9 

ooOO \ ^b \ 6'^ \d 



n 



<ir 



] [ 






<r 



page 33 



MODIFICATIONS 



VERTICAL SOLAR LOUVERS 

Modifications can be made to the basic passive solar concepts 
to make better use of collected energy. An interesting modification 
is the vertical solar louvers, which combines the concepts of direct 
gain and thermal storage wall. Vertical solar louvers are a set of 
rectangular columns situated directly behind south facing glazing 
and oriented in the SE-NW direction. By orienting the columns in 
this fashion the morning sun is allowed to enter and warm the house, 
and the afternoon sun will store energy in the columns. Exhibit 20 
is a schematic of the vertical solar louver system. The vertical 
solar louvers has most of the advantages of the direct gain and the 
Trombe Wall systems, while eliminating most of the disadvantages. 

The Trombe Wall will give even temperatures, but does not allow the 
use of alot of possible storage in the interior floor and walls. The 
Trombe Wall, that circulates room air, also has no convient means 
to clean the interior glazing, which can deteriorate the system 
preformance. The vertical solar louvers allow the use of thermal 
storage in the walls, and still have easy access to the glazing for 
cleaning. The primary disadvantage of the direct gain system is the 
wide range of temperature fluctuations, which the vertical solar 
louvers system will not have. 5J 

The vertical solar louvers are designed to function so that in 
the morning when the house is at its coolest the maximum amount 
of solar energy is directly admitted into the house, and then around 
noon switch to a Trombe Wall for storage of solar energy and to 
prevent the house from overheating. The louvers will admit over 
80% of the morning sun directly into the interior. As the day 
progresses less and less of the sun’s direct energy x^ill be allowed 
into the house. At noon from 40 to 70% (depending on column spacing) 
of the sun energy will be intercepted by the louvers. In the after- 
noon all the sun’s energy will be intercepted by the louvers for 
storage, after the house’s interior has been heated up. The vertical 
solar louver system should reach maximum comfortable temperature 
at about 11:00 A.M. and last until 5:00 P.M., as compared to 1:00 P.M. 
to 5:00 P.M. for other hybrid passive systems. The maximum surface 
temperature for the louvers will be reached between 4:00 P.M. and 
5:00 P.M. , as compared to shortly after noon for the Trombe Wall. 

With the maximum surface temperature being obtained around 4:00 P.M. , 
it will give a heat pulse at the back of the louver between 2:00 A.M. 
and 5:00 A.M. (depending on the column material and thickness) when 
it is most needed. The vertical solar louver, therefore allows the 
house structure to absorb and store heat for a longer time at the 
maximum comfortable temperature. The vertical solar louvers should 
improve the storage effectiveness of the internal structural elements 
from 50 to 100% in December and January. 5J 

An advantage of the vertical solar louvers over the Trombe Wall 
is the front surface temperature of the vertical solar louver is 
lower. The vertical solar louver^ s maximum surface temperature is 
close to 75 F as opposed to 125 F for the Trombe Wall. The lower 



page 34 




\/^V^\6aI -S'o Iav- L Q ^ 



EXHIBIT 20 



page 35 






surface temperature will give a lower temperature potential for heat 
loss to the exterior, thus resulting in less heat loss (up to 50%). 

This means that less expensive night insulation can be used and still 
maintain good performance. The actual storage mass temperature of the 
vertical solar louvers should be about the same as the Trombe Wall. 

The problem of excessive floor temperature, that can occur with a 
direct gain system, is also avoided. 5J 

Jim Bier, the designer of the vertical solar louvers, has built 
a house using his concept. The house has 1,300 square feet of living 
space, has two stories^ and is made of 8 inch concrete blocks tamped 
full of sand. The floor is 4 inches of reinforced concrete on a 
4 inch open gravel base. The house is well insulated, even to using 
triple glazed storm windows, but no overall heat loss coefficient 
was provided. The house has 300 square feet of south wall, 50 square 
feet of clerestory and 70 square feet total in east and west walls, 
which are used to collect the sun’s radiation. There are 100 tons of 
insulated masonry and sand in the house to provide storage for overnight 
and on partially cloudy days. The five vertical solar louvers are 
12” X 32” sandfilled blocks, the full height of the south wall 
angled 45 degrees east of south. The south wall is double-glazed. 
Insulating drapes, stored on the NE interior side of each louver, are 
used to prevent nighttime heat loss and to control unwanted insolation. 
The system was designed to give from 50 to 85% of the houses needs by 
passive solar energy. There is no actual performance data available. 

The cost of constructing a house, using the vertical solar louvers, 
should be about the same as a conventional block house, but no com- 
parative cost data was provided, 5J 

Combining and modifing the basic concepts of passive solar energy 
is often done in efforts to maximize the use of solar energy. The 
greenhouse is often combined with water walls of Trombe Walls, or in 
the Doug Kelbaugh House. Jim Bier’s vertical solar louver is a 
creative mofification that makes use of the advantages of the direct 
gain and Trombe Wall systems, and eliminates the disadvantages. 
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HYBRIDS 



Perhaps the biggest source of disagreement with the people involved 
in solar energy is the definition of when a system is passive or 
active. There are some people who contend that if any source of 
conventional energy is used the system must be classified as active, 
and others who try to use a coefficient of performance figure to 
distinguish between the two types of systems. For this paper there 
is a catergory of passive systems, hybrids, that use conventional 
energy with one of the basic passive solar concepts. The conventional 
energy will usually be a fan to force warmed air through a storage 
medium. The reason hybrid systems are used is to overcome the over- 
heating and underheating problems of passive systems. Examples of 
hybrid systems will be^ given to show how the concepts of passive solar 
heating can be combined with conventional energy to produce an effective 
solar heated house. 



UNIT 1 FIRST VILLAGE 

Unit ill First Village, in the Santa Fe, New Mexico area, is a 
2,300 sq. ft., two story house that uses the greenhouse concept, 
with a rock bed storage to do the space heating. All rooms face 
into a t rianglar-shapedj 20 foot high greenhouse, on the south side 
of the house. The side of the greenhouse facing south is made of 
glass, while the remaining two sides are made of adobe, for thermal 
storage and to separate the heat collecting greenhouse from the living 
spaces behind them. There is about 400 sq. ft. of south facing 
thermopane glazing mounted at a 60 degree angle. The adobe is 14” 
thick at ground level and 10” at the upper level, to provide the passive 
thermal storage. The remainder of the house is well insulated. The 
rear wall is rounded to provide less resistance to harsh winter winds. 4 

Exhibit 21 is a schematic of the house’s hybrid passive solar 
system. As can be seen from the schematic, the system is made up 
of a passive part and an active part. The passive part uses the green- 
house principal to store energy in the adobe walls. The heat stored 
in the adobe walls will work its way through the adobe and be avail- 
able at night when needed. The unshuttered greenhouse in this system 
also acts as a buffer zone for the living space, in that it is allowed 
to fluctuate, usually between 110 F to 45"^F. The thermal wall will 
average about 73^ F. The active part is made up of two horizontal 
rock beds located under the living room and dining room. The beds 
are two feet deep and 10 feet wide, with one being 19 feet long and 
the other one is 15 feet long. The rock bed has sufficient storage 
capacity to carry the house through a couple of sunless days. The 
hot air is pulled from the top of the greenhouse, by two 1/3 hp. 
fans, and forced through the tw^o rock beds. The air is then returned 
to the greenhouse. The heat in the rock beds will radiate through 
the 6 inch concrete floor and into the living spaces. The temperature 
of the floor along with the rock beds will range from 85^F to about 
78^F. The rock bod will go down to 70 after a sunless day. Backup 
heating is provided by baseboard electric heaters. The domestic hot 
water heating is provided by a two panel flat plate collector array. 4 
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Temperature data collected over a four dav period, lowest 
outside temperatures normally experienced in the Santa Fe area, 
illustrated the thermal stability of the house. The temperatures 
ranged from about 57^F to 63^F for this period. During most of the 
winter the internal living spaces where able to maintain temperatures 
in the upper 60* s. There has been no monitoring of conventional 
energy useage to try and determine what percentage of the heating load 
was supplied by solar energy; but utility bills have consistently 
been less than $10.00 per month, usually being the minimum service 
charges. This was true even under the most adverse conditions. There 
was initially an apprehension that the sloping greenhouse would cause 
excessive heating in the summer, this has not been experienced. The 
owners have been very pleased with the house’s performance. 4 

BRUCE HUNN HOUSE 

The Bruce Hunn house has 1,955 square feet of living space and 
is located in Los Alamos, New Mexico, which has 7,000 F days with 
plentiful sunshine. The house has a two story Trombe Wall with 300 
square feet of effective area. The system does not use a natural 
convection loop to circulate the air into the living space, but a 
blower to circulate the air between the glazing and the wall into a 
rock bed storage. A three-zone, forced-air distribution system, 
with a natural gas auxiliary furnace is also connected to the rock 
bed. There is also 140 square feet of windows for direct gain solar 
collection. The house has a heat loss coefficient of 8.5 BTU/^F/day/ 
sq. ft. which is a well insulated house. There are plans to install 
interior insulating shutters on most windows, which will reduce the 
heat loss coefficient to 7.5 BTU/^F/sq . f t . The window shutters should 
improve the systems overall performance. ' 

The Trombe Wall is constructed of 12" open slump block completely 
filled with cement. This is considered the optimal thickness for 
concrete. The outside surface was painted a dark brown stain that 
has a measured solar absorptivity of 0.91. The wall was covered 
with double glazing mounted 2 inches from the Trombe Wall. Since 
the wall has no vents for a natural convection loop, the only space 
heating is by radiation. In the active mode the hot air in the space 
between the Trombe Wall and the glazing is forced through the rock 
beds. The blower initially circulated the air at a rate of 500 
cubic feet per second, which as will be discussed later, is too great. 

A differential controller, which operates on the temperature diff- 
erence between the top of the wall and the cold end of the rock bed, 
controls the charging process. The rock bed has a 3 foot by 12 foot 
face, with 6 feet in the direction of flow, and is 12 tons of gravel. 
The rock bed is insulated with extruded stryofoam, 2 inches on the 
side and 3 inches on the top, and rests on a concrete footing. The 
three zone forced air distribution system, with a natural gas auxiliary 
furnace, is connected to the rock bed. When a zone needs heat, 
returned air is ducted to the cold side of the rock bed too the warm 
end, and finally through the furnace. The furnace is set so that 
when the inlet air temperature is below 80 to 85^ F it will come 
on. 4 
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For first winter* s operation the system was in the passive mode 
only. The passive system under most conditions was capable of main- 
taining the inside temperature at 60^ F or above. As long as the Trombe 
Wall was charged this temperature was satisfactory to occupants. It 
was found that when the wall was not charged, during cloudy periods, 
the house was uncomfortable unless the room temperature was 65“" F 
or greated. Temperature readings within the Trombe Wall indicated the 
outer surface would reach a maximum of 152"' F around noon. This peak 
temperature would than move through the wall , reaching the inner surface 
from 6 to 8 hours later and be about 90°F. Temperature data taken 
oven a years period shows that the Trombe Wall gives a very stable 
inside temperature*. The wall performance, being only 12 inches 
thick, will suffer when there is prolonged cloudy periods, which 
temperature data taken has verfied. 4 

The active mode was turned on November and December, of the 
test year, to try and activate the rock bed. It was found that the 
inlet air temperatures to the rock bed were too low to charge the 
bed adequately. The active mode was therefore turned off. The 
reasons suspected for the poor performance of the active system 
are as follows; 

1. cold air leakage into the s-ystem 

2. direct losses 

3. too high an air flow rate 

Steps will be taken to remedy these problems. The blower was slowed 
do\^m to 200 - 250 cfra, because the design value of 2 cfm per square 
foot of collector area was based on the optimun value for an active 
collector, and is too high for a passive collector. Since the wall 
v/ill absorb only one-half to two thirds of the incident solar energy, 
a value of 1 cfm per square foot of collector was considered more 
appropriate. As mCny of the air leaks as possible will be sealed. 

The air flow across the Trombe Wall will be reversed to shorten 
the length of warm air ducts. 5K 

The system's first year performance was monitored to determine 
the percentage of the total heating load supplied by solar energy. 

The system was capable of providing 60% of the total load, without 
the active system operating properly. When the shutters are put on 
the windows^ and the active mode is working ^the system should be 
able to handle a good percentage more of the total heating load. 

The cost of material and labor attributable to the hybrid passive 
solar system is about $5,500. 00^ which is less than 10% of the total 
construction cost. 5K 

MARKED RESIDENCE 

The Markle House in Vermont has 1,100 square feet of living area 
with a heat loss cf 17,500 BTU/hr. Particular attention was paid to 
reducing the inflitration rate during design and construction. The 
windows on the east^west and north sides are triple glazed. Since the 
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oimers wanted a south view, the Trombe Wall concept was not feasible. 

The main living area is an open double height area with glazing placed 
at conventional height^ and up near the top to direct the light to the 
back of the room. This was done to try and alleviate the glare problems 
associated with direct gain systems. The hybrid system is a direct 
gain system with a forced air rock bed storage. The auxiliary heat 
is supplied by a wood fueled stove. 5M 

Exhibit 22 is a schematic of the houses^ heating system. The 
house faces south 17.5 degrees east. Solar energy comes through two 
banks of night shuttered glass, glass doors main level and clerestory 
glazing at the top. The solar radiation is absorbed by a dark slate- 
type floor^and by a full height south facing wall of dark green marbel 
rubble. These areas act as the direct gain portion of the system, 
and to heat up the air in the living room. The hot air will raise to 
the top, and when hot enough will be forced into the rock storage by 
a blower. The rock storage is behind the marble wall^and is 16* long 
by 6' wide by 12* high filled v/ith three inch round stone. If the 
room requires heat from storage, and a differential thermostat indicates 
that the storage is warmer than the room, three dampers activate 
and the cool air pooling in the living room is drawn up through 
the storage and circulated throughout the house via conventional warm 
air ducts and registers. This type of system allows the solar energy 
collected to be distributed through the house, even to areas that are 
not in close proximity to the radiating panels. 5M 

How the wood burner stove is used in this system is very unique. 

The primary problem with the way conventional wood burner stoves are 
used is that they are essentially air tight, so that they can burn 
for a long time without being reloaded, l^/hen the stove is operated in 
this manner it burns at a lower temperature. This type of operation 
will cause creosote to form in the stove pipe, because the cool chimney 
will condense the steam from the wood and combine with the volatiles, 
which is a problem. To overcome this problem the wood burner stove 
was placed in an opening in the marble wall, backed right up against 
the rock storage bin. The wood burner stove is then fired so that high 
enough temperatures are maintained in the stove pipe to eliminate 
the creosote problem. Some of the heat from the stove will go into 
the living space, but most of it will radiate into the rock bin storage 
and the marble wall. The stovepipes, with heat exchangers, also will 
radiate to these spaces. The hot air created by the stove will raise 
and be put into storage. 5M 

The performance of this house has been difficult to determine 
because it was not used full time, it is being used more like a vacation 
house. With the house being used in this manner the shutters are 
not operated properly, when not occupied, resulting in the rock storage 
system not being fully charged. The small amount of data gathered 
does indicate that the system is operating as designed. The stove was 
designed to carry up to 46% of the January load. There was no data 
available on the systems cost. The designer of the house indicated 
that a transparent insulation, to be covered latter, could be used in 
place of the mechanical shutters if the house was not going to be 
occupied full time. The designer also has plans to retrofit the house 
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with an automated signal to indicate when to start the wood burner 
stove. 5M 



SOLAR ATRIUM 

The solar atrium is a room which is surrounded on the north, east 
and west by other rooms of the house. The cost effectiveness of the 
solar atrium concept is the primary concern of its developers. A 
test model residence, using the solar atrium concept, is going to be 
built. The house will have two stories with 1,800 square feet, and 
will be a well insulated wood frame house \>r±th a peak heating load 
of 37,800 BTU/hr. 5N 

Exhibit 23 is a schematic of the solar atrium concept. The south 
facing wall has 600 square feet of glass in the wall and roof. The 
glass is 1" insulation glass, with inside louvers which have a total 
heat loss coefficient of 0 . 085BTU/sq . f t . ^F. The atrium is situated 
in the house so as to be adjoining as many rooms as possible and to 
permit free air movement:. The atrium is constructed of material 
which will give it significant thermal mass. The floor is five inches 
of concrete with a dark ceramic tile finish. The chimney is made 
from concrete block and finished with dark ceramic tile. There is 
also an ornamental pool of 300 gallons, which also provides thermal 
mass. All this thermal mass will be exposed to the sun's energy. 

This type of system will allow the benefits of a Trombe Wall, without 
its visual barrier. There is a rock storage unit under the atrium 
to store any excess collected energy. The rock storage is 200 cubic 
feet with a thermal capacity of 4,000 BTU/^F. The back up heat is 
to be supplied by an air to air heat pump. The atrium can also be 
used as a living space. 5N 

The sun's energy will enter the atrium, charge the thermal mass 
and heat up the air. The hot air will raise to the top of the atrium 
where it will be circulated to heat the rest of the house. WTien the 
houses temperature reaches 68^F the air will be drawn to the rock 
storage. At night the shutters will close, and the thermal energy 
in the atrium or the rock storage, will be used for heating. When 
necessary the heat pump can be used. These different modes of 
operation provide a great deal of flexibility in maintaining comfortable 
interior temperatures. 5N 

Since the house has not been built, as of yet, no actual performance 
data is available. Design calculations show that there should be a yearly 
savings of $183.00. (Fy 1977 dollars.) The investment to be amortized 
is $2,453.00. Calculations over the life of the system gives a pay- 
back period of about eleven years. This is considered the maximum 
amount of time to recoup the initial investment. As time progresses 
this payback period should get even shorter. The solar atrium gives 
cost effectiveness in heating and cooling that are equal to or superior 
to existing conventional systems. It has been designed to be convenient 
to operate by building occupants. 5N 
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LINDEBERG RESIDENCE 



The Lindberg House will have 2,100 square feet, with a heating 
load of 4.15 BTU/DD/sq. ft., in Minnesota, which has 8,400 degree days 
of heating. A hybrid passive solar system is situated so that it will 
be protected from winter winds and have large summer shade trees to 
help summer cooling. The house will have three levels, the garage and 
basement, the main living area and the bedrooms. The second and 
third levels are designed to give the maximum natural circulation of 
warm air through the house and the thermal storage area. The windows 
on the north, east and west sides have triple glazing, with double 
glazing and movable insulation on the south windows. There is also 
a skylight in the roof to give added collector area, which has a 
movable insulation system. 5F Exhibit 24 is a schematic of the houses 
heating systems. 

The hybrid passive system uses the direct gain concept to collect 
the solar energy^ and a rock storage system to store excess energy. 

The south face has 275 square feet of double glazing, with movable 
insulation and 75 square feet of south facing skylight , with movable 
insulation. The floor has 63,600 pounds of precast concrete for 
absorbing solar radiation. There is a metal absorber plate located 
above the rock storage bin and positioned to accept the solar radiation 
from the skylight. The rock storage bin will have 40,000 pounds of 
rock, with the air circulated through it by the furnace fan. 5F 

One of the problems with passive systems is the limited range 
it can operate over and stay within interior comfort conditions. 

The bedrooms will be put on the third level so that during the day 
•this level will be allowed to go in excess of 80^ F, thus allowing 
higher storage temperatures. Temperatures in excess of 80^''F on the 
third level should not put the second level beyond the comfort range. 
The metal plate collector will be used to further increase the storage 
temperature by the warmed air from living spaces being, forced by fan 
between the collector panel and the wall and into storage. The furnace 
fan will draw the air through the rock storage bin and either return 
it to the space through the concrete air floor^or recirculated via 
the collector loop. The precast concrete floor is an air floor, which 
allov;s flow both parallel to and perpendicular to the direction of 
the plank. Turbulence is created by supplying air to the Center 
location of an area and forcing the air to flow perpendicular to the 
length of the ‘plank, to floor grills located at the edges. The heating 
system operates by the following modes of operation: 5 f 

MODE //I Circulate air. This mode operates when an uneven 

temperature distribution exists to circulate air from 
the high interior spaces to the lower spaces and north 
exposed rooms. 

MODE //2 Solar Storage. I'Then comfortable interior temperatures 

have been achieved and solar energy is continuing to warm 
the spaces beyond their needs the excess heat will be 
drawn through the interior solar collector, and into 
the rock storage, returning via the concrete air floor. 
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MODE //3 Heat Charging only. This mode operates only when the 
interior spaces are at maximum comfort levels, \>Jhen 
this occurs the air will circulate through the solar 
collector rock storage loop. 

MODE //4 Heat from Solar Storage. When there is no solar radiation 
and the radiation from storage is insufficient to meet 
heating needs, the furnace fan will withdraw the heat from 
storage by the same circulation as Mode //2 . 

MODE //5 Auxiliary Heating. ^-Then the hybrid passive solar 

components can not meet the heating needs the auxiliary 
heating system can be used. A gas furnace is to be used 
as the back up system. 

The gas furnace uses its o\^m supply ducts. This is a partial 
duplication of the distribution^ but it was needed to prevent charging 
the solar storage with the furnace heat, making the solar storage 
ineffective to the lower temperature solar heat. A bin method analysis, 
estimated that the passive system could supply about 65% of the annual 
heating needs, which is pretty good with the available collector 
area. This design shows that a hybrid passive solar system can provide 
the major portion of the heating needs of a large well insulated 
house in a cold (8,200 DD) climate, and give relatively conventional 
comfort standards. There are plans to monitor temperatures, strati- 
fication of temperatures and heat flows for this house when complete. 5F 
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TECHNICAL ADVANCES TO AID PASSIVE SOLAR HEATING 



TRANSPARENT INSULATION 

The passive solar concepts collect energy through glazing 
to be absorbed in a thermal storage mass at low temperatures- With 
conventional building materials these concepts have draw backs, the 
glass will lose heat by radiation, the thermal storage will have to 
be massive because of the low temperature^ and the system will have 
temperature swings through out the days operations. Some technical 
advances that have been developed to help eliminate these disadvantages 
are the transparent insulation, the optical shutters, the movable 
insulation and the phase change material. The heat mirror will 
allow short wave solar radiation through, but be an insulation to 
thermal radiation. It is a thin transparent optical film which can 
be applied to glass or plastic glazing material. Depending on the 
application^ the thermal losses can be reduced by 25 to 75%. Normally 
multiple glazings are used to reduce thermal losses, which will 
also result in a reduction of solar radiation transmitted through 
windows. As glazings are added the transmission will continue to 
decrease. The ideal heat mirror for use in passive solar applications, 
should be as shown in exhibit 25. 

HEAT MIRROR 

The heat mirror coating will transmit solar radiation, .3 microns 
to 2.5 microns, and be highly reflective to long wave radiation. 

Some examples of characteristics that can be achieved by the use of 
a heat mirror with windows are given in exhibit below. 






BTU/sq. f t . /hr/ F Trans 

E=.2 E=.05 



Single Glazed 

Inner Surface .72 .63 .78 



Plastic Film + 

Air Gap Retrofit .31 .26 .75 



Double Glazed 

Air Gap .32 .27 .77 






cs 



Double Glazed 
Mid-Gap Plastic 

Film .21 .17 .66 
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In comparing the performance of windows with heat mirrors to multiple 
glazed windows, it was found that the heat mirror was about equivalent 
thermally to adding a glazing plus air space. Condensation and frost 
on a heat mirror surface has the effect, thermally, of eliminating 
the heat mirrrr. In cold climates heat mirror on single glazing 
would not be adviseable over double glazing due to the frequent 
condensation and frosting that would occur. 5S 

The cost of the heat mirror to the consumer is very important 
to how much gets put into use. The demand will affect how much is 
made and at what cost per square foot. Some estimates of cost to 
retrofit existing windows to consumers, assuming high production rates, 
are $.75 to $1.50 per square foot, if homeo\Nniers installed , and $1.50 
to $2.50 per square foot if professionally applied. Heat mirror on 
new windows should add from $1.50 to $3.00 per square foot to the 
retail cost. With the cost of adding additional glazing in the range 
of $2.00 to $4.00, it is apparent that heat mirror coatings can be 
potential competitors. The primary factor in decreasing cost is 
whether there will be sufficient demand to cut production costs. 5S 

SOLAR MEMBRANE 

The solar membrane is another type of transparent insulation. 

The solar membrane is made from a plastic film that has a unique 
combination of properties possessed by neither glass nor other 
plastics. The solar membrane has a transmission of solar radiation 
of 95%, which means less reflection losses than glass or other plastics. 
The transmission also falls off less rapidly with increasing angle 
of incidence, than it does in other materials. The solar membrane 
has a transmission of thermal radiation of only 13%, which is lower 
than all other plastics. Solar membrane characteristics allox^ many 
layers to be used and still be effective for use in solar energy. 

If four layers separated by a 3/4" air gap are used, the heat loss 
coefficient will be 0.21 BTU/sq. ft. /hr. F for upward heat flow, the 
worse case, and still transmit 82% of solar radiation. This com- 
bination of solar membrane can give a flat plate collector saturation 
temperature of 400 F without concentrators, vacuums or selective 
blacks. The solar membrane has characteristics that can make it very 
useful in the field of solar energy. The solar membrane is available 
at 35c per square foot in larger than 1,000 square feet quantities. 

This is the same as single strength glass, but because of ease of 
installation would be cheaper in overall cost. If the solar membrane 
is purchased in 200 square feet or larger quantities the price per 
square foot is $3.00. 5U 

OPTICAL SHUTTER 
CLOUD GEL 

The optical shutter is a material that will vary its transmission 
of solar radiation as a function of temperature. A plactic has been 
developed which will function in this manner. Cloud Gel. Cloud Gel 
will allow solar radiation to pass as long as the material stays below 
the set point. At the set point the plastic v^ill turn opaque white, 
thus reflecting all the solar radiation. I'Hien the material cools down 
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it will turn transparent again. The purpose of this type of material 
is to help pervent overheating that can occur in passive systems. The 
switch from transparent to white and back again will occur instan- 
taneously over a three degree celisus temperature change. The set 
point can be tuned to any value between 0 degrees and 100^ C, by 
adjusting the proportions of its constituents. Cloud Gel is inexpensive 
to manufacture, making its use in passive solar energy desirable. 5D 

PHASE CHANGE HEAT STORAGE 
THERMOCRETE 

Another material that can be effectively used in passive solar 
concepts is the use of phase change materials for thermal storage. 

The unique characteristics of phase change materials is that a large 
amount of energy can be stored at one temperature. The temperature 
will be that at which the material is going through a phase change. 

A problem with most material is that it will remain a liquid when 
cooled a bit below its freezing temperature, because of the lack of 
a seed for the crystal to grow from. Amother problem with most materials, 
in the liquid state, will be a container to keep things from getting 
messy. The problems can be overcome by incorporating the heat storage 
material into the very pores of a foamed concrete. The concrete 
will function as the seed to initiate crystal growth, and also serve 
as the container for the heat storage material. A concrete has been 
developed with phase change properties, Thermocrete, Calcium chloride 
is the phase change material for melting o: room temperature. Other 
salts can be used to produce melting temperatures between O^C and 
100 ^C, Thermocrete can store twenty times the energy as ordinary 
concrete. The compressive strenght of Thermocrete is one thousand 
pounds per square inch, approaching ordinary concrete. 50 

PARAFFIN BASE 

Another type of phase change material is the paraffin base 
material. There has been tests conducted to determine if a paraffin 
base material could be used effectively in a passive system. Results 
of test to date have indicated that the paraffin based material could 
be used to collect and store solar energy. The materials developed 
can store up to 250 BTU/sq. ft. Another benefit of using a paraffin 
based material is that when it is liquid it is transparent, giving 
it the ability to function as a window. There is on going work in this 
area, to try and produce an even better material. The paraffin 
based material will be inexpensive. 5X 

POLYHEDRAL WALL 

The Trombe Wall, with phase change material used for thermal 
storage, can use a polyhedral (honeycombs) plastic cover instead of 
glass. The polyhedral plastic will transmit more than 60% of the 
solar radiation and act the same as an inch of thick foamed plastic 
(R-4 or higher) . The polyhedral plastic will have to be. made from 
a plastic that will not be attacked by solar radiation. The polyhedral 
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plastic will be inexpensive, and when combined v;ith a phase change 
material gives an effective passive solar system. There is further 
experimental work going on in this area. 5 

HEAT MIRROR, CLOUD GEL, AND THERMOCRETE USED TOGETHER 

The heat mirror, Cloud Gel and phase change materials can be stacked, 
as in exhibit 26 , to produce an effective passive solar collector. 



On a sunny day the solar radiation will pass through the heat mirror 
and Cloud Gel and be absorbed in the Thermocrete. The Thermocrete 
will start to melt, when all the material has melted the temperature 
will raise causing the Cloud Gel to turn white, if set temperatures are 
correct. At this time the Thermocrete is completely charged, and 
prevented from overheating. The Cloud Gel will not be affected by 
outside temperatures because it is shblded by the heat mirror. The 
Thermocrete will thus slowly freeze, but stay at one temperature. If 
the Thermocrete melts at 70^ F and it averages 30^F outside, and there 
are two layers of heat mirror, in a 24 hour period only one third of 
an inch will freeze. Calculations show that on an average January 
day in Boston, there will be sufficient solar radiation to melt tv;o 
thirds of an inch of Thermocrete. Thus with this type of system 
it can be seen that even in the terrible New England winters it can 
maintain a comfortable 70"F.50 



Movable insulation is often used in passive solar systems to 
prevent the loss of heat when the sun is not shining. There are 
many kinds of movable insulation, ranging from very simple to complex. 
Four different types will be described to give an indication of what 
has been develpoed in this area. The night wall panel is simply a 
lightweight rigid insulator chat is attached directly to the window, 
when there is no solar radiation. The panel can be held in place 
with magnets or some other convient method. The night wall panel 
allows single glazing to be used, in a passive system, which will 
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allow greater solar energy into the house for use in heating. The night 
wall panel can also be used to keep solar radiation from overheating 
the house. Tests conducted with night wall panels, including edge effects, 
showed that an effective thermal resistance of 3.12 can be obtained 
from a one inch headboard, including the value for the glass. 5cc 

BEADWALL 

The beadwall uses tiny styrene beads that are moved from storage 
to fill a gap between window glazings when needed. A vacuum cleaner 
motor is used to transport the beads. The motor is automatically 
controlled so that when the sun is shinning the beads will be in 
storage. Exhibit 27 shows a beadwall system. There was no thermal 
resistance value given for the beadwall. 5CC 

EXTERIOR SHUTTERS 

Exterior shutters can be used many ways, and is an effective 
and often required insulation for passive solar systems to operate 
properly. The thermal resistance obtained will depend of the material 
used and how all the edges are sealed. An added benefit of increasing 
the solar radiation collected, can be obtained if the shutters are 
designed and used correctly. 

SELF-INFLATING CURTAIN 

The self-inflating curtain is an automated device capable of 
providing an efficient movable insulation for glazing areas, especially 
for Trombe Walls. The curtain is automatically put in place when 
there is no solar radiation. The curtain will then inflate from the 
hot air raising. The hot air will be intercepted by vents in the 
curtain, creating an increase in pressure and inflating it. How 
thick the curtain will inflate will depend on the mechanical stiffness 
of the material and the gravitional forces acting on it. The curtain 
will be stored in a roller above the wall. As the curtain is rolled 
up the air is evacuated through the side channels and the slots at 
the bottom of the curtain. The curtain consist of a number of layers 
of thin flexible material of high reflectivity and low emissivity. 

When the system is inflated it will reduce the heat loss through the 
windows. An important consideration is the area of glass to be 
covered. If a large area of glass can be covered by a single curtain 
it will give the minimum infiltration rate, by keeping the linear 
feet of crackage to a minimum. A small curtain tested had a thermal 
resistance of 9. This was considered low because it was a small 
curtain, which will give too high of a infiltration rate. It was 
conservatively assumed that a curtain that would be used on an actual 
system would have a thermal resistance of 12. The cost of a curtain 
20 feet by 17 feet is $4.28 per square foot installed. To determine 
the payback period a curtain was used on a Trombe Wall with double 
glazing in a 9,000^F day climate. The curtain was assumed to have a 
thermal resistance of 9, which is low. The payback period for Aspen, 
Colorado, with an electricity rate of three and a half cents per 
kilowatt hour and an increase in cost of power of 15% per year for the 
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next five years, is 4.2 years. Considering the improved performance 
and the cost of the curtain, this is an effective movable insulation. 

HEAT PIPE 

A technique that is being tested is the use of heat pipes in 
combination with a water wall, concrete wall and paraffin wax wall. 

The characteristic of heat pipes that is useful in passive solar system 
is its thermal diode effect, allowing the flow of thermal energy in 
one direction only. The heat pipe will take heat from one end and 
give off heat on the other end only, with the reverse flow prohibited. 
The heat pipe used in a passive system is shown in exhibit 28. 5 

The absorber plate, with heat pipes attached, will collect solar 
energy. The heat pipes pointing up through the insulation into 
the thermal mass will transfer the collected solar energy to the 
thermal mas^ when the absorber plate is warmer than the storage. 

At night or when cloudy the absorber plate will cool doTO, but 
because of the thermal diode effect will not drain the storage. This 
system will not need movable night insulation, because there is a 
permanent insulation material between glazing and storage. Another 
advantage of this system is that whenever exterior conditions allow 
the absorber plate to cool there will be insulation in place, unlike 
other systems. In test conducted the water wall gave the best per- 
formance, but the paraffin wax wall and concrete wall performance 
were close. \^en a water wall with a heat pipe was compared to a 
water wall without a heat pipe, it was found that up to more than 
twice the useful heating output c‘ould be obtained. The single glazed 
water wall with heat pipes had an output that was nearly twice the 
output of a simple, no night insulation, double glazed water wall. 

The use of a heat pipe in passive solar systems will need further 
work to determine the optimum configuration and number to use. The 
heat pipe appears to be a very reasonable devise to use in combination 
with thermal storage walls. 5 
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WINDOW 



WINDOW DESIGN 

Windows in passive solar design is very important. To get an idea 
of how important windows are in the useage of energy, there is about 
one third of the nations energy consumed in buildings (residential 
and commercial) . Of this energy consumed two thirds can be affected 
by windows. Therefore, the proper performance of window is very 
important in the field of energy used in buildings. The six energy 
control functions of windows are passive solar heating, daylighting, 
shading, insulation, air tightness and natural ventilation. Recognizing 
these six functions, one can begin evaluating how a window’s per- 
formance can be improved. The following strategies, in windov/ design, 
will be discussed concerning improvement of window performance; 
site design, exterior appendages, window frames, glazing, interior 
accessories and building interior. Site design is important because 
by altering local solar, wind or air temperature the thermal perfor- 
mance of a house can be improved. An example of a site design that 
will improve a v/indows performance is the windbreak upwind of a house. 
This stategy functions because wind direction in an area Is usually 
the same during seasons. With the winter-* winds coming from essentially 
one direction, a windbreak can reduce the infiltration rate and con- 
vective heat losses. The infiltration rate is reduced because the 
windbreak will decrease the wind pressure on the joints and crakes 
of the windows. The convective heat loss will be less because with 
less air movement across the window the boundary layer will be more 
effective in its thermal resistance. If a planted windbreak is going 
to be used, it should be located downwind a distance no further than 
one to one and a half times the building’s height. A fence windbreak 
should be located closer, and allow part of the wind through. By 
allowing some wind through the fence, the turbulence of a solid fence 
will be evened out. If the wind’s direction changes for summer, the 
cooling affect will be maintained. 5BB 

The use of exterior appendages is very valuable in window design. 
The appendages can give daylighting, natural ventilation, increased 
insulation, and shading, depending on how they are used and what 
effect is needed. During winter when the sun is shining the appendage 
can be opened allowing light and solar energy into the house. At night, 
or times of no sunshine, the appendage is closed creating an air space 
and along with appendage itself will be about equivalent thermally 
to adding another glazing. During summer the appendage can be used 
to shade the sun’s energy and allow natural ventilation. Reflected 
ground light will be admitted to give some daylighting. The exact 
type of exterior appendage will determine how effective each of the 
factors are in improving the windows performance. 5BB 

The frame design is important to window performance for several 
reasons. How much lineal feet of framing there is will affect how 
much infiltration is allowed into the house. The larger windows 
will let in more solar energy in comparision to the infiltration 
rate than a smaller window, but the amount of conductive losses will 
be greater in the larger window, unless night insulation is used. 

It i^ for these reasons that the sizing of the window area, along 
with window perimeter, should be carefully designed to give optium 
energy gain. The direction a window opens will be important in 
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allowing air to either enter or leave the house. Also of importance 
to window performance is the tilt of the frame. By tilting the 
window frame towards the ground the summer load can be reduced and 
the winters solar gain largely unaffected. The summer load is reduced 
for two reasons, the amount of direct solar energy allowed to strike 
the window, and the amount of solar energy allowed through the window. 
With glass the solar transmissivity is greatly reduced for angles of 
incidences greater than 57 degrees. In summer, with the sun higher 
in the sky, a greater than 57 degree angle of incidence can be obtained 
with the correct tilt. Also to be considered in the tilt is to insure 
that less than a 57 degree angle of incidence is maintained during 
the winter. If the above tilt can be obtained the amount of solar 
energy allowed through the glazing will be essentially unaffected 
during winter, because the sun is lower in the sky, and greatly reduced 
during the summer. 5BB 

The glazing is the last chance to stop adverse climatic forces. 
Single glazing has been greatly used throughout the United States 
and is the least effective in controlling climatic forces. \{hat 
can be done with multiple glazing will be discussed later. An often 
overlooked glazing material is the glass block. The glass block 
has some unusual properties that can be effectively used in energy 
conservation. The U value of glass block can be as low as 0.4^ if 
a double air cavity is used. The larger the face area is the better 
the performance will be, because of infiltration. The glass block 
also has substantial mass, which can be used to store some heat, 
which then can be radiated into the house at a later time. How much 
light is transmitted can be controlled, along with the direction the 
light comes out of the block. Therefore, the solar energy could be 
allowed to penetrate further into the room and onto more mass. I^Jliat 
kind of glazing is used will effect the amount of solar energy admitted, 
the amount of daylighting obtained, and the amount of insulation 
provided for heat loss. 5BB 

Interior accessories are numerous, and varies greatly in how 
they effect the performance of windows. One advantage to interior 
accessories is their accessibility. An interesting accessory that 
can be used on any south facing window is the Ark-tic-seal. The 
system consists of a guide frame with three film shades; one heat 
absorbing, one reflective, and one clear. During winter days the 
heat absorbing and clear shade are lowered. The solar energy will heat 
up the heat absorbing film. I'Jhen the air is warmer between the shade 
than the room temperature a vent opens, allowing the air into the 
living space. I'lhen the air cools down lower than room temperature 
the vent will close. At -night all three shades are lowered providing 
an effective insulator. The window can also be used during the summer 
for cooling. This accessory, by collecting and distributing solar 
heat at the window vice the living space, alleviates the problem 
of fabric fading and glare from direct sunlight. At night, this device, 
by providing a good insulator, greatly decreases the windows losses. 5BB 

The building interior will greatly determine of what use the 
admitted solar energy is to the house’s heating load. How the window 
is used in the basic concepts have already been covered. An idea 
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that can be used in just about any house, is to place some thermal 
mass by the window to store what solar energy is admitted into the 
house. Water can be used for this very effectively. The above 
strategies give some idea of what should go into the design of a.: 
window for a house, for them to be effective in energy control. 5BB 

THREE GENERATIONS OF WINDOW DEVELOPMENT 



The glass used in windows have gone through essentially three 
generations to date. The ’’first generation” is the flat clear glass 
used in most houses. If this glass is used as a single glazing it 
is not very effective for use with solar energy. It will allow 
the most energy through, but also has the lowest U value. Multiple 
glazing can be used to increase the U value, but the amount of solar 
radiation allowed through will be decreased. The ’’second generation” 
is' the use of heat absorbing and glare reducing tinted glass. Their 
function is to reduce excessive solar brightness , like putting 
sunglasses on a building. They can also have low U values. There 
is also some special low absorbance glass to give high solar energy 
transmission. This kind of glass is usually used in solar collectors. 
The ’’third generation” glassware those with reflective characteristics. 

A thin reflective film is applied to the glass. This type of glass 
can have a range of reflectance, absorbance and transmission properties. 
The successive generations of windows show the increased consideration 
of the solar energy in glass advancement, which should result in better 
performing windows. 5T 

ENERGY TRANSPORT CONTROL IN WINDOWS 
A new concept in window systems designs, uses a larger than 
normal interpane separation, with Venetian blind-like convection 
inhibiting and radiation-controlling arrays f Exhibit 29 is a 
schmatic of a between-the-glass convection-radiation control system. 
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More than one array can be used in a window . \Ihen more than one 
array is used, a higher R value is obtained and greater operational 
flexibility is achieved. The arrays can be arranged so that the 
maximum amount of solar energy is admitted into the house, and so 
that the solar energy is directed to the thermal storage within the 
house. In tests conducted with this type of window system, U values 
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of 0.1 BTU/hr/sq. in the nighttime mode, and 0.2 BTU/m/sq. ft./"^? 

in the solar acceptance mode have been obtained. Even smaller U values 
appear to be achievable. This type of window system can make the 
window, usually an overall deficit to the house *s heating load, an 
asset. This window system is a good example of what can be done to 
improve the performance of windows, a very important part of passive 
solar heating systems. 5R 

At present multiple glazing is the usual method used to improve 
the U value of windows. The primary disadvantage to using multiple 
glazing is that as the U value decreases, with added glazings, the 
precent solar radiation transmitted decreases. Exhibit 30 
shows how the precent solar radiation transmitted and U value varies 
for different number of panes. Exhibit 30‘ indicates that multiple 
glazings can be effective in controlling heat loss on the windows 
not used for collecting solar energy. As previously covered, there is 
work going on to try and improve the characteristics of windows for 
use in a passive solar system, I'Jhen the disadvantages of glass are 
overcome, with the advantages maintained, the performance of passive 
solar systems will be improved, 5 

WINDOW QUILT 

Window coverings can greatly improve a windows thermal performance. 
As mentioned earlier, there are many different kinds of window covers. 
Most coverings are designed to work with windows associated with a 
passive solar system. A concept that can be used on any window is the 
window quilt insulating shade. The window quilt has five layers, 
consisting of two layers of polyester batting separated by a vapor 
barrier of metalized polyester film with a decorative covering. The 
mounting hardware is designed so the the shade can be moved up and 
down and have the edges sealed. VJhen the shade is all the way down 
the perimiter of the sahed is sealed. The shade performance is as 
follows ; 

1. Window Quilt applied to double hung single glazed window 

a. U value 3.15 BTU/sq. ft,/hr/^F without window quilt 

b. U value 2.30 BTU/sq. f t , /hr/'^'F with window quilt 

2. Window Quilt applied to air tight fixed single glazed window 

a. U value 1.03 BTU/sq. ft./ hr/^F without .window quilt 

b. U value 0.23 BTU/sq. ft./hr/'’F with window quilt 

The payback period will depend on how consientious the user is, 
but should be from 2 to 8 years depending on weather conditions. A 
nice feature with this insulating shade is that conventional shades 
can also be used with it to make an attractive window. 5Q 
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THER>IAL DATA REQUIREMENT AND PERFORMANCE EVALUATION 



Perhaps' the weakest area in passive solar heating technology is 
in predicting a systems performance. There are simplified techniques 
and computer programs which are designed for specific passive solar 
concepts, but when combinations or hybrid systems are used the 
techniques have to be modified, which usually requires a working 
knowledge of both computer programming and heat transfer. To further 
complicate the problem, when passive solar concepts are combined their 
performance is not always the sum of the individual concepts. This 
lack of workable methods to simulate the performance of passive solar 
systems has resulted in passive solar heating getting a bad name, 
because of inadequately designed building being uncomfortable. The 
siuation is improving becasue of experience gain from houses with passive 
solar systems. Therefore, a standardized method of evaluating passive 
solar systems should be developed, so that systems designed \j±ll 
give at least adequate performance and preferably the maximum per- 
formance. The procedure may require a computer, but it should be 
a canned program that can evaluate any of the passive solar concepts, 
combination of, and or hybrids by simply changing parameters. Another 
problem v/ith simulating passive solar systems is the need for hourly 
weather data, which is very time consuming. There has been a procedure 
proposed by D. W. Low to greatly reduce the required weather data. 

The procedure reduces the years hourly weather data to three 

"weeks” of hourly data. This procedure has been tested and can 

produce results that are 99.6% accurate at about 6% of the cost of 

100% accuracy. A description of the procedure is presented in reference 1. 

Conceivable if a canned computer program as described above, could 

be developed along with reduced weather data, the optimum passive 

solar System could be determined for a particular area. 

STANDARIZED EVALUATION PROCEDURES 

The National Bureau of Standards, at the request of the Department 
of Energy, is putting together a document that will standardize the 
evaluation of passive solar systems. The scope of the documents 
is as follows; 

1. Establish pertinent performance factors. 

2. Determine the type, accuracy, range and frequency of required 
measurement. 

3. Identify appropriate instrumentation. 

4. Recommend data analysis and reporting techniques 

From the results of the project the evaluation of the passive solar 
system should be able to; 

1. Determine energy savings of fossil fuel and electrical power. 

2. Determine the fraction of the building’s hot water, heating 
and/or cooling load contributed by a passive system. 

3. Evaluate passive components thermal characteristics 

A. Determine the building’s comfort level 

5. Determine the extent of occupan interaction required for 
functional performance or other operational requirements. 
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6. Determine the system degradation for the measurement period, 

7. Obtain insolation and other ambient environmental data for 

designing and evaluation purposes. 

For further details on the project refer to reference 5V. This kind 
of work is the kind needed to get passive solar systems more fully 
accepted, and should result in easier designed systems that provide 
acceptable comfort conditions to homeowners. 5V 

SIMULATIONS 

SIMPLIFIED TECHNIQUE /CROSLEY HOME 

A non computer method of predicting the annual performance of a 
passive solar system with minimal heat storage has been performed on 
the Crosley Hone, by Andrew M. Shapiro. The method uses average 
monthly sun and temperature conditions from available weather data. 

The month is then broken down to full sun days and non full sun days. 

It is considered that on full sun days there will be an excess of solar 
energy available. The average monthly and daily heat loss is calculated 
using standard procedures. For a full sun day one starts at a given 
hour, and calculates the heat lossed and gained during that hour. If 
there is any extra energy available the surplus is carried to the next 
hour as heat stored. All heat storage in the house is treated as 
an isothermal mass at the temperature of the interior of the house. 

When there begins to be a surplus of energy the house will begin to 
raise in temperature. The greater temperature will cause a greater 
heat loss, because of the greater difference in temperature. The 
increase in heat loss is tabulated separately from the heat loss if 
the house stayed at 20^C. For the next hour the heat loss due to 
increasing the interior temperature is subtraced from any surplus 
energy available, whatever energy remains is added to the next hour. 

This process continues until the heat storage reaches full capacity. 

At full capacity any additional solar energy is not allowed into 
storage (dumping) . Losses due to any increased interior temperatures 
(overheating) or dumping are totalled separately. At this point the 
only sun heat being used is that which is keeping the house at 20 C. 

If a person did not want to dump as much solar energy, additional 
thermal mass could be added to the house. After the sun goes do\<ra the 
heat stored in storage will be released until room temperature is 
reached. The heat from storage will be equal to the heat loss for that 
hour, which is the heat loss at 20 and any heat loss due to the 
interior being at a higher temperature. One full sun day is done 
for each month using average mac-zimun, minimum and mean temperature, 
and average dear-day radiation. The average monthly gains and losses 
over the heating season are tabulated as follows: 

1. Average heat loss in absence of sun and average gain from 
sun are figured from the average daily figures, for loss and gain, 
multiplied by the number of days per month. 

2. Full sun day totals for overheating and dumping are multiplied 
by the numbers of clear days in the month. 

3. The results of 2 subtracted from the monthly average gain from 
the sun will give the useable sun heat. 
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4. The useable sun heat subtracted from the heat loss in absence 
of the sun will give the needed supplementary heat needed. 



Exliibit 31 is a graph of the above procedure used on the Crosley 
home for a full sun day in December. Although this is a simplified 
method to simulate a passive system, it is obvious that it will be 
time consuming to perform. Considering the cost for the time to work 
this method makes the computer simulations more attractive. The problem 
with computer simulation is that they are as of yet not flexible 
enough to be used on any system but what they were designed for. The 
computer programs can handle some modifications but it will usually 
involve modifying the basic program. This kind of modification will 
require computer programming and heat transfer knowledge. A computer 
program that can give results for many different passive solar systems 
by modifying input parameters is very much needed. 

COMPUTER SIMULATION /PASOLE 

PASOLE is a computer program to predict the performance of a 
Trombe Wall or a water wall. This program was select to work with 
to gain the experience of working with a canned computer program, 
and to see what kind of performance could be obtained from a Trombe 
Wall in Corvallis, Oregon. Weather data for the program was collected 
for the month of December only, for the year 1978. The winter of 
1978 in Corvallis was a severe winter. Therefore, the percent solar 
load from the computer run will be close to the minimum that could be 
expected. A copy of the weather data used is in appendix 1. Tom 
Wilson, from Oregon Department of Energy, was consulted to find out 
how to use PASOLE. Tom Wilson has used and modified PASOLE to its 
present state. 

A copy of PASOLE is in appendix 2. The program can be run from 
cards or a terminal . The terminal was selected for tv/o reasons; 
to learn how to use the terminal and because the terminal uses less time 
between runs. PASOLE has been constructed as the source program 
with PASOLGO as the object program. Before PASOLE can be used two 
files have to be created. A weather file and a command file. The 
weather file was created from cards in the following format; 



COLUMN 




DESCRIPTION 


UNITS 


8-9 




Year (19 ) 


— — 


11 - 12 




Months (1 - 12) 


- - 


14 - 15 




Day (1 - 31) 


- - 


17 - 18 




1 Hour (0 - 23) 


- - 


20 - 23 




Total solar radiation 
on a horizontal surface 


W/M 2 


28 - 29 




Wind Speed 


M/ sec 


31 - 35 




Dry buld temperature 


deg. C 


Hourly weather data 


was collected and put into a file 


named weather. 


The weather 


file has 


to have data for 24 hours a day. 


The first 


hour of the 


day has 


to be zero hour with hour twenty 


three as the la 
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hour or the program will not run properly. The command file can be 
easily created from the terminal, and has the following command inputs: 

CONTROLITIP (next line gives minimum and maximum control r.ode 

temperature for night and day) T-Min_nite, T-max-nite, 
T-min-day , T-max-day 

D LONGITUDE (next line gives difference in time longitulv and 
actual longitude) 

HOURLY 10 (next line specified starting and ending days for 
hourly output) MMDDY, MMDDYY 

LATITUDE (Next line gives latitude in degrees) 

MONTHS (Next line gives the number of months the simulation is 
to be run.) 

RESNIGHT (Next line specifies R-value of night covering insulation 
for wall. Units are 1/BTU/^F/sq. ft.) 

STARDATE (Next line specifies the starting date of simulation, 
MMDDYY format) 

U LOAD (Next line gives the load value of house normalized by 

the collector area, house load/collector area (BTU/hr/^F/ 
sq. ft. of wall).) This value is usually about one. 

VENTS (Next line gives vent status, 0 = no vents, 1 =- vents always 
open, 2 = no reverse flow, 4 = Thermostatic vent control. 
Default = 0 

AUX COOLING (specifies whether auxiliary cooling will be used 

On the terminal a number of command files can be created so that 
runs with different parameters can be easily run. An example of a 
command file created on the terminal is as follows; 

Note; Underlined material to be entered at terminal 

1. Log onto terminal 

Control A 



2. 


/EDIT, COM (Com is 


command file name) 


3. 


Begin text editing 




4. 


? A 




5. 


Enter text 




6. 


? / CONTROL TMP 




7. 


? 60. 75. 65. 


75. 


8. 


? D LONGITUDE 




9. 


? 3.4 




10. 


? HOURLY I 0 




11. 


? 120178, 123178 
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12. 


? 


LATITUDE 


13. 


7 


44.55 


lA. 


7 


MONTHS 


15. 


7 


1 


16. 


7 


RESNIGHT 


17. 


7 


3 


18. 


7 


STARTDATE 


19. 


7 


120178 


20. 


7 


ULOAD 


21. 


7 


j. 


22. 


7 


VENTS 


23. 


7 


2/ 


2A. 


Ready 


25. 


7 


END 


26. 


End text editing 


27. 


$ 


EDIT, COM. 


28. 


/ 


SAVE, COM 



Now the command file COM is created and can be used in PASOLE. The 
data for the commands can be loacted anywhere on the next line for 
that command. 

With the weather and comjnand files created PASOLE can now be run. 

The following steps are used to run PASOLE and have the results printed 
on the computer printer. 

Note: Underlined material to be entered at terminal. 

1. Lot onto terminal 

2. / SETTL, lOQ (sufficient for one months data) 

3. $ SETTL, 100 

4. / GET,\>7EATHER,C0M 

5. / GET, PASOLE /UN = (file // where PASOLE stored) 

6. / FTN, I = PASOLE, L = 0, B = PASOLGO , OPT =2 

7. / SAVE, PASOLGO 

8. / ATTACH, IMSL/UN = LIBIU\RY 

9. / X, LIBRARY = IMSL 

10. Library = IMSL 

11. / TITLE (PFIL) Put in title desired for output 

12. / PASOLGO, PFIL 

13. X CP seconds execution time 

14. / ROUTE, PFIL, DC = PR 

15. Route, Complete 

16. / BYE 

At this stage PASOLE will be run and the results printed on the printer, 
with the input of the weather and command file used to give the 
output. One month of data will take about 43 CP seconds. The ULOAD 
and RESNIGHT conimands where considered the most important as far as 
system preformance was considered. The following values for ULOAD 
and RESNIGHT were used to obtain the precent solar load from PASOLE. 
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ULO.YD 


RESNIGHT 


PERCENT SOLAR 


1 


1 


3.82 


1.5 


3 


7.89 


1 


3 


11.53 


.5 


3 


20.55 


1 


9 


18.61 



The actual printouts for these runs are in appendix 3. The above 
results show that the percent solar can be changed greatly by changing 
the values for ULOAD and RESNIGHT. The lower value of ULOAD in an 
indicator of a tighter house as far as infiltration is concerned. 

The higher value of RESNIGHT means there is less loss of heat at night. 

To maximize the performance of a Trombe wall, ULOAD should be a minimum 
and RESNIGHT a maximum. 

The PASOLE program can simulate a Trombe wall or a water wall, 
but to change from Trombe wall to water wall requires modifying the 
base program. There is an on going attempt to create, from PASOLE, a 
program similar to F-chart. This program will be of greater use for 
designing passive solar systems. A description of PASOLE is is Appendix 4, 
it is essentially a thermal network solution. 

EVALUATION OF PASSIVE SOLAR HEATING SYSTEM PERFORi'IANCE 

Once a passive solar system has been built, how does the user know 
how well the system is performing. Bristol L. Stickney has devised 
a simple method to evaluate a passive solar systemfe performance. The 
method requires only a pair of maximum minimum recording thermometers. 

At present the method is applicable to winter heating only. The following 
are definitions of symbols used in the evaluation process. 

ITMAX-Indoor air temperature maximum 

ITMIN - Indoor air temperature minimum 

OTMAX - Outdoor air temperature maximum 

OTMIN - Outdoor air temperature minimum 

ITR - Indoor temperature range 

ITMAX - ITMIN = ITR (equ. 1) 

OTR - Outdoor temperature range 

OTM\X ~ OTMIN = OTR (equ. 2) 

RR - RANGE RATIO 

ITR/OTR = RR (eq. 3) 

OTDY (Delta Tee Daytime) The indoor-outdoor temperature difference 
during day 

ITMAX - OmVX = OTDY (equ. 4) 
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OTNT (Delta Tee night) The indoor-outdoor temperature difference 
during the night 

ITMIN - OTMIN = OTNT (equ. 5) 

OTAV Delta Tee average 

(OTDY 4- OTNT) /2 = OTAV (equ. 6) 

ITAV - Indoor temperature average 

(ITllAX + ITMIN)/2 = ITAV (equ 7) 

OTAV - Outdoor temperature average 

(OTMAX + 0TMIN)/2 = OTAV (equ 8) 

The thermometers should be placed in a shaded area, inside and out- 
side, and reset so that no previous measurement will be recorded. 

They should be set up in the morning betv;een 7:00 A.M. and 9:00 A.M. 
and left for 24 hours. At that time ITMAX, ITMIN OTNAX and OTMIN 
should be recorded. The data can be taken for any day, but a cold 
sunny day will tend to give the most meaningful results. A cold 
sunny day will tend to drive up the indoor temperatures and tax the 
limits of the storage mass. 5W 

Once the four initial temperatures are obtained, the remaining 
eight parameters can be derived from quaions 1 through'' 8. ITR and 
OTR establish the limits of the indoor and outdoor temperature 
fluctuations. RR establishes a reference number that indicates the 
quality of the buildings response to these fluctuations. OTDY, OTNT 
and OTAV indicate the building response to raises, drops and fluctuations 
in outdoor temperature respectively. ITAV is an indicator of the 
average comfort level indoors and OTAV indicates the average severity 
of the outdoor. The absolute accuracy of these parameters are not 
necessary. Once the eight parameters are obtained interpreting them 
is difficult. Bristol L. Stickney has devised a chart to easily 
interpret the meaning of the parameters. Exhibit 32 is the chart 
devised by Stickney. 5W 

The RR and OTAV values are used to determine how the system is 
performing. The OTAV value is plotted on two parallel lines offset 
to give the limits of fluctuation on the OTAV scale. The OTR scale 
is used to set the upper limit of the RR values. The OTAV scales 
and OTR scale will line off a comfort zone, zone A, zone B, zone C, 
zone D, zone E and zone F. By plotting the point where the RR value 
and OTAV value cross, the performance of the system can be determined. 

The best performance would be the center of the comfort zone, but as 
long as the point is within the comfort zone the interior temperature 
stayed within comfort standards (60 to 75^) khat is indicated if 
the point is within one of the other zones is as follows: 

ZONE A - RR too high, OTAV is too low 
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The thermal contact between the storage mass and the collector 
is inadequate 

Add more storage mass and/or more sunlit storage surface area 
Add more insulation to improve the inadequate thermal 
envelope of the building, (i,e. night shutters, roof and wall 
insulation) 

ZONE B - RR is too high OTAV is O.K. 

More §torage mass and/or sunlit mass surface area is needed. 
Venting overheated air and/or shading the south of the building 
may be necessary. 

ZONE C - The RR is too high. OTAV is too high. Same as ZONE B 

ZONE D - The RR is O.K. OTAV is too low. 

The building needs more collector area. 

A better thermal envelope is needed. 

Add insulating shutters and/or more in the walls and roof. 

Remove any southern shading devices. 

ZONE E - The RR is low OTAV is O.K. 

There is no solar gain. Either the building is being heated 
by storage discharge, or a thermostatic space heater is in 
operation . 

ZONE F - The RR is O.K. OTAV is too high. Same as ZONE B 

/ 

This system can be used on any type of passive solar system in any 
location, and is very easy to use. This type of evaluative procedure 
is needed if passive solar systems are going to be generally accepted. 5W 
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The thermal contact between the storage mass and the collector 
is inadequate 

Add more storage mass and/or more sunlit storage surface area 
Add more insulation to improve the inadequate thermal 
envelope of the building. (i.e. night shutters, roof and wall 
insulation) 

ZONE B - RR is too high OTAV is O.K. 

More gtorage mass and/or sunlit mass surface area is needed. 
Venting overheated air and/or shading the south of the building 
may be necessary. 

ZONE C “ The RR is too high. OTAV is too high. Same as ZONE B 

ZONE D - The RR is O.K. OTAV is too low. 

The building needs more collector area. 

A better thermal envelope is needed. 

Add insulating shutters and/or more in the walls and roof. 

Remove any southern shading devices. 

ZONE E - The RR is low OTAV is O.K. 

There is no solar gain. Either the building is being heated 
by storage discharge, or a thermostatic space heater is in 
operation . 

ZONE F - The RR is O.K. OTAV is too high. Same as ZONE B 

/ 

This system can be used on any type of passive solar system in any 
location, and is very easy to use. This type of evaluative procedure 
is needed if passive solar systems are going to be generally accepted. 5W 
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CONCLUSION 



At this point it would normally be approriate to present 
why passive solar heating is the desirable way to go. To do this 
would mean trying to point out how passive systems are better than 
active ones. This type of thinking has created two sides, when there 
should only be one. Initially separating the solar concepts into 
passive and active was probably good, in that a concentrated effort 
could be put into each side. As stated by Richard S. Levine, it is 
time to start using an integrated approach to solar designing, one 
that would take the best from both passive and active systems to 
come up with a 100% solar heated house, except for pumping power and 
controls. Our present state of design is like the first powered 
vehicle, the horseless carriage. The horseless carriage had an 
internal conbustion engine on a vehicle that use to be pulled by a 
horse, which was an awkward arrangement. With time an integrated 
approach was taken and the Model T came out, which is what is going 
to have to happen to solar design. An example of how a non integrated, 
approach can lead to not the best design is the F chart program. 

The F chart program is very popular and accurate within the range it 
was designed for, but it is only designed for active systems and 
gives accurate results only for solar systems that will provide 
up to 70% of the heating load. These assuptions are obsolete, 
because it is possible to design a cost effective house that receives 
all it’s heating, cooling and hot water needs from solar energy. 

It is therefore time for solar design to get out of its infancy and 
start advancing towards an integrated approach. Once this design 
becomes accepted the two sides can work together, and start to quickly 
advance solar design. The Raven , Run Solar House has been designed 
from an integrated approach to supply 100% of the heating from solar. 
This house has active and passive collectors side by side on the south 
wall to optimize the solar collection. The integrative design concept 
goes beyond maximizing solar collection and storage to what the house 
is made of, how it is constructed and site conditions. At present 
the biggest problem with integrative design is the design itself. 

There will have to be simplified design procedures developed before 
solar energy systems become generally accepted. A possible hindrance 
to integrated design being accepted is that passive systems have been 
thought as for architects, and active solar systems for engineers. 

This kind of thinking will have to be rejected, with both engineers and 
architects working together. Together is the key word, everyone 
working towards a common goal. 5A 

There has been some advancement of materials that can improve 
the performance of solar systems. These materials can be very use- 
full, but unless solar design becomes generally acepted there will 
be no market for them. At this stage then, it x^ill be important 
to take the steps that will begin to make solar design accepted. 

The more accepted solar design becomes, the more incentive there 
will bo for improving materials for solar systems. A weak incentive, 
by the governments, is the tax incentives. The problem with tax 
incentives are they encourage systems that supply the minimum 
percentage of the heating load by solar energy. Wien just the 
opposite should be encouraged, the use of larger solar percentage units. 
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The one area where an integrative design would not be applicable 
is in retrofitting existing houses. The active systems are probably 
best suited for use on houses already constructed. Passive systems 
can be retrofitted, but their performance will not be equilivalent 
to new construction, as it can easily be v;ith active systems. 

Solar energy can supply a major portion of household heating needs. 
How much of the nations household heating needs will be supplied 
by solar energy will depend greatly on the direction taken by solar 
designers. Exhibit 33 shows how available energy is used in the 
United States. 2 The flow of available energy shox^s that the house- 
hold and commerical area is a very bad user of available energy. If 
solar energy systems could become generally accepted it could consider- 
ably improve our nations use of available energy. This should be 
a considerable incentive for both government and the private sector 
to push hard for solar energy use in households and commerical 
buildings. At present there is a good start, but solar design is now 
at an important crossroad, if the integrajtive, approach is not taken, 
solar energy systems will not be very successful!. 



/ 
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ALL VALUES ARE AVAILABLE ENERGY. VALUES REPORTED 

ARE lO^^BTU. EXHIBIT FROM PAPER" AVAILABLE ENERGY CONVERSION 

AND UTILIZATION IN THE UNITED STATES" G.M.REISTAD 
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PAS0Lf-:(JUrPUT,lN^UTtTAPcL,TAPE2) 

Ol.!?.?-} CHANGE GECLi.NAricU FORMULA TO 0. AND. i. INSTEAD 

OF L^iiL '"0 ";HJ_h 

01»(Jti»7D aJO jtVtRAL PAr%MNETER COH'lM‘jOS» ■\£Ai\ KAN GE 

A G S I G N I'l N r OF G r <. H u A 1 0 P K A N E F E K S 

-- 12.E9.73 NGOiFr FD ’03E—A 1 r; J g PEE D ' I N M/G' - ChANGES'IN 

FIAwD .node. A„D in CAlOUumTIONS (G.^- ..NSTEAD OF I 3 ) 

DON NON / JLADl/NVfNFtNT *NNA \ ,<Fl3Q) »iC»wV(pj) i'_r (siu) »i_T(GG) t 

1. (\lL»rl.NGiKi.'-*t xF NAAtiic..\'\fU'i.i o»\hUAuw »KCAwC tNOc^F t*sCo 10 
C JNnuN / 5lmJ 2 / M' iP.E,MFiNE,ir:NE,JATE,JMr,D«fl,jAf2,NJAr<i2), 

1 NO , NCI , Yr;, N(Nu, .nGAY 1, INDA I E 

OONNuN/ S LAUG/F(pG),TJ(5j) ,D?N(50),;> (G3),SP(5a),00^0(5G,5Q) , 

I GCuN ( 50 ) ,uFADO ( p 0) ,E.K<F ( EG) , FN ( :> Q , 5 0 , 2) , FD (5 J , EG , 2) 

— c 0 N N c N / j L A 0 1 / A ( G 0 , 3 0 ) , G < 5 0 ) , s 0 SHOT 1 5 c r, 3 u GK c 2 c E a ) ■ r SO CO r? rt 5 urr 

1 G0o0.N2(50),uVHi fJi/CL»DAHT»., AG;_»QMZfOlN3»OT‘\Af< 
G0.N^1C^^/3LMub/G0NKTt0srf^N^Q^^II^^J^^>lO|/'^,Aut 
1 PI»7xN::A»0£Gf\,AJ»T3A5Et0Dj r.'lAA»FNiFJ»T3A.K,FACI.KENAIN 
DONnGN/ 3l^0o/...GKD ( 3 G» , 0 CON ( 3 G ) , 40 1 2 ( 1 0 0 J , SuC 12 1 ( 1 0 0) , 

1 34C122 ( 100 ), T CMi,<, TCNAX, FS NI ti.N, TCMliiOt rO.lAXN, TCNA.^J tTCOOi. - - — 

0ONrlu.N/’3LM07/rPEKJAYfi3C0f<»r AN3*0CI NT 

CONNON/ 5^AO i/iOA Y , IOPEnYX, InDmY, < Jh T , OHCON V , hL AT , 0. CNG, 

1 HiNJEXf 7 CL r fuE3 t OUSDES «S1N0LC« GINgAT tOOSLAT tOSP 

2 |Nl^3^^1 tOIGUl ,KgH(JT, KDCOl., NOgHI , NO SH2 

Cu ,'lNUiN / 3 LmOG / T A J ( 5 0 ) » VELO ( 5 0 ) » 4.N 0 ( G ) f .RH Q ( 5 0 ) » N J3 » 13 0.. { 1 5 ) t 

1 FIw.r(15) ,WA2xN(l:>) ,AG.Z(l5) ,NG^Z(i5) ,.UA,IAH3(10), 

2 AuFh( 15) »Ui-lMiNG(l5l »0GxPR(13) »TG^Z(lp) »JA3G(1^) »IS0uX(24) 

C0MNuN/3LAQxG/4jNNAAfr^.DxFF(?)tEA»KHC»TZEKJT — - 

1 NCuNN, IIOO.N (130 ) , I23GN (130) »U(i0'0) ,AOON (100) ,UG( 100 ) ._UM (lOQ) 

CON nOiN/gLmOII / TILT G»AZlNO».KiGO»COGTO»GiNTC»FSGOtFSG3»Fs£S Nil 

DON NUN / T KGNSE/VolF < kCGN (50 ) *.K\/£NT « IAIR» JAI.R»KLHGTh»RHOSR* 
lASP.\MT»FCSrx»JwA»TKGGAT»JCN»KWALL 
GO HN ON / I OFu GS / i S .RCP .R> x 3 1 » 102 1 1 PR I NT H , IPK I NT 0 » iPRI NTH » xP.K S N » KHEOP.R 
CONMO.n/ gxNEO/W<i»REA (1G3 ) « IH<mKcA ( 100) 

G 0 N I'iO'N / F I N E o / i N X G ri T » F ^iO RN t i c. V GjN t 0 E G T 1 » Oc. L I t 7 X tN E ■ ■ ■ - - - 

_ 0 1 H £N S I D N L OA Y_( 2) M i (12» 2 0) ,NGY( 12) ,SU,NTOT (20) 

INTEGE.R day , C A Y1 , Q A Y2 , Y R, CA FE , UNIT S 



DA TA( N0AY( M ), H =l,12)/31,28,3i; 30, 31, 30, 31, 31, 30, 31, 30, 31/ 
DATm NM a a , P I » DEGR AO /5 0 , 3. Ih.1 3 ^2 OP-+ » . 0 174 332 93 / 

- PRINT real time 

4 A L w" D A r E 1 ( A RG t) •' T" CACr~Ct;'3 0 K rniiG^) 

PRIf)T 2999, ARG2, AKGl 



OiT A..N INPUT OAT*. 

CALi. INDmTA 

— —INirxMi. C A 1 .GU L A T I 0 N3 
G ALi. FRl nE 
Kc Rpi = C 



RE Hi NO THE WEATHER OATA FILc. 

i%.cWINu 1 . .. 

QVH i'=G\/ :l=oaht=4acl=dd=goh= 0 . 

QdSuM = QI SUN- 0 , 

00 220 I=1,NM«,< 

S-OS.-^Cl ( I ) =SQSr%02 ( * ) =3QCON 1( I ) = SOCuN 2(1) = 0 . 

220 CONTINUE - 

00 225 J=1,NCCNN 
225 30C121 ( J ) =S4C122 ( J ) =0 . 

-- fO,-< hCkYI gays O.NLY (nGAYI .£s5 THAN 32) 

IFLol=0 S IF (nOA Y 1.GE.32) GO TO 220 
IFwol=l 3 f,NU=l £ lGAY( 1) =INJmY«-NDAY 1-1 
IF ( uOxY ( 1) .l-E.NDA Y (HOD ) GO TO 223 

lOA Y (2) =LCm Y ( 1 ) -t.oA Y (NO 1) £ lU A f ( 1 ) =N0 A / ( MO I ) 5 NM0=2 

223 CONTINUE - '■ 

N0=n 01-1 £ H7iNE=0 £ NCALC=0 5 NSTEP=3 

START FI.1E STEP .OOP 

DO 1200 NON 7h= 1, (,MO 

HOrMO-H S IF (MO. £U. 13) MO = l 

rCCOOL = 0 £ I F ( ’\£i N X . L7 . 0 . 0 1 ) GO TO 23 7 

Kv^O*jw“l £ XF(rlO.GE.HOSHl.MND.*NO.w£.MOSH2) KCOOL“*l 

0AY1=1 £ IF ( month .£Q. 1) CAYl = xiNDAY 

0Ar2=N0A Y (NO) S IF (IFLSl. EU. 1) D A Y 2 =lD A Y ( NON TH ) 

DAY wCGP 

00 1100 DA Y=DA Y1 ,OA Y2 

FI .lt=-OEL71 5 NFiNEl = NriN£ 




lL)Air = iD; Y n $ IF (lOAY.Gr. I0AY=1 

CmLu 3AYLY 

u rwLi=b. 

IF ( jATE. £U. iOl ) iPRi-NTHr 1 5 I F ( J A 7 E . ECi . I o2 ) I PKi Th = 0 
ir ( .Ok, (IP^i.ua.EQ.lJ ) pKiri7 hQGO,CA7£ 

IFUHxIMTH.EQ.OJ GO TO Z\b 

J' pKlf<T hEAutKi> FOK HOU,x_Y pKlf<7Ii<L> 

IF(KH‘-L_«EU«0)PkiNTA0Q't 

iT (kkAi .EO.i) PkInT ‘♦004 

IF ( KkAi.^ .EQ .2) PRINT 4002 

245 JJNIINJE - - 

% 

/“““■■ T I -iE It^GHEMENT LOOP 

■— 00 IGOO ITHE = 1, NT! MEl 

u7Ifl£ = I 7 lME-1 S IF (i_71'lE. EQ. 0 ) LTiM£ = NTIME 
.17 IrlE = N7 IRE+ 1 

0£l1=G£^T1 S 7IN£ = T INE^OElT i T 1 R£a = 71 24 . / TPEROA Y 
<IOniNG = 0 

aHSuM=QHSUM40ELT»^HZ S OISJH=uI I uH^OElT^OINC — 

•KNJ-l S iF<7Iiit,\.G7.rR0RN.ASJ»7Ifl£A.L7.r£./EN) KNC = 2 
KSriUT = -l S XI- ( KND.cQ. 1. ANC. rEGNI . G7 .0 .0 1) <3r(UT=l 
; ■<3HU7 = :<SHur*’KOOOL • 

i SET MmA and min temperatures for THIS TIME PERIOD 

IF t ( f IME a .UT.THGrN ) .AMD . ( TiMEX ,lT .TN l^HT ) ) GOTO 255 
250 TSMiN: TSHIN n £ rCMAX=7GMAAN S GO TO 2o0 

255 7GM4.N = TSHIN0 S TG MAa = TGMAxO — 

26 0 GJN.'INJE 

rG3jL=TGHAA S TCMAXl=TGMAx 



CAuu SUNSRC 

IF (kAUXGL .EQ.O .ANG .TAM3 .GT. TCmAa) TGMAX = TAM3 

IF (isIG.EQ.-l) T(I^)=TGMIN 

XF(RIG.£Q.+1) 7(iC)=TCMAX •- - - 

1F(kI..£Q, 2) T(^C)=TGHIN 
I2d 0 = G 

'2o0 SuNTINUE 

I2dj=1230+1 ' 

XF (iZaO. GT.5) GO TO 1675 . . 

N S 7 c P = N i T E P <• 1 

; irERnTiUN LOOP FOk TEMPERATURE OEP£NC£NCE 

00 -jOG XTER=1,ITMAX -- - - — 

NCAlC = NoAlC Pi 

; GET uP'if iP,GONC,T (FXXED) , SCON - GOUlO 3E F(T) 

sal- prop 

iF ( mTIME.EQ. i) GO TO 800 

) GALGU-A7E COEFFICIENTS 3(0) FOR HOCE EUJATXONS 

00 tOG IV=1,NU 

1 = LU (lU) 5 SU.MKr = 0. — - 

UO SOC JF=i,NF 

J=lF (OF) 

100 3'0MK.r = i0MKT ♦CUNu ( X » 0) ‘*■7 70') 

3(IV)=C?M(I)*T0(I)/0EL7+FAC»(S(I)+SUM.<T) + (1 ,-FAC) *GPHt I) *QFAC0 (1 ) 

IF ( KCAlG . £U »2 ) GO TO 400 
JO 32u 0V=1,NV 

0 = L rf ( 0 V ) 

A ( I V , OV) =-FAC*CCN0( !♦ J) - ■ 

IF ( X .EG, 0) A ( W, ON) =GPM ( I ) / 0 Ei.T pF A C * ( 3 C ON (I)-SP (I ) ) 

320 Continue 

--40 0 CONTINUE ■ ^ 

XF (kCh-G ,EQ. 2) GO TO 460 

: SJLYE EOUmTIONSX 3JMC A( I , 0 ) *T (0) ) =3 (I) ; 1=1, N'T 

I FOR T(0) ” ANSWEkS in 3 ARRAY 

THE MAT\IX A iS MuOIFI^O af THE iMSi. SOiROOflNES. 

; XF TmE M^TklX A 13 NOT SYMMETklC, THE FO-l.uWING GOOE MUST 3£ USEI- 

I CA-L L£UTiF(.->,l,N0,NHAA,3,0,WKAREA,IER) 

I - ■G'jRRc.NTLY, symmetry IS o-VEn HnE.'i COnJ ( i , 0 ) - 1 S— E Y MM ET R IG A nO~ T HE 

I Fj-uUHINu 15 APPLiCAouE 

CmLu \/COTFS(A,N02NMAX,A) _ ^ 

CA — u — £0lE(A,N0,a,l,N.^IAA,u, I>*,KMREA,Wr\AxLM,iEr*;.) 

1 CHEoK Tj eEE if THE An3Hcp,S HmOE CONVERGED YET 

N£ Rk = 0 

JO ‘T'tO iv = l,((>y 

1 = uvil\/) £ rGj£jS = T<I) 5 T(X)=3CV) £ ERRT ( I) = TGU EbS -7 ( I ) 

I F ( mcS { ERKT ( I ) ) . G f . TulT ) -NE-c,r<=NEkR+ 1 ■ — 

440 CONTINUE 

I F ( kLkR . E u« 0 ) GJ TO 520 
GO TO 5 u 0 

: SXMP-IFIEO CmLCUuATIUN (SEE PRxtME) ' 



tail 






I A. = ^F < I J ) 

30 -oG l\J = l,UJ 
i = u >/( iV> i r (i) = 0 . 

JJ tsC Jv= 1,NV 
I-(NiHUr) V c t -*0 *T , <4 62 

462 r { I J =r ( I ) +Fi. ( I V , JV , I\) * 3C J/ ) - i oaTO-^St] 

46 ri:j=r<i)4Fj(i;tjy,iA)*3(jy) 

4 3 0 :j:;.'iNo£ 

63 ro 520 

500 CJ < f iSUE 

En'J UF ITERrtTiUN LOOP. IF Hi JET HERE, W£ JIO 
The Flli. NUH^c.r. OF ITERATIOHS HiTriuOT GOfWERGE.NCE. 

IF ( ITMA< .GT .1) <ERR = KE<R+1 

52 0 C 3 M i I N J E 

— “ — hE A T i. N G OR COOlING OF CON'TRUw iNOOE 
QJIi»=C._£ i r ( .<iC . tU . 0 } GO TO 560 

333'tuiT = l»NT . 

I - - i (IT) 

54 0 03 if. = JGIN«-CONu(iC , I ) •» ( r (IC) -T (in 

03l;(=0CIN-3F ( :C) ♦ T ( 10) - s (rG) 

iF(AlL«EU.2) GO TO a 0 0 

5o3 COHflN'OE 

IF ( <3Etr .EQ .0 . AOO. I 230. EQ .2 ) GO TO 800 
CmLu control S u£uT=FRAC*GE.T1 

IF((MjELT.Nc»C) go TOeOO - -- . _ . 

3ELr=j£_Tl 5 IF(KI3HNG. EQ.O) GO TO 800 
GO TO 280 



CA.CJlATIONS for OoTPJT, SEXr TIHE STEP LOOP 

30 0 CO.'iiiNOE 

CMi_syU_Mr£ FuRWA^^D Ai'(3 REVERSE FwOWS ALONG EACH CCNiNECTION 
30 oOa J^I^NCuhN 

0..120 = Q: 12 ( J) 5 Il = IiC0N<3) S I2=I2C0N(J) 

3C 1 2 (J)=C CM 0(11, 12) ♦(T(ll)-r (12)) 

IF ( NTINE .EO . 1 ) GO TO aOe 

JhLu avjNlnT (U\/120,3 g 12( J) ,FAv^i,'3Ei_T ,0oQul.21t0S0C122) 

303121 ( J ) =S QC121 ( 3) +0S0C121 . 

S33122( J) =SUC122 ( J) +3CUC122 
303 CO.NilNUE 



COHhUTE COOlING 0EGREE“TIH£ 
iF(,s^jACL.iNE.G) GO TO 806 

■■ - U30h = 0E.T''( .5*(r (lG') + rO(IO-n--rCMA?rlD— f-TF-tOO-JH-TirT-OT-)— OO-OH^Or 

3jH=3CH40ijOH 

306 C3NfiNUE . 

C03PUTE HEAT FLOW IN,0JT OF EACH NOCE, OF AC 0 ( I ) , QCON ( I ) 

30 aOO I T = 1 ,NT . . - .... ... 

X=wf ( IT) 

i*5RCO = QSRC ( I) S QCON3 = QCON(I) 

US,\C ( I) =S ( I) +OP( i) (1/ 

aC ON ( I ) = 0 . 

li3620JT = 1,N7 . ._ . _ 

j=L r ( JT) 

320 aC JH(i)=QLON(i ) J-COND( I ,J)*(T(3)-T(I)) 

tjE N~1.0 6 Ii"(OPh(I).GT.O.) OctN — i^PtH(i) — - — ..... ,. — 

^FLCj(I> = (uSAC(in-uC3'-i(I)) /OEN 
TO (i ) =7 ( i) 



xF ( HTIME.EO. 1 ) GOTO 300 

InTcGRATE total heat IN/juT F.^O.H^SOURGE EQUATxOM^ 

oAi_L GGNIN T (iGRCCfCSRC (I) ,FmCI»0Ei_T»0 3 QS RC 1,GSQSRC2) 
oQSR: 1 (I ) = S0SR31 (I ) »-03a3RCl 
4Q3.\32 (I ) =iOSKC2 (I ) ^33Qa<C2 
I.(TlG\AIE hEAT IN/JUT CO" TO rLOW ACROSS COMNECTIONS 

oL^L SGf'JiNT(vjiC3HC,0./UN(I) ,^*mox,l)c.lT,3SJCO<'i1,OSQv^ON2) 

jOCCiN 1 ( I i =SQC UN 1 ( i ) <• 3SQ3 ON 1 

aQL3'i2(I)=S 0:>0 n2 (I) ►0SO3j.n2 

9 0 0 C 0 . J f i N U E 

IF (.iTiNE.Ea. 1) GO TO 910 

:>3H nE.,71NG, CCO.InG -- UV 3t VE., T I L A T lOM, UA 9T AUXILIARY SYSTcH 

QCImT = 3 ElT* (FACi -^QC IN + ( 1. -F^Ci )*-GL I NO) 

Ir(i\I3HNG.£O.0»3K.KCELl .ME.O) G3 TO 402 
' ••• ur ■^HL = XLf1Ai N S iF(KIC.Eu.O) uF’<ml= FhAC 8~QG IN T = 0C INT *UF.RAC 

9 0 2 C’JN r i NU E 

IF ( uCiNT •GT .O^G^A.NO.TA'-IJ. Gf . Q VHT^Q ;hT^QCINT 

ir ( uCiN r .lT .0 . 0 . AmU .T^>H 3 .lE . TC JD - ) U V C'_ = 0 vC L >• oC IN T 

IF (uC.Nr .GT .0 . 3. AuU .] AHd. wE. 7 o.-liN) U«HT = OAHT KJCINT 

— 



c 

c 



P;<I.NT HjUr.LI' Fi.3W VALUES 
Xr^*P*vISTH*Ei.i#u) oO rO *910 
Ir (.\W^L. .cJ.2) 30 TO 9’35 

p-:iiur ■» j 0 : , r ime, i ( I ) , s iJWA) . 3 c jwa VI) , ( u:i 2 ( j) , j=i , j) , 

- 1 i r (I) , 1 = 1 , 6) 

GO r 0 "j 1 c 

90 3 H.EIr.I ttOOi, TIP.E, 5 (1) ,S { JWA) , 5 ( JWm vi) , aC 1 2 ( 1) , UC 12 ( 2) , 

1 ( aol2 ( J ) , J-7, U J , ( T( I) ,i = l, 11) 

9100CIi)u = OClN 

IF (.CCiNo.EQ.a ) GOTO G2Q 
J£_T= vciAiN*UEuI 1 S GOTO 2d0 

— 92 0 CO.MFIG'UE 

or r)u = j r r)o vr (d -r ams 

c 

C EN'J GF TIME STEP LOOP- 00 NEaT TIHE INCRE.-iE.NT NOW. 

1000 GONriNUE 

□ TWu = urwC/NTlMEl - 

C P<Ili7 totals F0\ THIS DAY 

IF ( *P^.lNTu. EG. 0) -GC TO 1100 

P^I hT 3 J Od, OATE, UAY ,H0, Q^HT , QVCu.UaHTiOACu., C0,QSP , OECtOTWO 

IF_ (li.'iOPK. £0. 0 ) GOTO 110 0 

P E i f< T 3 0 u fc 
uO 1020 IT=1,NT 

I = ur (IT) - - - - 

102 0 P^.IuT 30 07 »I , SOSRCl (I ) , SaSR02 ( I) , SGCONl ( I ) , SOGON2 (I) 

^110 0 GONflNUE ___ 

G — SUMS FuR SUHMAkY PRINT 

H=iuNTrt 5 HOY(H)=MO S OSUMi ( H , 1) = JO S OSU.-I,-) ( M , 2 ) = 0 . 
aSUHH (H, 2 ) =u 5UHH ( N, 2) vSOC 121 ( 1) vSuC 122 ( 1) VSQ0121 ( JWA VI) V 
1 SJolZZi JWAvl) 

QoU.iM (H , 3 ) =Q7HT vQhHT 5 QSUH 3 ( H , A) =Q VCu v GAOL 
aiUHH(M,E)=3. 

00 1130 J=1,NJS 

I=ISO- ( J ) 

1130 QSJ:IH(M, 5)=GSUH.'l (H,5) vSQSkOI (I) _ 

Ui.J,lh(H,&)=SOwi21(JC?) '■ ■■■ ■■■ 

IF ( JC.'. . GT .0 ) UoUMH ( H, 6) =QSUHH Ui , 6 ) v$ QC 1 21 ( JO W ) 

Q5JMH ( H , 7) =aAGu S IF ( <AlU Gl . EQ. 0) aSUH.'i ( M , 7 ) =OQH - - 

6) =UISUH 3 USUIHCH, 9) =JH:>UH 
IF ( iPRlNTH. EQ. 0. ) GO TO 1200 

PR I »T 30 17 , f1 » M'J, Q VHT , Q7uL » 0 AH r'rGAGtrflOO 

PRi mT 300b 

GO llcO IT=1,NT _ _ 

I = i. T ( i T ) 

1160 PRI.iT 30 07,I,SOSRC1(I) ,SQSR02(i) ,SaCONl(I) ,SUCON2 (I) 

^1200 CJNriNUE 

0 - — -SUHHmRY PKINT-UF- HCNTHLT-SUHS 

N1S=3 _ _ _ _ _ 

GO 1210 K=1,NQS 

SJH I G r ( <) =G3UHH lNHO,K) 

SJNH1 = JSUH.-l ( 1 , <) - — -- 

OJ 1210 M=2,NNJ 

SJ.-1fi2 = 0iUHH iM, K) 5 QiU H M ( H , K) = S JM.t 2 -SU HM 1 S S'OH HI =SJ MH 2 

-■ IF { Ur.I r 3 .£Q .2 . A'lO.K.GT . ir-QSUM.'t ,K)-=aSU«M (?t, iO/lOOO .' 

1210 CJNiINJE 

IFtNMjAOL.EO.O) PRINT 3020 

PRINT 3 0 1 9 t M T I HE » N S T E P » N C A 0 , K E RR 

JO 1220 ■H=l ,(iHO . . 

H^ = .'iOr ( H ) 

IF < uE-H -1 (H, b) . lE . C . ) QSU.HH ( 9 ,5) = 1. 

• PJTnCL= 100 . ^( 1 .-QSUHH (H ,3 ) /uSUHi'HH, e )-) 

1220 PRIuT 30 21 ,.1L t 10SUH.H( H, L) »L= 1 .NJo) , PCT30L 
II" (oU''’‘iTjT Id) .Lt-.C.) SU9TDT (o)-l. 

PJTdC. = 1CO.*(1.-Su.HTOT(5)/SJiHTOT(d)) 

P,^Il-.T 30 27, ( SUHIQT I K) , < = 1 , N OS ) , P J TS OL 



IF(IP-<.iH.EQ.O) GO TO 1300 

I PRIfJl overall SUHS of nOuE and €OftNFCTiON‘H£AT~FLOWS 
PRItH 3023,QVHT,QVCL,QAHr,O4CL,JO 
PRIi^T 300b 
JO 1230 IT=1,NT 
I = L t ( IT ) 

123 0 print 30 07 , 1 ,SQdRC1 (I ) , 3QS:R02 < I ) , SQCUNK I ) , SQJJN2 II) 

PRINT 30 2b , I J , II CJN < J ) , IZLJ'-n J) , bUC 121 ( J) , S GJ122 ( J) , J= 1, MCONM) 



'130 3 Cu.SriSUE 
oj ro 1700 

1675 ^ 0 0 5 » M G » 0 A Y I I T I HE » K i C » F rtA C » f\t HA I N 

^1/00 CGNIINGE ^ 

3---"FO.-v.HmT5'*** — 

'2690 F0^,'i».T ( /•PA30LE - PASSIVE SOuAri EnE PG Y ■» , 3X A 1 0 , 3X A 1 0 ) 

3001 Fj-^riATl/^ £533 NO GONVEn.G£NGE A7 iTlH£-=*I3»'^« GAY=*I2j 
1 * , Hv-Nrh = *.2,5«>*NEr.^=*-I3/5'.*£.-<rsT < i) =*/( 10£12.‘+) ) 

30 02 FO'iHAI(//5X’-TiH3*,vA’-OAY»,3x*,1C*,5x*T(IC)*,7X‘‘UCiN*, 7X'^QSRG ( 10) * 

1 ♦Aa*OOON(IC 3 *«2X-FFAU>/E12VA,2I6, AE12."-t7F6.3) 

30 03 FJ<HA7l/AA*x*,3x>»KF*,5»:*r*,l Ia^Gj.Xu * , iA^UCJN *,dA*S*) 

3J0>t Fj-sHkM 216, h£12. 4) 

30 0? FO X.1 AT (/2a* Gate*, 2a *0 AY*, 3X» HO*, 5a ^QOHT’tdA^'JOOL*- ,6A*QAHT*, 

1 6A*G‘<Cw*,3A*uJ*/io,215,5£12.‘+/5A*06P*,9X*jE0*,9A*0rV(j*/3El2.4) 

3006 FJXHKr t/‘+A*I*,5A*b0i,RGl*,oX*300Ho2-^ ,6X*3QC0Nt*,oA*S0CGN2*) - — - 

30 0 7 FoxHh T ( i 5 , •+£ 12 . A ) 

3013 FJR,-1«I(/3a*:T*,tX*I*,5a*T*,1 1a*CPH* ,9X *3* , 11X*3P* , 10 X*SC0M* ) 

•3015 FORMAT ( 215, 5E12. 4 ) 

3017 FjR.NAi(/* HO.NTH i NO Ex = * I 3 , 5 X * HON T M OF Y£AP, = *I3/ 

1 5A*uVrir*,OA*UVCL.*, 3x*QAHr*,dA*CACi_*,3X*00*/5El2.4) 

3019 FOXHATi/* 3UHHA\Y -- H f i i t = * 1 5 , 5 a * I, STEP-* 1 5 , 5 A* NC ALO = * 15 , 

A 5A*KcRX = *i0//3A*HJ*,2‘,*0£G uaY*, 

1 3A*QjUT*,3A*Qrf£A7*, 3X* QOOO . * , 2 a + a3 0LAR* , 3X *0_0 AO* , 4 X*QA C*. * , -- * 

2 3A*G3I MC*, 2a*PJTSOl*) 

3021 FOxMmT ( 1 5 ,dFd. 0, Fo. 2) 

3323 FGRHArc//* TOTALS*, 5x*JVMT *TO A'*Q V G-C-*TTi X * Li AH T *7^ X*-irAO t. *T 3 X*"JO*/ 

1 7a,5E12.h) 

3025 FJRHmT ( / 4 a* J*, 3A*I 1 *, 3 < *I 2* , 5 a*SUG 1 21*, 6A*SQC12 2*/ (315 , 2E12.4) > 

30 2 7 FjRrlA7(/* S UM S* , 6 F 6 . 0 , - d . 2 ) 

3 0 2 0 FJxHhI(/* no AU«.iuIARY COOulNG - UACL IS INT ( T ( IC ) -T CMAX ) OT* ) 

•,0 00 FORriAr("4,0,Fb.l,10F:?.l,6F5.Li) 

40 01 Fjri.lMr ( F4 .0 ,Fb . 1 , 11F5 . 1 , 11F5 . 0) 

40 0 2 FJ.a«AT(* TiHE 5W SGI SGO GCl QC2 QG7 QCd 009 QCIO*, 

1 * LiCll G012 QC13 ri T2 T3 T4 T5 T6 T7 73*-, 

2 3X*T:, 710 Til*) ' 

‘♦OOh FJ.\.iAI(* TIHE iW SGI SGO QCl QC2 UC 3 QO-+ QC5 . 006* , 

1 * Q07 QCd ri 72 73 74 75 76*) 

-+003 FOHHAT(/F DATE=/I7) 

4005 FOR.-lA7(//* S3S3 TiHE STEP SuaOis/IOED 5 Tf\I£S - PROBLEM*, — 

1 * oTaPPEO ££SS*//6A*M0*,dA*L.AY*, 6A*ITIHE*,bX*.<IC*, 

.. ^,^X*FXAO*,s.X*r<EHAIN*/5I1 0,2 F10.4) 

£Nj" ^ 



F 0 R W4 \ U _ A L40 _ RE'O t iaS £ _F-_ 5 H _ INT EGRA T 10 n_RO JT IN E 



30 3 f\0 JTiNE 3 G('ilN7(Yl,Y2, FAC, 0Ei.TA,AF!_U3,Arl I NOS) - 

j Y1,Y2 ARE rtElGHTS OF h FU NOTION* FAC IS THE 

3 lt(7cGb\“,TxL»ii A0c.RAG.fcNG FACTOR (i(Of\HAcLT *5 FOR 

3 ■■ ■■ r.HIo PROGRAM-- I.E., 7-1 E 7 R A PE 70 1 u A L KUwE; 

3 Y1 iS CJ).SlJEi\£0 lEFT Or Y2 (Trie rOeO^' OAlOc FOR 

3 TiHc'nlSE ^INTEGRA r ION) AND F'nC Ib rriE^'WEIGHT OF THE 

3 NEW v'lfcuJE* uEuFA IS THE HORIZuNThw QISTANCE OcTAEEN 

3 THcTi^iJ PCI NTb, 

3 hP_u5 is TmE FOSITIOE A, REA FOn, ThcSE hOINTS, 

3 AHIiiUS IS THE NcGATIi/E AREA OcTAEcu THESE TWO POINTS. 

% 

QoA=( 1.-FA(3)*Y1 ■£ QQ3 = -AC*r2 

IF (Y2*Y1.LT.Q.) GOTO 100 

s 

3 OOTh Y2,Y1 ARE THE SAME SIGN 

Ah1=CclT A* ( uQdfOUA) S uA2=0. 

GO 7 0 2 0 0 

3 Y2, Yt ARE OF OPPOSITE SIGN. GET TriE AREA OF THE 2 TRIANGLES 

100 AA1 = LEL r A*UUA* *2/ lUOA-QUB) 

A42=0cLi A*0Q6**2/ ( U03-QQ4) 

3 SET T.hE PoSITiYE AnJ NEGATI/E PARTS INTO A1 ANO A2 

20 0 A P L 0 S - A 1 A A 1 ( A A 1 , A A 2 ) 

AHl,,'Jb= AMIN i ( hAI , AA2) 
ixE T uRN 




Mu i §^N A T E_ y A 7 E _ b(J 3RGU T 1 NE 



3U-3<\OuriNE UAVEl (mRG) 
GAlu C-MTEiARG) 

Rt r URN 



C C Jfi 1 SU jKGU n NE 

CE====£:£r=i::5i:rf=z==r===rr====f£f=:E5fr£f=r3=£rE=£ 

c , 

sjSfvojr k^^E co^^^KOL 

c j I M c \ / j L M g 1 / : . V , F , N- r, N M _;yr< r r 5 tr rrir; L 7 F c 5 a r,-zFx 9 cr jte t 1 5 tm 

t Ki.i-»n\^tKxC fXr 'l^X»N£rv.“\fiJ*’iIT otKAuA'Ct- tKoAwO»NJc.wT »KCO.'iU 
C0'i.Jo% / AuAj(X/f<rj.MEfMr*HEf^rLHc.?jMi E*‘JMf»LlHirifjAf2|ti-*At (X2) > — 

CJi.nON/3LH2.i/T ( 3 O) »T0{5Q( tCPr){£0) »S(:jG) »SP(30J »CuriG(30»50) t 
1 SGofi{3Q)»0FAo0(30J»c.a\r(::0)jFiJ(£Q»5Q»2)»FG(:>j»30t2) 
COM.iOX/3uAJh/«(pC,50) ,3(5-3) ,GuJaC1(50) ,34 -.'.02 (5 0) , 3 4C0N1 ( 50 ) , 

1. 34uON2(5(i),Uvrii ,U'/CL,3AHT,4>-oc-,4HZ,U1.4w,Qr.’^AN 
v4*iMC),/ Owk^up / oU**lxr ,Q^k4 ,4Cx'(G,F #,vj,F 
1 PX ,r*.Mz.X,GtiGr', , TGm^E, OJ ,7 MAa, i.iXN ,T4A:'^,FACX,\Ei‘iAIi4 
C'3 ‘iMv^S/3LAuo/>.3Xo (53), J30N (50), 43 12(133), sac 121(103), 

1 S4ol22(lGu) i C-iXNu,rCiMX4,roMAAJ,rCC0u 

GJ.'lii'j', / 3LMjF/T?c.f\0AY,33CoN,rA*lS,UvjiNT 

GvJHiloN/ SuAJo / xl)AY , lL)PE"kY£, XNl)m Y,KUATA,QrlC0iN'7 ,AlA7 ,3»CN&, • 

1 fSkitUEA, 1 Qu 1 ,LfE>./ foG^JE^ ,SX'*4 uc.\^,SxN*^T ,30SwAT,asP 

2 , rl rs.3 \ iT ,uICC'T ,>^_riUi, o3L,.'i0oril,M 03 H 2 

Coi.iCr,/ 3LA0-J/ TA j (50) , UEuj (53 ) , Qmo ( 5 0 ) , RrtO ( 5 0)-rr(a3-,-ISC;rttr5 

1 riwT(15),WAXl((l5) ,Ab.Z(ip),;4G^Z(15),NJA,XAMi(10), 

2- HLFA(15),OriKnS(15)»OS£?\(lp) ,TGLZ(15),JA53(15)»1 SO.. X ( 2i+ ) 

COi.'io.Y/SLAOlO/ajNMAA, xDXFF (5 ) ,E \ ,xH0,T4cP. C, 

1 NCuNN , i ICON ( 133 ) , I 2 CuN ( 100 ) , U( 130 ) , AC J.4 ( 1 30 ) , UO ( 103 ) ,_UN { 130 ) 
oOi^JNON/ 3l.a 011/ rLLTC,MZXi^li-/,(^H0O,w034 L»,oXNt o,F ^^\/,F •^k,'j,-\c3NX, ^ 

1 JGP,KNa 

CJMMGN/r FsOMSE / VJlF ,PCUN (53 ) ,i<3EnT , XAXk, JAXPv.,f<wNGT H,phOSR , 

1 AsPxaT,FG3R,JWA,TRCGaT ,4CW,KwAlt C ’ 

CO. "i i’l o N / kO F lGS/XS.x CP K, -01, xO 2 ,X.^\lNTH,IP.RIi'>r0,IPP.Ii\T'i,XPR3M,<H£0PR 

G J I*! -ION/ r X flE i / T N X Gh T , t .iO xN , T £V £N , J E L T 1 , 0£L T , T X •'■)£_ „ 

C 

C 

^ XNTEGER CAY , UAYl , 0AY2 ,YR, CA T£ , units 

C TESTi FOR TE.-IPERh TURE LXHITS ON CONTROL NOOE 

c <XG -1 FOR r(ic)<TCNiN wi f HO j r-uGiN r-a-GiN> 3 

C RIG 0 FOR TG;iXN<r (IG) -fTGMAx WITHOUT QCXN; 001.4=0 

C KXG 1 FOx r(lC)>TC.NAA wiT.HGJT OCIn; QCIN<0 

XF (a.xGHNG,£u. 0) F.£.'1AIN=i. 

F.RAo=REHAXN 

GO to (100, 200,300) ,K 
100 IF (UCIN.GE. 3. ) oG TO dOO 

XF(X^EiNi •Ea,*4*AN13,V0LF,Gr"*0'*)'” C'3 T 3 S 0 0 
^^\^,o = GG-^(U/ (CulNu^UClN) S O^x.N — G# S nXC=0 

KF(IC)=1 S ,iF=NF-l i N7=N7H 5 GO TO 503 

20 0 IF ( . ( IG) .lT .TCNXN ) GO TO 220 

IF ( T (xG) .GT .TC.RmX ) GO TO 2.+ 0 

GO TO 60 0 -- - - 

220 FxAu=(TGi'lXN-TO(IG) ) /(T( IG)-T0 (IC> ) S T ( IC ) =TC .HI 4 
KIC = -1 S GO TO 2t)0 

"'25 0 F<a u= { T GMAx-T C ( IC ) ) 7 ( T ( IC) - T-0( iOF ) — S -T-R'IC ) = TCrtA-x— 3-<'IC=-l 

2o0 KF(xC)=2 £ NF=N-+i £ N3=N7-1 S bO TO 500 

C 

300 IF (oCIN.lE. 0 . ) SO TO 300' 

F\Ao=UCINC/ (GolNO-QCIN) S QGI4=Q . S KIC=0 

<F(xC)=l £ NF = 4F-1 S (,V=.N\/+1 ■ -- 

C 

500 CuNTlNJE 

■ IF ( r RaC . lE . 0 . 0 . OR .F RAG . Gf-.-l j O )- FRAC^ . 0 1 " 

FRrtu=FRAC'''(xEHAlS S h.E.MAiN = REMAlN~F.RAC 

KICh4u=1 S xF I <0ELT ,EQ. 0) RcTURN 

uO b 03 IT = l, h'T 

x=i-T(IT) S Ih(I«£G«IG) go TO bOO 
T ( I ) =10 ( I) +F.R«CM T ( I) -r J( I) ) - - 

600 GoNlINUE 
kE T uxN 

dOO C-JNiiNUE 

.<! G .H.Nij= 0 S RETURN 



SOO.xUUtiNc 140 AT A 
C*3ASE=10u. 
b ^ I'i J C Ou — C. o T 
C ♦EaGEPT 

C*INOEN'f = 3. 

C^NULlbT 
G »STAT = 10 . 

C£ £=£££££-£££££££££££££££££££ £££££££££££E£55-£===p= ^ 

CoRi J0N7 5 l Au 1 7 .W , NF ♦ (NT , NMh X > .<F (5 0 ) , i o t w ^ ( li 0 ) ? uF (5 0 ) » LT { 5 3 ) > ^ 

l?XIofir'JAAjt<FK.<>UMxTi>tK>>JA^L»KUMu.OfX->t^L-3 >KCJ»iU 
(jO.i.iUN / 3 l a 02 / .][ I )iF » M r I ilE » I r i HE f J M r E > J A Y ; J A Y I « G A Y 2 f »MDh y ( i 2 ) ) 



q L. 



i MO , i10l , YR, MMU , N3>; Y 1, I'J JA TE 

0 JMMO.S/ 3LAO J/ T I ) , 70 ( po ) , ^PM <3a ) , 3 ( aO) , <5u) , ON J (3 0 , 5 0 ) , SOON 

1(5 0 J *v«hwuo(30) iE,-\r<T (pJJ ♦rN(50»5Q»0) tP0(pu«3 0>0J 

GJMMON/OumJw/ a I p 0 , ?G) , J ( -30 ) , i ( 5 0 ) , S OS RoO ( 5 0 ) , S GOON 1 ( 5 0 ) , 

IS'.Jvji'tOlpU) *OVf(T»OV'^L*G-»Hr * lAot_t''«fiZ»OXNO»Or-<Af'i 

^^OMrlON/pwrti^p/ jOMRi t*\M G i P ^ t 7 x a t 

10wjiV"j| I SASEfuJ* I MAXtTHx'J *7 cMxIiS 

u JM.nON/ jlA J o/ OpRO( 50) . aOGN(pJ) .OOlOdQJJ .300l2i(iQJ ) ,3QC120(10Q) 

!♦ 7 s/TliUj i oMmA t7oM*^^0f70'ix'iiJt I oMHAiij 7 loOOt. 

0 J MMOiN/ p (_ A 0 7 / 7 Pi In Jm Y , S3 GOiN . r M i 3 , uG i N7 

bJMriON/pLAu'P/ aOhY t iuPt“^Y \ « iNlImY t RO«T «t QHi^ON V tAi.AT *0_ GNG» RING ca t 

liJ^l ijEa »u0S0iu»Sir^0EC»SIr^LH7 toooLAT »OSPfnRGKx7 >QiOJT > 

2R-»Ho7 tisOO^u tMjOnl jMu5H2 

OOMMCN/pLAOp/ rAO(pG) ,VEOCrtpO)VJHO(5arTTHOfpGrTr'rJOTtOOUlT5>TTru7 

idpj f.iM/iMdP) ,mGw 2 (15) »NGi_2 (15) ,.dA,IAM3(ia) ,muFh( 15) tOHANG(lp) » 
20p£?Rd5) ,7oL.Zd5) ,OAoS (15) ,iSOwA(24) 

OOHMuS/pLmOIO/ 00iO3Aai\DIFF(p )>EAt.\HO»7Z£.'; 0»NOONN,I1CUN(1QQ) i 
tl2 0u,N ( 10 0 ) ,0(100) , mGON ( lO'j ) , Ou dOO ) , UN (10 0 ) 

GOMi*iGN/ SLAOII/ 7 xu7 0,mZx3u, xM Oo,OOSTo,SXM i v^jFSS'j, Fov?C,^£sNX ,J*vP , 

KNO 

OJMrlON/7F.OMJ£/ VOLF , ROOM ( p U ) , K V EN 7 , I A XR, JA XR , R_fiG7 H , RHOSP, A S PRA 7 , 
lFOSf^,JHA,7RCGMr,JGR,KWALL 

CJ1MGN/ICFLGS/;SF,GPn.,:0l,IO2,..PRIi-.7ri, IPRI n7 0 , XPRI m 7 ,i , X PR S K , <H£OPR 

C J NM C N/ 7 IME S/ I NX GhF , 7 MO RN , 7 £ V tf( , U£ L 7 i » JcL T , TlilE 



OXMcNSIJN EPSGLdQ), EPSGO(IO) — 

lN7£GcR NAmARAY(3Q) 

XN7c.GiR DmY ,UMYl,JAY2,Yr^.,CAi£jUNX iS 
0 m7h SEG7ION 

OATh (NAMAR^ Y ( I) , X= 1, 2) / fG0()7R0w TMPF, /riOURLYIO*/ 

0a7m (NAMH.*\MY(x),i — 3,*^) /^S7Mr\iw'Mrc.^,^h.ATx7ULiF/ 

0A7m (NAMARmY(X) ,X-p,6)/F0i.0N&I7Ou£Ft?RiSiNXGH7F/ 
Oa7m ()(hHhRAY (x) , X — 7 ,o) / ^OLOmoF, AmUaCOOuXNGF/ 

Oa7m iNhMARAY ( i) , 1=9, 10) //VEN7SF, /MON7HSF/ 

0A7A NUMNAMO /IC/ 



OEFAUL7 values OF PROGRAM VARXA3LES 

Oa7a ALA7/35.1/ 

OATh IiS0A7E/Cu1373/ — - - 

Om(m IGl,X02/ii,C/ 

0«Tm XPRIN70/1/ 

OATm iPRlN7H/G/’ 

OA7a 1PRXN7M/1/ 

OmTa XPRSM/i/ 

0 A 7 « iS \GP k/ 0/ 

Oh7m KmJXLl/0/ 

Oh7h MiiDFR/1/ ■ 

J*,7h .<VENI/0/ 

OATm RWALL/2/ 

■■■ OA7m .'iOAYl/10/ 

Oa 7m NMO/1/ 

0A7APESNI/J./ 

OAT M ^C■^Xi’i^^, TtjMAXf), TGMING, 7 CMaxO/6p •, 75* , op • ,7p •/ 
OATh T£ VEN/ 17./ 

0 A 7 M r M w K N / 7 . / - ■ ■ — ■■ ■■ — ^ 

OATh 7iNiGHT/23./ 
uATm GL0Al)/1./ 

•JNX|'S =1 i NriM£=2^ Sir MAX-5 

lOPtL- YR = 3o5 5 <OArA = 0 S RIN0 £a=1.p 26 5 £X=0.p-^l2, 

R HO — 0.3 S FmC— p.P £ r f-»^/x — ').p S rOuT — 1.0 

rXu7C-0. S mXIM0=0. £ RhOC=0.3 

<GA^C = 1 i ,<G£w7=l S <C0NU = 0 

QXCU7--1.E6 £ MpSril=13 S MJ3H2=0 

RH0j\ = R..'.G7H = FCSR=0. 

Om7m (RJIFF(N) ,'J = 1,p) / • lo, •2‘», .2p,. 33, «3p/ 

Oai'a is jla/ 2->-^0y 

00 ilO x=1,NMmX 

uO 110 J=1,NMAX 
GGNU(X»J)=fl. 
liO CONTINUE 



REAO IN VALUES FROM FILE FOR PROGRAM CONT^Ou 



■•REWiMO-2 

120 continue 

rEAJ (2,130) NAME 
IF (EGp(2).NE.0) go 70 130 
130 FORMmI (AiO) 

00 140 ^ OuP= 1 , NJHNAMS 

■ LM .. P PAV/.1. 






r : O 






1 •+ 0 C G .'(I i N O E 

P'tiNl ^fXii'HPur ^G.'IMANO not f^ECOGf^IZEO^tMAit 
S r •Jr' 

^150 GJTj (l30,17G,ir3fl73*ir9fl61,i62,163,loV,l65),GU0P 

if 

c R£Aj coNKOG f,GjE T£H?£-^Arj<e LiNrrs 

IbO f\ciMvj (2i*J rCMx*'Mf*G''^^^NfTwMiNJ>rG*'iMXO 

GJ TO 120 

G t\£ALi X ^AwJE Gf iNSUwATION J.N TrtEkHAL WALu 

Ibl i\cAu(2»^) i\ESmI 
GuTJ 120 
C 

C ‘^EAo VA^ue FCP 'LOAO/AKE-A 

lo2 xEAj (2,*) ULOAj 

GOTu 123 

G auxiliary COUuING specifieo. set flag to 1 

163 XAGaGl.= 1 — 

G G r u 12 0 

c 

C reau vent OUNTKOL-FLAG 

16L REAU (2,*) KVENF 

GJT^j 120 

C 

C .^EAj NG’-tSErv, OF MONTHS “Ox SIMULATION 

165 REAu 12 ,*) NMO ' 

NOAri=300 
GuTu 123 

c - 

C REAu OATES TO START A(0 5TO? HOORuY PRINTOUT 

170 REAuCa,*^) 101,102 

GO lO 120 
C 

0 REAu STARTING DATE OF SIMULATION - - ■ 

173 REAu(2,*) IN0«T£ 

GOTu 120 

C REAj uATITuOE OF SiMUuATION 

176 REAu (2, «^) AlAT .... 

GOTO 1 2 '3 
C 

C READ THE nEuATIUE lONGITUDE - 

179 r\EAu(2,*) QLCNG 

GOTO 120 

C GO REST jF set-up FOR THIS ROUTINE 

160 CONTINUE 

C 



190 
C 



200 
C 



210 



GO TO (190,200), UNITS 
CON riNUE 

•UNITj 17. 3TU ,HR,FT ,GEG-- 

S3JuN=1.7132£-9 5 TPERjAY=24. S T3ASE=o6. S QHC0NV=.3172 

TZErO = uoO. S QONMAX = 317.2 

G u TO 210 
CUNT INUE 

•UNiTo 27 S.I. UNITS ( J , 3 , .M , 0 EG -G) . . 

S 3 G uN = 6 . o 7 0 £-0 S T PER J A Y = 6 6-, 0 0 • £ T3ASE=16. S OriCON'y = l» 
rZ£.v.0 = 273. 5 UONMAA = iaOC« 

'CONJ i.NUE - - 



— STmRT SPECIFIC .MCOEL -- SOUTH MASS WAuL 



1) - PEk UNIT A.-lEA OF Gl 



TRJMuE nALv. (KV^ALL 2) OR AATuR NAul (K'WAlu 
NS Gw ,N,./N-ZER0 FO-^ massive EATER I JR WAllS 
SET riuuE.. PARAMETERS 

N G L - — N J li o E r, OF P R I M A ,R Y S 0 u h R SOURCE ( M ■) S S W A u L ) GLAZINGS 
•HoHT-- sMw htiuMT FOR T ME r M JUi RCu u a T I ON^C ai_ C (FT) 

U^OmC-- STATIC ROuM McATIfiG uUhJ CUE FF IC I ciN T PER U/il T SMR GuAZING 
1 T U/ iMK/ 0 EG • F /S Q * F T ) ■• THI^ uOES .-tOT iNGwUCE HEAT uOSS THRU 

MASSivE WAuLS* 

P.E SISTAi-iCE uF 5'IH NIGHT INSULATION (0£G.F/(3TU/HR/SJtFT)) 
THEkMCCi R-^ULAT I ON VEuT G.ONTi%Ou FLAG (SEE PROP) 

SfiW Tu RuOH U (OTU/ir. /0£u«F/u>3«rT) 

AIR 'VOLUMETRiC mEmT uAPA..#ITY (3T‘J/0£o*F/CU«Fr) 

iClEeERATIun uUc to j^aVITY (FT/S/S) 

■ aIr thermal CjNJu: T iVir r ( jTu/H^/ JEG .F /F T) 
tilt cf SMrt Glazing ok of smw if 3are (oeg.) 

iLL AZj-MUTh (OEu.) 
iLL AduORdTiVITY 

,AZiNGTMlCRr^ES5PERi-AY£R(FT)_ 

GLAZif^G THErHau GUNuCuriVITY ( 3 TU/HR/oO .F T/FT) ' 



G 


_( 


G 


3 i *1 u R 


C 


Rc S N I • * 


C 




G 


J J r * • 


G 


C A I 


. % 
U 




0 


r fi c M I \ ~ 


c 


t i L i * - 


c 




c 


H F - W’ 


c 


roL-- G 


c 


T H vj L — 



i 



25 0 



AI ^uA=>-- 
EP 5 H , EA i 
Tr<Cu>* r — 
F X H o .*. I 
-1 

- 

A S P N M T 
FCi.x - 
CJIo-- V 

M 7 J H S “ ~ 

G u M Z i N 
ASOmG-- 



GLAZ ING 
Or i iS Nc. P 



GPACiiWj J£TrtEt.^^ GLAZiiW^ji) (FT) 

A,EPSG-- I.P. £ 'IHiiS ^w'lry UF HALL, S<<y,_«NJ 
T •<>. Ni '1 i, o j I V I T y Or JuhT.iiG OiG IiGnEk. oUprAi.c 
ZCGTau jPEGUwAR r.c.F_EJ7uP, '/Er-iTICAL W«l.l (5£E SUNSKC) 

~ kEFl« LEuorn/ Gi-hZ» nExGrtT 
F£FLLGrj./xTT ^ 

~ »£NG 7 i/WI 0 TH 
HMLL/f\£FLc.CTGR 7 I£i'l £>-( 070 .% 

EK7 OiiChM%GE OOEFF. 

AREA uF ONE ROW (TOP UR 60770ii) OF J£r(7i PER UNI7 SMH 
Af^E A 

GL • AP. EA 



LAZ 



ONE 

'W A (. L / G L M Z ING 



AIR 



UR dOTTOiD OF J£r( 7 i PER 
SPACE FlOH area per UNIT SRW 



220 



230 



24 0 



,NGL = 2 
N S G w = 
USTrvN 
Ti>TG 
r R 0 G L 
A £ RG A 
EPSn = 
RwNgT 

FOG (\- 
G J I o = 
00 22 
EPGoi 



5 £ CAIR= 0 . 01 a S GRAV= 32.2 £ rHCAIR=Q .015 
. 5 AZH= 0 . S AwF= 1.00 S 7 G L = 0 . 12 5 / 1 2 . 



S HGHT= 3 . 
C 

= 1 . 

= 90 
= 0 . 
r = 0 

Q . 9 

H = 1 

. 5 * 

0.6 



ij j. 

(13 



,25/12. i IFCNGw.EG. l)AIRoAP=0. 

5 £PSG=Q.9 5 EPS«=0.9 S 7RCOAT=1.0 
.0 £ RHCSR=0.a 6 mS PRAT =5.0 

(KLiNGTH + l .-SORT < RLN G T 2 <■ 1“. ) > 

^ A 7 OAG‘“ 0 .Q 3 S >ac> 0 .(G= 0.02 

G = 1,NGl 

) =EPSbO ( IG) =EP3G 



IF (KGA7A.EQ.0) GO TO 230 
FOR owEA\-GA*' WEATHEk DhTA 
OwOIiG-0 . 

■CONTINUE 



GALCUcATION:, (SEE DAYlY) 



<0£..7 = 1 S KCONU=0 S KCAl.C = 1 
IF (KCA..G.E0,2) .<C0NU = 1 

Cm^C(J'.h7E Thermal network parameters 

GO TO ( 240,250) , KHaLL 
C 0 ij f I N J E 
WATc.% WAlL 

ASSUMED TO tiE SAME Ai<£A''AS-‘GtA‘ZING 

SInGuE NJuE (NODE 1) FOR 5MW - Gi_AZ 
NOOE 2 IS %OUM AIK -- NOOE 3 IS 



NGS START AT 

CUTSIOE MM3IENT ... 



NODE 4 la Wh i_l/ Gw A Z I No AIR S^AUE 
CP'l.-v-" aMi^ HEAT_CAP«CITANCE (STU/OEG.F 
XAlr.*** WmUL/GlAY.. AI'^. SPAUc NOOt. — ■ 
JAI%-“ 7 HEnHUU IRCULmT I 0 r( CONNECTIGn (XAIR 
IG~~ ROOM AIR f.juE (CCNTROw NoUE) 

JC? — CUNNECTiON OF ULCAO 

JWA~~ GuNuEGTION SETwEaN SM« SURFACE NOOE 

seuUEiNJE of connections th'^u Glazings, 



PER UNIT SMH GlAZ. A PEA) 



TO IC) 



AND lAIR - STARTS 



HSEo-- .(UMSER OF oERItS waLw ScGMENTS (GEE MAGONr^Y WAi_L) 

CPMR=.,? S GPM(1)=CPMR - — - 

NCOhN = NGL+o 5 IF (NGL.lE .0 ) NOON. (=3 S JO = NCONN 
N(JoE = NoL+’ 4 5 Ir(NGL.ui.0)N(UuE = 3 S IO=.NNCOc 
IAI(\=JAiK=4 S IC-2 ’S'JO P='3 S J'tf A—5~'£ nSEG ="1 
GO (0 2o0 

CON i InUE . _ . 

MaSOi'iRt HAi_u 

ASoUittO iO d£ SAME mREA AS oLAZiNG 

NSEG<-2 hAlu iNGCES x.N 5E R x E 3 _S Tm%T I NG AT_.NOuE 1. NU.OES 1 A.N-J 
i(SEo<’2 «kE .‘‘lAaSLcoS SURFACE NuJEo. NStG^a lS .ROOM AI.R (CONT.ROl.) 
NaEu+4 IS OuTSIjE A.MSIEiNT , MSEG^S IS WAll/GLAZ. AIR SPACE. 

Glazings start at nseg + 6- — 

T.HOuN-- iMASCnP.Y ThEiRMmu 0ONJUCTI7ITY ( 3 TU/H R/ OE G. F/F T ) 
rdIwK--WAuw THICK NEoS (FT) 

70loP~- MAsONRY VOlUMETRIC MEAT CAPACITY ( 8 TU/0 wG • F/ CU . F T ) 



THGu.N = 1.0 S ThICK= 13 ./ 12 . S 7 Ov.aP= 30 . £ 

C J.( i = Ns E G* ( THo ON/ ThICN ) a 0 ON 2 = 2 . ♦ C C.N 1 
CP M(\=v'UwSP*’ T hI CK a CPM R 1 = CP.'Ih/NSc G 



N 'S £ G = 4 



26 0 



NCOiiiN = NS£G^'iNG(_*'7 S IF(NGw.La.G)N(vONN = riSEG+4 S J0 =nC0J1N 
iNNOu£ = ,*iSEG + NGw S IF ( J Gl . L £ . a ) NiNU 0 E = ^( 5 E & ^4 a IG=N,t 0 CE 
XAI’\ = ns-G + 5 £ J A I .%= r(S E G 4 5 £ ix. = NSEG + 3 S JCP = i'')Sc.G 4 L 

continue 

FOR jjlH TYPES OF WALu 
JW;\ = iNsEG+o £ JCh -=0 
N JA = l S lA.M PI 1 ) = OOP 

NJS =1 S ISJL( 1)=1 

GO.Nf(cCTiJN Parameters 

ilCoN(l)=l S I 2 l 0 .n( 1 ) =I AIR^l £ ACOM( 1 ) = 1.0 
IF UjGu.LE.O) 12 C U.N ( 1 ) = I A .MS ( 1 ) 

J 2 = NSwGH' 

00 270 J= 2 ,JCP 



uN( J) “U0< J) =wONl 
17 a CONTINUE 

U.u JG^-1 ) =UO ( JC^ - I ) = JSr RN 

UN ( JCPJ =UD( JCP) =ULOAJ 

iF (Nol.Ll.J) go to 290 
ilC^N( jMlr\) = lMlF i 12C0N ( J \IF,) = IC 6 AGON ( J« IR) = 1. Q 
IiCuNtJ4i,^-H)=l S I2:,uS ( JAI G-t-U =lHiR S AGON n =1.0 

IloLyN(J..i(\>2)=lAlR 5 i2COn(JA^R+2)=IAIR*-l S AGON(JAiR*’2)-l.’3 
Jl = i^SEG^•<J 5 J2 = Ji+Nou-l 
OJ 2dC J=J1, 02 

ilOOrao)=J-2 i I2CON(J)=J-i $ ACON(J) = 1.0 

■230 "CONTINUE 

i2CuN(J2)=IAMd (1) 

290 CONTINUE 

— nooe paraheters 

00300i = i>NNOO£ - — - 

f(l)=7Q. S KF(I)=1 
300 CONTINUE 

I = iMMO ( 1 > ■ 5 -KF C ) =2- 

IF (KwAtL.Ea.l) GO TO 320 

I2 = NltG<-l _ _ 

JO AlG 1=2,12 ■' ■ 

310 CPM (iJ =CP.iRI 

320 CONTINUE - - • 



---CA^CuoAT E GUNCTANTi-FO^ - TEnPERA TUKE-OEPENOETt T-UtO— ITF-PROP- 
i\COi«(i) = C0IS^M00AG’'30Rr (GRAV'^riGrtT )*3 dOO» 

RCOn ( 2) =2.*oAIR i RCCN ( 3; = 1 ./ (3o00 .*ASOAG ) 

RoOkCt) =1. / (AiRjAP/THCAIR+Tjc/T riG GL ) 

I CONOTAnTS rOrt nAOI«TIOt< U^S 

RCO.N ( 10) =39CGN/ ( 1 ./ EPSN n ./£P:>A -1 .) 5 IF ( NG L . uE .0 ) GO T037 0 

RCJN(lC)=ioCUN/(l./EPiPt-l./EPiGI(l)-l.) 

Ir (Nou.Eu.l) Go TO 3G0 

NG = rtGL-l 

00 33o IG=1,NG 

330 RCO. N(.G<- 10) =S3C0n/( 1./EPSG0(IG) +1./EP3GI (IG-H)-l. ) 

3-+Q RCOim ING. +10 ) =SdCON/ ( l./EPSGO (N3 u ) <-1 ./EP3A-1 .) 



— SOwAx 5COKCE iiNFORNATICN 

TIuj (1) =Ti’.TG S wA2:m(1)=AZH B ALrA(i)=ALF 

Z ( 1) = 1 . 0 S N'G L Z ( 1) =S Gl' 3~f -ouZ (l')"=TGt:: 

0A3o(l)=l.-EAF(-1.215*E;^*TG.Z^l)) 

OHAiNG 1 1) =Q3 EPR{ 1) =0 . 

IF (,oG.^ .Lt.O ) GO TC 370 

— OPriwNAL EATf^A FO^ MASSIVE EaTERIOK WAwLS -- - 

N:>Gr(~“ NUfi3£ R_OF viAu.- SEGMENTS 
rlAli\ — outside air FIlM COEFF. 

U-<MW“ ROOM To wAlL~U’ < iMCN.UOES'~A.Nt — IrfTcRIOR ^IM SUcAT'iOH) 

RwIinS — RESISTANCE G” ANY ExTErIOr WAL_ IisSUwATION 
T'iCrt”- ThER.'lAi. 0 C iN J UC TI V I T Y OF MATERIAL 

Tn.<r(-- riMLu TnICKNESS 

VS°A~* WAwL MAT. VOLUMETRIC rlEAT CAPACITY 

Art-— Arch/ SM'OGwAZINGArEA 

OCW-- WAcL fU OUTSiCE 0 G.NNEO T I ON. HEAT FLOW THRU 00 W IS AQOEO 
TO That thro JCP to get total OcOAO. 



iNCOHh= JO+NSGW+3 B MiNCOE = IO +N3GW + 2 
HaIr=-. 0 

URMrt=1.5 B rwINS=0. B 0 WOU T = 1 ./( RW I NS + 1 ./ hA IR ) 
THCrt = 0.3 S THl<W=1.0 i VS^+i=30. 5 AH = 2 .a 7 
C P M rt = E rt *■ T h< W + V oP W S C P M w I = C P M A / N i G W 
GJi=foGrt*TriCrt/THKW B CC2=2.*CC1 

- -- 0=00+1 S I ICON (0 ) -IC--S I2CON (0) =10+1— 3 -ACOM (-0 ) =AW 

UO I 0) =UN ( 0) =UrMW 
OJ2 = NSGW +1 
00 390 O0=i,0o2 

j = 0u + 00 + l t> 1=10+00 

xlC0N(O)=i B I2C0N (0 ) =1+1 6 AC0N(0)=AW 

cLx=CCl i IF) 00 . £0. 1 .OR. OO.EO. 002) CCI=CC2 
UO < 0) =UN ( 0) =v^G I 

350 - ■ CUM I IMUE ■ - - 

0 = 0j + i''iSGW + 3 B I=IO + iNSGW+2 

IlCoNC0)=I $ I2GUM ( 0) =I AMa ( 1 ) B mCON(0)=AW 
UN ( 0) =U0 (0) =UWOOT B 0CW=O 

Il=xU+l B I2=I0+NSGN+2 
JO 360 1=11,12 



T (i.) =70. 5 KFd) =1 

J.F (I.Ea.il.JK.I.EJ.12) 3U TO ibO 
mF-I ( I ) =Ct^HnL 
36 G _G 0 N 1 1 Mj E 

E.JJ £-ATE\IOf^ CAlC. 

370 caNri.suE -- 

IF (<HA-L.Ea, 0 ) GO TO *10 
IF (,<G0NU.NE.1) go to +1G 

— FCJ-<_uGfjSrANT COEFFIOIEuTS U^^» UQ 

T.\Ec=7o. S TFUiNC = KHUr-N(T-\EF»TPv£F» TZ EkO ) . — 

OAi.'\=l. S UuiKC=0. S niI\ = -».0 
XF INoL.gT.O) go to 3oG 

UN ( i) =U0 (IJ =HAlHfRCO,N ( L 0)’^r~UNC 

IF UEX NX .G r . 0 . 0 1 ) UN(n=i./ (l./UNd) 4-RESNX) 

G J r G G X 0 
330 GONTINUE 

UN( =00 (JhIP, )=UCIR0 

ON ( jWA) =Uo ( J wA > = U.N ( JWAf 1) =00 { JWA <■!) =UAIR •— • 

UN ( 1) = JO ( 1> = kCON ( 10 ) ’■T' UNC 
Jl=JWA+2 S J2=jHAfNGL S <X=1Q 

IF (NGL.EU.1) go to 400 

JO 090 J=J1,J2 

R<n = .<>1+1 S I1=I1C0N(J) S XZ = I2CCN(J) 

— aGiUnE S^ACG EETnEtN Gi_AZInGG TOO SmAi_u FOR CON VEC TI ON 

UN(J)=uO ( J) =RC0N (KaJ *-rFU.NG<-RG0N(4) 

300 CONTINUE - - 

400 CONTi-NUc 

<A = iCttlGL 5 J = J2+1 5 II=I1C0N(J) S I2=I2CON<J> 

- UN ( JJ =00 ( JJ =RGON ( KX ) ♦ T FUrJC<- 3 AiR 

IF dESNl .GT . 0.0 1) UN ( Jl) =1 ./REGiNl 

AiO continue _ . . 

— OP^aOUE INSULATION ON OUTSIDE OF SOUTH HASS HAui. (RSHWl) 

XF (kSMwI.lE.0.01) go TO ^20 . . .. . 

NG Oi'«iN = NGONN •*•1 S fGN 0 _/E = M NO 0£^- 1 
I IGvjN UWm) = IlUOiN ( 1) =NN0D£ 

J = N C .N’ N S 1 1 G G N ( J < = 1 ■ S I 2 CO 'tlO ) - NN 0 0E~S— ACON t J 7 -Iri 'J 

UN ( J) =U J ( J) =1 ./r.SHW I ' 

<F(NNGuE)=l S T(NNUO£1=70. _ __ _ 

ISOu(l) =NN0J£ 

>♦20 CONTINUE 



£n0 FROBwEH DATA 



-- PRINT DATA 

IF ( H » L • £Q « 1 • 0 K • K WA LL . EQ. 0 ) PRINT >+40} INCATE»UwOAO 

xP AN'^AwL.c. 0.2) PiX.lrN 't ^ 0 f INuATEfULOAO 
F.iiur ♦flu, CPHCIU 

IF (NiiEOPfA.EQ.OJ GOTO‘+3 0 - — 

PRINT >, 60 , (J , IICUN ( J ) , I2C0S ( J> fUO ( J) f UN ( J) , AC ON ( J ) , J= 1 , NC ONN ) 
PRINT 47 0 , (IAH5 ( J) ,J=1 ,NJA) 

- Print 430, ( J, I3UL ( J) , NO;. Z ( J^, TI*_T iO)-,WAZiK-J) ,-AL-FA(“J^, AGlZRJTn 

IOHAhG(J) ,USEPf^(J) ,TGuZ(J) ,0A3S(J) ,J=1,NJSJ 
PRxNT 4 90, iN H >•( A, iC, NTxHE, N SiG, N!10,NCAY1, I TH A X,Ui’T ITS, KOATAfK WALL, _ 
1RCauC,Kj£LT iKCONU ,NSuW,\AUA.*L,NJtir^T ,JuP ,JUfi ,JWm,x/\xR, 

2 J M X r; 

PRxnT 50 0 , QH DON J , Ai_ A T , G L 3 NG , Rx NO c. A , T HORN , T E VEN , T ZER C , QO N HA A , RriO , 
rCwT,TG-’'ixriN,TuHlNO, riLTC,MZIHu,RHO^,pAuI, iGilNXO,THCON,THICK 
2, i/J»SP, JLuAJ,R£jNI , RLtdTH ,RHOsR ,ASPR«r ,FC3K, CP HR, PSHWI 
IF ( N 3 G w , i. E • 0 ) GO to 4 3 0 ’ 

PRI..T 510, kWiNi , UWGUT, TriCH, TriKW, VSPW, AW 
430 CONTINUE 



H',0 

45 0 

460 

470 
43 0 

490 



I" OR HAT ( / >t 1 H 'W « TE P. WALL ■ 
16nJ_ut<D = F6. 3) 

FOR.'IhT (/42HTkJHS£ WAlL 
1dHJ>.0w0 = F6. 3) 



Ul'lX I O M U i 






lDnJb.UwJ = rfa. o ) 

FORHi-T (/2'4H conductance G'jNnE 
12HUo, 10X,2HUN, IGa, -fN>^>^0N/ (3x6, 
FjRiMmT (/26h mHSIEnT 7E.JPE\-<TU 
F0\NAT (19H SOi-A''. hEAT o0j:R0ES/4a,/ 
1, pHrSKZI*'i ,3 a, ‘■,HmLF>a, 3a , ♦HAnuZ , 2X, vriU 
2>tH0MGS / (3l9,2F7.2,F7 »h,F7 .2, >.F7 .>+) ) 
F J RH A T ( / VA , -rH NH A A , 6X , 2 H I C , 3 A , 5 Hi. T I 
1 pHi^UmY 1 , 3/( , 6HI T MmA ,3a , p:iUr(I T u ,3 a , 5 
■♦ / 9 i 3 " 7* " 



- PEi\ UNIT WAl... AkEA - IMDArE=Io ,5X , 

U.N S/>*A , 1 Hj , 3 A , 2 n I i, 3x , 2 HI 2 , 5X , 

. 4 ) ) 



''E connect . 

■ ‘ 6, 3E12 . 4) ) 

TURE nodes/ ■,X,2HI = ,1uI5) 
£S/ 4 a,/HJ _ISCw,_lX,;.hf, ' 

2 X , 'J d U H A N >j» 2 a, “HJScPr, 

T I He. ,hX , *+HN3cG At 3HNHG , 3X , 
5 rH K J A T A 



tui.Z, 2X > 4H I xlT ,2a 
,3X,“+HTGlZ, 5X, 



/ iAif^HNlriALL i3Af 5hKCAu3 f3X»pHK0bLl »3A»5HKN./iJNUi‘^Xf^ S G W > 2 X ^ 
dmK^U^CL t iX > pH K ^yEMT 
/7ib 

/pAf3riJGP#5X?3HJCW,5X» 3Hj7iA f4HJAXR 

/ p ^ 0 ) 



( 

I 



I 

I 






I 



I 




I 




i 




1 7a » r h 7 £N/ o El Z* -r / 5 a » 5 H TZ EkJ , / A , 0 mU JiN.i Aa QA t t ♦ ZHE a t 1 ]a , 

23HrMC»9AtA>rir0k-T/6£lZ«^/2< fontw l^,^.4t6A»6HT0iI’iJto/\ iJ^^IlLTCtZAt 
3:>n7^ZlMC»7Xf^H^MOt-»3AfAjrAoI/oEl :? a f o>‘^ \ C ^ ^ kJ ^ oX »3iTrt33;i ^7 a t 

/ A»:?HyGwSPt7At3HJLJA^, 7A,'^nrs£o^I/c£ IZ.^^/^XfoHALiNurHi & At 
55H^vnG:)rt7AtoHAjPr. m 7 toAt •♦"iFOiKtOAf-^nZP^x, 3 in,x:5‘^hI/jElZ* ->) 

510 F JR't^T ( A » 5HPWIN3 t 7X , 5riUW jUT ,7a , 4HTHCW, JA , ,i A ,^H*/SP;r,3X , " 

IZrlAn/oElZ**^} 

kE T URN 

END 



SJG/\CuriNE_ 



CONNCN/iLAD 1/NV, Nr,NT,NMAX ,:<Fl5U} , IC TZV (5 0 5 0 rri'T (50-^ 

1 KiL»JiNbtKiw» , ir^AA ,NE%(\, UNiT o, \ AJaCv. ,\.^AuD ,Kjc.wT ,KC0N(j 
C u^IrlOi^,/ 3L.ADZ / tiiAfictNTx'ic., ^ri.fl£,uMrE,UAlf,»jMYltuMfZ jNDmY (1.Z) j 
1 . i 0 , M J i , Y ro N N J , N L A Y 1 , I N| :. ; T £ 

CO i.ioN/iL A Ji / [(£ C ) , TO (5 j) , CP.i (5 0) , S (5 0) , 5P ( 5 j) f C0N0( 50 , 50) , 

1 o C u N ( 5 0 ) , 0 r A G u ( 5 0 ) , E \ T ( -^ 0 ) t r N ( p 0 , 0 , Z ) , F '^ ( 0 , 5 u , Z ) -- — 

COMiiO.N/dLMj^/A(:?0,50) t3{5G),Sj3\Gl(5Q), iOOPCZ (50) , SUC JNl ( 50 ) , 
i :>aoGNZ ( 5u ) , Qvmj ^ QVC:. , OAHT , OACu, UhZ , QiNC, UTRAN 

COMilUN/ iLAJ5/::.UNKT, UCIS* ,QCI riG , FAG , FRaC-, 

1 PI,riNEx,0£GNMu»,T5A:>r. tODtliAAtTiilN, T 0 A R , FA C I , R EM A I N 

GJMilv>'</3wMLio/i^ox*-/(pO) ,0CoN(?0) fOolZ (IjO) tOOolZl (loJ) , 

1 SOulZZ I 100 ) , TG’IIN, TDMAX , TO HlNN, T G i 1 1> J , TC MA aN , TCMmaD ,TCOOL 
C J li M u N / L H 0 7 / T P L K G A Y , 5 3 C 0 N , 7 A h 3 , Q C 1 M T 

• - GvJMMQm / 8LmD 6 / iOAY f iGPE\YR tXNu^Y ^ ROAi AiQHGOr<7 , ALAT tDu GN G j 

1 RjlNC^A , TOw i t ut.Gf COEoEv jEo t-blNuAT , uDSGA T j OgP 
Z ,Mr\G\I7tUlCUTt/<Grf‘JT,KOGDwtiG:)HlfMOEriZ 

CGiriG.\/iLA09 / Taj I 50 ) , VElj (50) ,UmO( 50) , RHu ( 5 0)', N JSvIS Cir( 

1 ri^T(15),WAZlM(lp),AG«Z(13), inGuZ ( 15 ) , N J A , I A Mi ( 1 G ) , 

Z H^rM(15)tOnAMG(15)tOSEP\(l3),r^LZ(15),OAiG(i5) ,ISO^X(Z‘i) „ 

COMMuN / SLmO 10/Qj NMh A, KJ IFF ( 5 ) , tiA , H 0 , T Z £r^O , 

1 NCuNNt lGOi< ( 1 00 ) ♦ I ZCGN ( 100 ) , U( iuO^j^ACON ( 10 0 N Uj ( 10 J ) , UN < 1 00 ) 

GUrirlUi^ / iLAOll/TlLio ,mZX»Tv ,RMUo,*i/UC> i F^^G ,F GUs^ ^'AEGNXt — 

1 JGr'tXND 

CGMMUN/ 7ROM jE/ VOlF , RCG^ (5 0) , K7E.NT t lAIR, JAXR , RLNG7H, RhOSR » 

* 1 A S H A H 7 t F C G R , J W t T R C 0 A 7 t J G‘ W t K W A G G 

CU M MON/ i. C Fu Gg/ is rc.C Pf \ , i 0 1 , I UZ , iPRl N 7 H , I PP>.I NT 3 , 1 PRI NT M t IPR 3 M , KHEOPR 
GDMMuN/:ilM£U/<^NAr\£A(iOG),iH<M,REA(lGu) _ _ __ 

C JNMun / Tl;iE«i/INiGH7 , THORN , T £ 7 £N , 3 £l 71 , DEL T , TIME ‘ 

QI HcN::)! ON iN VO ( 2) 

Q _ _ . „ 



1,'UiiGEk- UAY tDAYl--,-OAr2vr:-iiuA-rE-ryrM-T'J 

rs 

'' TC.'1iN = T:MIHN £ TG.^Aa=T3MAXS 

i1Ol-INC;TE/ia00G S iMi: W = I.MGA7E/iQ O-MOl* loa 
IF ( 10PE% YR. EQ. 3=6) NDAY(2)=29 

JEL f l = r? E k J A Y / FL 0 A r ( N T I Me ) ■ - - 

Xr* ( A6*EQ*Q) iTMAX—l 

r> 

"" ao IDO J=i,NCuNM 

Il=xlCCN(J) 3 i2=I2CON(J) 

CJNl. ( 11 , 12) =GONu ( 12 , ID = JN( J ) *AGaN ( J) 

10 3 GJ.'iilNUE 

00 233 i=l,.HMAX 

SGO.J(i) = 0. - - - — - - 

DO 2G0 J=l,fiMAX 

20 0 oGJr, (l) = SGC;DI) ^CONO(I, J) 

, 

IOAY =1M0A Y-1 

iF {.-lOi.wE.i) GO TO 500 S H2 = H01-1 
30 -rbO M = i,-I2 
^50 iJA < = j.GAY»-NUAi' (M) 

5u0 GO.SjINUE 

% 



SINwAT=S1N(OEGRAO^mLAT) 5 DuglAT=COS(DEGKAD'^A^AT) 



IF ( n.T: .LE 
GJS rC = CGS ( T 



0.) GO TO SOI 
lTC'^uEGRAu) :> 



Fci S G = • 5 ^ ( 1 • ^“03 3 r C ) S F3G3 = *5' 



SI nTG=S.N (DEGRADE 7 IlTC) 



o 0 0 C J N i I i''i J £ 

-- IF 7£MP uIMITS A-^E WITHIN 0*1 Jtu 
•• FUR Awu TIrlE aNJ KIC IS SET Tu Z 
iMNuU = A5S(TGHAAJ‘"TCHlN3) S iA f , 3 > ) 
‘ ' IF (oaNuj*GT •0« l•OR•dA^^ON*‘o^ .O. i) 
KxG-Z S KF(£C)=^ S T ( 1 3 ) = T C Mi N)N 
605 C3N r INUE 



( 1. -COSTS) 

T (IS) 



, i \ j. o j IS 

M S M F u A b 
mSS (T C lAXiN-T 
GO 1 0 o 0 5 — 



FIXED AT 
rlINN) 



TCMIN 



» COuGcGT AND ZDE^niFY MODES 

h7^!iF = NI = 0 
30 oGO I=l,NMaX 



/oi> 



;ooc50 



,< = <K(1) £ IF {< .;JE. 1 .,iND 

ir Gj rO D'+O 

GO fO (ol0,6^0),K 

C 7E:1? NJjEi 

£ ^ J i:< J) =1 
■IP MUGEo 
£ 

S 



► K*'I£#E) Go TO G'+O 



£ GO TO 3 OO 



610 N y = N V *• 1 

C PIxEu TE' 

62 0 iN r = N F 1 
06 3 N f = uT 1 
o 4 0 C 0 'I I i N J I . 

C- “O'vJ N T rvO L i»UGc. PUT AT ENO 

ur ( ) =uV {NV^-l ) =L T (NT 4 - 1 ) =iC 5 NT = l 4 r»-l 
IF UF { i : ) .EO. 1) (W = NVH £ IF (KF(i C) .cQ. 2) 



LF ( f,F) = I 
UT (NT) =i 



iiF = NF 1 



C 

0 

c 

G 

C 

c 

c 



-- KCAEC”IS~Z~FOF 'SI MP-E -JOUL ATION-Irr-FH: OH " ALT' CONC (iN J )-," CP'-t (I ) — 

G?(i)f ‘iNU OElT A-it QONSTw'ir FO.\ mLi. TI'IE nESUwTInG In CONSTANT 
-- m(I»J) Ar^kATS Fur NIGmT ANJ Ortf FOk CONTROL Tcik FI\£j AND 

— CjNTRjl Tc-NP VA^^AJlc (FOUR SETS uF C OEFF I C lEfiT S) — 

-- Cm^Cu^ATE COi(STAnT CCEFFiCiEN TS F(1,J> (FN(I,J,<) OR FO ( I , J , K) ) 

-- IN - T ( I ) =SUHl F( I , 0 ) *6 ( J > ) ; 1=1, NU 

-- <oe_t 15 zero for ino partitioning of tine inckenent when operating 

-* flOOc uF CONTkOi. node CHANGES ~ ROLlT MUST 3c ZERu FOR KCmLC OF 2» 

— KOJNU IS ONE FOR' Ur; ANO UL) CON'OTANr FOR ALL*" CONNECT I ONE 

-- KCOnU most 3E 1 FOR KCALS OF 2 

IF ( KCAL: .NE .2 ) GO TO 10 00 

K0EuT=0 £ <COMU=l 

00 -iOG KNC=1,2 - 

00 ZCQ J=1,NCCHN 

I1 = I1G0N(J) S I2=I2C0iN(J) 

Go TO (DpC,bS'p),KH0 
650 U(J)=UN(J) S GO TO 60O 
6j5 O(JJ=U0(J) 

6 dQ CJNoai, I2)=CONG (I2,11)=J(J)*«CUN(J) 

700 CJNilNUE 

00 720 IT = 1,NT - - 

11 = i-T(IT) S SCO.N(I1)=0. 

•Jj 720 JT = 1,NT 

12 = .T(JT) £ SCON ( Il)=SCON(rir+CONO( IlG I^r 

720 CONfl.NUE 

i.F (KF (IC) .EU.l) N'y0(lJ=N7 £ IF ( KF ( I G ) . £0 . 2) NyO(l)=NV+l _ 

NVO (2) =N0O( 1) -1 
00 dCO Ia=1,2 

Ni/Oa = N\/j (IX) - - - . . 

00 oOO <=l,NVOX 
DO 750 I V = 1 ,N VOX 

3(IV)=0. 

1 = uV(Iy/) £ rAK = FAG S IF ( CPM ( I ) .wE . 0 . ) FA< = i. 

0 0 / 5 C J y = 1 , N V 0 A . 

J=Ly(jy) £ A (ly, JV) =-FAR*COND (I, J) 

IF ( i. .EG. J) M ( 1 y, jy ) =CP1 ( I ) / jEi_TltFA<» ( SCON ( I) - S?(I) ) 

7 5 0 C 0 N T I N U E - . . 



3( K) =1 . 

Cmll. lEOTIF (A , 1, NOOa, NMAX ,3, 0 .WKArEA , IER) 

00 6CQ J = l, NVOX 

GO 10(772,77a) ,<N 0 
772 F^l(J,^, IX) = 3(J) S GO TO dOQ 
77-r FO ( J,rC, IX) =3 ( J) 

300 Continue 

IF (rvhcDPR.EU.O) GO TO 900 
print 3711,KN0 tlA.NyOx 

00 3701 1=1, NVOX 

GO TL(3oyi, 3 dU2) ,KNG 

3091 print 37 03, (FN (I , j, ;< ) , J= 1, NVOa) S GO TO 3731 

3092 PrIhT 3703 , ( F j (I , j , Ia) , J= 1, NVOX) 

3701 CONilNUE 

37 0 3 FOrMAK 10£12,A) 

3711 FORnAr(» RN G=*I5 , 5X*IX = *Ip, 5x*NyOA = *I 5//* I F ( I , J ) , J= 1 ,N y ♦ ) 

^37l-,FoRNAT(I3,l\3El2.4) 



900 CONTINUE 
10 OJ CjNi'InUE 
Nc r URN 

ENO 

su-3;^cuf In£~pr5p 






C0MNUM7 dL>.Gl/Ny , NF , N'f , N mA A , <F (50 ) , iC ,L'y ( 50) , lF (5t3 ) ,•_ T (50) , 
1 KIurli«o,Kiu ,Il*HMA,Nc"'x>*\,Ui*'i* To,RA»JAGc,Nv'-4i-.o,f\Jc.. T ,KvOIU 
Go '1 rtON / 3 L « J 2 / t«T i Me , M r I H £ , j. T I Nc , 0 A r E 7 •-> A / , 0 « T 1 , J •< T 2 , r( 0 A V (12) 



1 MO ,MG1 , Yk, NMO , N uAri , iiOA TE 



G‘0 M M 0 N / 3 L AO 3 / T ( 5 0 ) , T 0 ( 3 0 J , 0 f'M I 5 0 ) , S ( 5 G I , S P _( d 3 ) , G U N 0 ( 5 0 , 5 0 ) ? 

1 o C o N ( 5 0 ) , U F A C 0 ( “J 0 ) , c R R i (GO) , F N ( :? C ,50,2) , F fj ( j 3 , 5 0 , 2 ) 

El ^ .U . MC -I 1 ■ Vi I 1 I T li I .■CJRR.IP Ml . N.l L n.N’1 I q. 



I 



u 




; c 



l SQouMZ ( ‘it) ) ,QVrir , J r . , JT.^Af) 

^ ^ ^ / A ^ - 'i ' i 'N r » . ^ i ^ I < ^ ^ • 

1^ If I iH“AfUcia.T.A4^, [ j*-i': t i> ,r^^Af F'ij.N ?Ti»^*\fFAGlt.'^cMAiN 

LiJiiloi / tUul«£ (l3lj)f5'-«GLcl.(l0j) y 

I ^ ^ '•J ^ 1 f- G ( 1 w j ) f T 1 ^ i » I j* j t \ f r ^ 1 iM' I y r L/ <i I r ^ ‘3 f r G M ^ X » ^ 1 1 c M M A J f r 0 0 u L 

o>J»Iiioi / 3L»-»G/^ / r P „ f j ^ X ■* t /-»j’lz3f'-4 0x TIT 

OG-'JjIo % / dL^iJd/^3-ffIJ^r. '^t^fISJM/fr\ij^T#%f 3HL»0**'^VyA^KT yO^Ci^Gt 
1 r^lu'^LAy r fwi-v t >-»J-:)->t.^ f-iI.N^Cv>f -)xi<^AT tCJu-A f 

f jxojr 

CJfinu^/ r ( :,o ) , J ( pO) f xJnJ (‘3 0) y (3 0 ) t N JS , IS OL ( 1 5 ) » 

I rXl«i ll?) 13 ) ft'lJv»Z F Id) ) f 

^ ALFAUd),G.-iA:.^{1p) yC3EP\(iD),To^Z(15),3A3S(lD) ,Ib0wX(24) 

G J M li u N / D k. M j 1 u / 3 3 ^ A X • .*\ S I F F ( 3 ) i £ X t ^ H ») » T Z £ .a 0 t 

1 NJuNN , : lCo:i ( 1 JU ) • I 2GUN U 03 ) , U ( IJ J ) , ACJ‘U LOG )-y UG( 100 3 tioiy) 

CJMHoN/ 3lAj 1.1/ riLlC fAZI-iC yXiGUtSuorCfSIor Of FioCfFwOO y f\£ S N I y 

1 JuF y/Sj^'iO 

C0Hr*CN/ TPsCflBE/ yivLFyr.OCN (50) y K;£NT tlAI.x, JA x ^ , X^N GT h , ^ hOSR , 
t ASPr^Kl ff-LSrsyJ^Uy r tJuHtKrtAuw 

oOM. ION/ I 0F(.GS / iirvCr*.-.,, 13 1, i02 , i?<iNr H , IP ^ INTO , IPRInTH, IP^S rt, <H£0PR 

0 j M Huii/ r irtti / r .'K OH T » r iio r.i< » r t ven » uel t i » jel t » ti l^£ 

INTE&Eri~DAY ,UAri , C A Y2 , T R, 03 T UNl TS 



— <V£NT IS 0 FCR NO VEf,rS OR VENTS ALWAYS CLOSED 

i<V£Ni iS 1 FOr\ V£^^^S ALWAYS OPEN 

''CVtNf iS 2 FOr^ NO •\EYErt.SE VENT FwOW (R'/F) 

-- RVENT iS 4 FOk ThE.xHujTAT^CA.LY 00nTRCl..£l) VENT (NO RVF) - 

-- VOuF IS VOLUMETRIC FLOW RAT£/GuASS AREA 

IF ( KLuNU .EU . 1 ) OU TO 75 ’ ' 

HAIR li OUTSIDE Al.~, FIlM COErFICIEiNT 
HHli\=2»G+-4»0'''VE^0(ITIME)/o«70bo 
0 ADO J.RANiPMfRANT OUFSxDE RESISTmNOE FOR NGL^Q 
IF(^.G_Z(l).EQ.O)h^IR=l./(l./MMIK^KGON(5)) 







VOLF=0. 5 IF (<VENT .£0,0 » GO TO 40 
0T0f=(T(iAir'J-T(IG))/(r(IAIR) +TZERO ) 

VJLr x = RCON{ D-^SOkT ( ASS (Or or ) r S IF ( CTOrvuT .■ i vi—'JOi.f X--'/ QlF^- 

GJ 10(10,20,10,28) ,<V£NT 
10 V0lF=V0lFa S go to 40 

20 IF(VOlFa.ijT.O.) VOlF=VOLFa S GO TO 40 
28 CONTINUE 

IFC.MTINE.EQ.I.OR.KIG.nE.-I.OR.VO^FX.lE.O. ) GO TO 40 - 

A T Rm = 0 . 

QJ 3C ir=l,NT 

---I = L1(1T) S IF(I, £Q.iAiRl“GO--T0-3C 

XT.Rm = XT\A<-CONO(IC ,I)»(T(IC) -T(D) 

30 CONTINUE 

Ar.RH = xTRM-3 ( 10) -SP (IC) T ( 10) ■ ■" 

VOlF=sTxA/AOUN(JAIR) / (T(I/,IR)-T(IG) )/RCCM (2) 

Ir(v0wr.LT»0») V0i-F = G. S IF(VOLF.GT»VOuFa) VOv.F = VOLFX 
40 OONTINUE 

-■■iOA.CulATE SPECIF! C-CONDJOT a NOE S--(UNv-UO)~FOR-TE.MPErAT-JR£-ANO — 
- Tlit jEP£iNjc.NT ^ Ci(N EC T IONS 

NGL-NGLZ(I) 6i"(M\jL»Gr,0) jOTJ4p 

“ NO CuUwEOTlk glazings 

Il = -^lC0N(i) S I2=I2C0M(1) 

UN ( 1) = U J ( 1) = H» xK + ROON ( 1 0) -‘RAuFN ( T ( 1 1) ,T ( 12) , TZERO ) 
Ir(NShUT.Ea.i)'JN(l)=UG(l)=l./ ( 1 . /Hm IR4- \ESNI ) 



N G G T 0 Z 3 

i — nG- gOwLEltor Glazings with .•vAoiKTiON/coN vecTx'cn gon-iectioms 

^5 G 0 N T I N U E 

VFluf^-ASS ( VCLF) 

UN ( JAIR) = UO ( JAl\) =RCON ( 2) V- uGW S UAIR=1. 4-. 21*RCON (3 ) ^ VFLCW 
UN (JWA) = UQ(JWA)=UN(JWAH) = U0 (JWA*-!) =UAIR 
I2=:2COn ( 1) 

UN i 1) = 'j j ( 1) =RGON < 10 ) *-RAD- N(T (1) , T { 1 2) , T ZExO ) 

Ji = ji'A + 2 £ J2 = Jii A 4,‘^ Gl £ KX-iO 
IF(„Gu.EQ.l) GO TO oO 
UO iO J=J1,J2 
<A = r,A4l 

Il=ilCON(J) 5 I2=I2C0N(J) 

UEv. rs = A-3S(T (ID -T (12) ) 5 TS A RS= . i M T ( 1 1 ) + T ( 1 2 ) ) H ZE R 0_ 

' HCrxT is T ml MlillMU'l (.TTAIFiASLE OOuVLCTION COEFF. IN '/LRIxCAL 

AIR SPACE - - - - - 

HORxT = l.G32 4( QEw TS/T6ARS)^“0.333 
H 0 U N V = R 0 0 N ( ) S I F ( h C 0 N V , L T » H G jR I T ) H 0 OM V = H C i< I T 
UN ( l) = Uj ( J) =RCON ( KX ) »'f,;ADFN( r ( X D , T ( 12) ,T ZER C) 4HC0N V 
50 GJNiINJE 
60 CON riNUE 

L <X = 104N3L £ J = J2H S Il = I10ON(O) £ I2 = I2C0N(J) 



M 



i 



I 

I 



( 



m 







c 



75 



Uj{ J)=Jj(j) = KCJ‘.<<A)»^4Grfi<r (ID 
IF ( NbnU T . f,£ . D .jJ TO 7b 
‘JJ ( J1 ) = Jli ( J U = 1 . /r^tiSN I 
G 0 iN 1 1 1< U c. 



I ( IG) ,rZ£R0) <-HAlF 



O 

0 -- oAwGOLk-rE U( J ) ,GJfiJ ( 4. 1 . 1 2) ,3G0N ( I ) 

C -- this ?h\T iHGJLJ at <tpr Fu.-i M.(r r-luOELi 
x=Ih^3( D i TO rTAHS 
OJ iOC J — ltNGut' 4'1 
il=xluON(J) i I2=I2CGN(J) 

GG Tu (9l*'33)»<iu 
91 U(J)-UN<JJ 0 GJ TO 95 

— ■ 9i •J(JJ=^C(J) 

95 CO.Judli 12) = C0iJu (12,11) = J ( J) *ACON( J) 
ICO cj^in.'Ujt 

ry 

IF ( <IGrt,NG.NE .2) GO TO 200 
•JO 120 J = 1,NJA 

I = ImK 5(J) 5 IF (,<F (I) .EQ.2) r(I)=TA.93 

120 CONTINUE 

* 200 GUNTINUE'" ■ 

C 

00 oOO IT=1,NT 



Il=wT<IT) S SCON(I1)=0. 

0 3 oOG JT = l,iNT 

I2 = l.T(JT) 3 S‘..ON (Il)=SCON (ID vCUNO( Ii,I2) 
80 0 GUNTI.NUE 
K.C. T uR.N 








EN3 ■ - - - .. 






F J r< 0 T 1 J N KA u F i>» ( T 1 1 T 2 1 T 0 ) 


- 


_ _ ... 


TTl=Tl>ro £ TT2=T2VTQ 

PA0FN=( rTl*-^2+TT2* + 2) *(TT 1 + TT2) - 

PE T u F N 
£N 0 


--- 


■ 


5J3kuJT:nE JAYLY 



Go M i1 U N/5 lA 01/N 9, NF, NT, S'iAAt'CF (50) ,IG,(_V (5 0) ,or (50) ,lT(50) , 

1 KaoHNo ,Kxu , li NhA ,Nc.‘^R, UNI I S,<\AUaLw ,K-^AuG ,K0c..«T ,KC0NU 
GJM.-1GN/3LAG2/NTI Nt , MTI9E, ITIHE, OhTE-, OA/, OA r 1 , GaY 2 ,N oAT <1 2) , 

1 MUjM^l, YF,NiiL,Nt.At'l,l'('0Afc. 

CoNMON/OLAJi/ r (5 0 ) , TO (5 0 ) ,CPci (5'0) , G (50 ) , S? ( 5.0 ) , CO N0( 50 , 50) ♦ 

■- 1 SCCN (50 ) ,QFagO ( 50 ) T ( 50) , F N ( 5 0 , 5 0 , 2 ) , FO (aO ,' 50 ,2) 

G JiOHUS/ OLMG-t/^. (;. 0 , p 0) ,3(5 0) , oG S■'^G 1 ( 5 G ) , GQGPG 2 ( 5 0 ) , GuCONl (50) , 

1 S'ToGN2(pG),Cv/rii,U'/CL,3AriT,GMCu,uriZ_,QINC,QTFAN 
C*j9/iaN / 3t»AGp/aU9,NT ,GGA-(,Guii4vJ,FHO,^t\AO, 

1 P I fTiMaAiOtGf'AG , T aAat, JJ , T■iM/^, Ti^^ilN ,T3 Aia, FAGI ,PtMAlN 
C0.9riU.-4 / 3 1 . AO o/v* it L ( 5 C ) , G C G N ( 5 0 ) , GG 1 2 ( 1 0 0 ) , oJC 12 1 (1 0 0 ) , - - 

1 ujwl22 ( 100 ) f tu.-'ixii, T’w-iarn'i, 1 Ci*ilNG, TU'**lAXii,Tut--lAXj ,TuOOU 

GG .•Iriu.N / 6 L AO 7 / T P£ KG A Y , S 3 COh , T h 9 3 , (i uI NT 

CO .iriCN / :3 !_aG 3 / i JA Y , I OPE^ Y < , I NOa Y , .^U aT A ,■ OHC-GN 9 . AL'AT ■, Ou CNG , 

1 r\ii(CaA, T&uT, G E C, 0 0 SOEC, 51 ii 0 EG, alNu. A T,C OS L A T,GG? 

2 ,Hi\C'4.xr juIGUi ,<.ariUT,\GtjGL.,.9Uar11,f’toSH2 

CU,9.-1G.\/ 3LAa-j/ ThJ( 50),7£-0(50), Gi-U (53), PHO (50) , N J3 , 13 Ca ( 1 5 ) , 

1 riuT(15) ,XaZ1M(1p) ,40.2(15) ,n3uZ(15) ,iiJA,iA>13(ia), 

2 mLFm( 15) ,Ofl.-''44j(lp) ,0oaP.'\(15) owZ (15) ,GA3a (.Ip) ,iSO'.A(2A) 

CO .lilUiS / 3i_A J 1 Q / GO N 9 m.( , r'vD I FF ( p ) , EA , r;rl G , T Z G » 

1 nGjn.N, 1 ICON (IOC ), I 2GGN ( 10 0 ) ,U ( lu 0) , AUON ( 1 0 0 ) , UDj 1 0 3 ) , UN ( 1 0 0 ) 

■ ■ GJ.iriuri / 3LA Jl 1/ Ti u T G , A ZI NC , kNU'J , GO uTG ,S i.NTC , F'GSG sue ,.<ES NI , 

1 JCr,<.N3 

CuMnCN/ T aC'NSE/ 7 GlF , ,\C0N (5 0 ) , K vE.iT , I A i.^, J A IP , PuNGT H,.PIiOSP , 

1 ASFr^.AT,Fu3^^,J^3M,T•xCGA^ ,JGH,KWA.i_ 

GGMrlCN/iGFLoS/IS.TCF.'\, iJ 1, IG2 » iPAiNTH ,I PPs IiNr0,I?.PIfiTM,iPPSf1,K.3£0?P 
oOiTtivT,/ ^xNEa/iMIorlT,^^l0^x'^,Tw7E^J,G£^-Tx>0EuT,TANE 
JI.NcNaiGN U,‘iGa(2a) 



' IN TEG£:^ UaY , OA Y1 , Oa Y2 , Y P, uA T£-,UnI TS 

IF (PLATA. NE.O) GOTO 1000 

3 PEAti DATA F OEx EACri 'jF 24 HOURS F.PC.N TAPE! 

% 

4 TaJ la 4M3IENT TEMPERATURE (DEG.C3 

3 7cLj is WiNJ 7E.0CITY (M/S)_ 

3 - 'QhO Iu SO-w.-4p PAG I A T ION ■ 1 N f tiN 5 I T Y"< W/SQ .9 ) 

IF (MTINE.GT.O) GOTO 130 

3 fM<t uArE UF FINJI.'TG ‘7AlJ£ 3 FO.R FI.pST TIME THROUGH 

lY P r = NC J (1MGA3 E , 100 ) 
j.NOhTT = I YrJ'* lu 03 0 10 0 v IN DAY 

50 PEAa(i,3G50) Y KMu , K UA Y , ,<H J UH , QHO ( 1 ) , 7ELu ( 1 ) , T AO ( 1) 



oO 



IF IITEMP-I NCATT ) 5>Q,ilQ,60 

P<i(.T ♦ , xST^^riNG GATE TOG EAk^Y pj^ FILE A , I.SJA TE 
G 1 0 A 

j XcAl; GATm. gate '-IUST .IATGh, O.A. ERRO-^ 

lOii Re Ac. ( i t iCbO ) r.R , K ,''0 f KG A r , RriJUr\ f Lnj ( 1 ), VEwG ( 1 ) t TAD ( 1) 

IF ( ( KM G . EU .MOr. AN J . ( K j Ar . E J . G AY ) . ANG . ( <M JUR . EG . 0 ) ) GOTO — riO — - 

PRiNl *>?GAY n»o/l.<7GH IN 'rt£ATht^c F i i.E^ • f R i K>IG » K GA Y , KhO OR 
SIGP 

c rEAu in 23 VA^.JES FOR IrtlS GAY 

110 GG 130 » = 2,2L 

130 R£Au<1,30d 0) Yr^, K.MO,<QAr, sHGGa,aiiG(u) , VEi_ 0(L) ,TAOa) 

OArt: = \MO-*'10Q00>KDAY^l0 0 4-YR 

c 

c 00 CONVEKSIOM Or-WEATHER-OArA— IF-GTSTNO— ASHRMr-TOnrrO 

IF (UNirS .EG. 2) GOTO 220 
GO ipQ Ih = l,2.. 

C GGNvE'.T TEMPERATURE TO OEG-- 

TA0lH)=TA0(Ih)*l.d<-3 2. 

li?0 CONiiNUE - - - -- - . - - - - 

220 CGNiiNUE 
C 

C COMPUTE NAXTMIMrAVG-TEMPERrruRE-ANa-OEGREE-OAYS 

TMAA=-l.t6 3 TMIN=1.E6 o T3Ak=0. 

GG2AGIh — 1>2 “t 

IF ( TAJ ( Ih) .oT ,TMAX) THAX=TAG(IM) 
xF(rA^(iH).;_T.rMlN) rMIN=TAG(iri) 

T d A R= T 5 A R +■ T A D ( I H ) - — — 

24 0 C o N T i iN G e 

GJEG-T3AbE-.5'^(rMAX UMIN) S DG=GC +AM AX 1 ( GOEG , G . ) 

- Ta AR=TauR/24 . 

C COMPUTE GEClINATIOM (OEC FG<MULA FROM G. AN C 3.) . ^ 

GcC-23.-f5>**‘bIM<D-GRAO*36']./3b5.’*“(Fi-uAT(IGAY) •) } 

SIiNuEJ = SIN(JEGRau^JEC) 5 C05G£G-G03(GEGRAQ^0£C) 

H<G\ii=5G. S IF<CtiC»LE#0*) 73 2d3 - — — 

H Ps C r\ i. T = ~ 9 0 • 

CH\wR=T AN (DEGRAG^OEC) /7 AN (OEG^AJ^^AlAT) S IF ( CHRGR . GT . 1 . ) GO TO 230 

1 = A C CS ( G‘MRUR ) / CEGRAO’ 

26 3 CON j iNUc. 

C — SOuMA CONSTANT _ 

L7<J iG(olut320) tJMXT^ 

310 uSP = 426*9o-i3* 30'^SIN(CEGRA3^3o0.^ <272*l^FLOAT (lOAY)) /365«) 

GJTC330 _ _ 

32 0 QcjP-13Uo. l-^2.5b^3IN(CEGRA3^3Da.* (272 .H^FlO AT ( I OA Y ) ) /365 • ) 

330 CON TIN Uc: 

RE Turn 

c 

3 3JM.1Y DATA FOR CONSTANT 5UMMY-G4Y CASE 

G 

1000 continue 

OmTm l4'lHAfTAiMl)NthRTMAX?H?\T»Mi.N/70*tDfl*tl3#t7* / — — - 

S L P i = ( T ^ li A A “ T M M i [ 1 ) / ( H r; T M A A - H K T M 1 1 'i ) 

S P 2 = ( T A H I N ~ T M M A X ) / ( H r>. T M 1 N 2 4 ♦ ^ rt.x T N A X ) 

3 ■' -QUriNT SOLA^; DATA FQk SJriNY 3 A Y ' CA wCUlATI 0 NS ."“V^/ SQ 

j4Tm 3H3a/ 5^0*,3. 152jl57.1t->63.I,/bi.lflJQ::>*if 
^ 12'j2»2>13l:^*otl'3i>7 •6>1319*6jl2 0 2«2jl3tJ!?»pj7!;)O*lj4Q3»lj 
^ io / . I ) 3 • 152 , s + 0 . / 

IF (MTIME •EU* 0) 3mTE=IN3 ATE S I F ( NT I ME . GT . 0 ) 3.A TE= OAT E+ 10 0 

33 iC2Q L-iiZn^ - - - -- 

H = w ~ 1 

TAMd = TAMIh+ SlPI^ ( HR-HRT MI N) 

IF (Mh .lT •HRfMi.N) TAM3=T AM AA^S;-P2^■(^r.R■^24v-^iRTMAX3 — 

IF (nR.oT .nRIhAA) TAN 3 = 7 AMAa<"SuP 2* (ri R- H .PT MAX) 

TAG(u)=TANU S VElQ(u) =2 .2352 :> Ori u ( L ) =QHOx I L ) 

_^1323 G'o'i^ilNUci 

GO TO 220 

SOpO hOXNAT (7AfI2tlA> I2t 1X>I 2t lX>I2tlX>FH*Q>4Xjr 2*0? 1X»F5*1) 



u J j .*% 0 u T I E :> U N S < C 



CjM^iUN / 3L.MGl/N\/fNFfNT ^SMAx> ^F (3 j) i 10 9 L.'Ji o j ) Ct?0) »wT (i?0) 9 

1 XlunNSf KIG ^ITM^X 9 NEP.KI Uf^I fj 9 XAJal^u fXUAwG tKuu...T ♦RCUNU 
C3N!'lui'i/3LMD2/NTlMEt MT-.lcf xTII'lEju^fE tOAY tL)4/ ifJAf2>Ni3AY (12) 9 
lMJfi‘4 01fTT\ii4M3«»'iL»MYltI^13HTc 

C 0 M N 0 N / 3 L M 0 3 / T i 5 0 ) t T 3 ( 5 3 ) ♦ C ^ M ( 5 0 ) > 5 ( 5 *3 ) » S P ( 5 3 ) > C u N O ( 5 u t 5 0 ) 9 

1 :;t C l.4 ( 5 0 ) t OF A C3 ( 5 J ) t E.RP T ( 5 0 ) , FN ( 5 0 » 5 0 > 2 ) , F 3 ( 3 ^ , 2 ) 

G U M ri 3 iN / J L M 0 / A ( :? 0 t 5 0 ) ? 3 ( 5 C ) > S U S \ G 1 I 5 Q ^ t S Q 3 R 3 2 ( :? 0 ) > S 0 CO N 1 ( 5 0 ) t 

i 5 JuUN 2(50) iQi/MT f3 7CLfQ^«'^T? C A ^ l » 0 Z ♦ U I N v > l-I T FA r-1 
C J MMUN/ j LA 35 / S UM K. T > uG 1 N j 3 3 i NO f FA G f F PAG , 

1 PI >I xMci#\#GEG.NAjf r SmSE^ ^U 9 T "1m/\> I tMiN 9 TSaR-j F MCi?Rc.MAo.N 
C3 MM-wN/ 3 lAOd/OSRC (5 0) t QG O^N ( 5 0 ) , QG 1 2 ( 1 0 0 3 , SO C 12 1 ( 1 0 0 ) , 



I 



( 100) , TGnI D, , r JfliwrJ, rC.-llNJ , Tu.iAXN, TCMmaJ .tcool 

JuAj//rp'KjATioOco'i»rM.is»'4Gi'‘)r 
CoMiiO.'^/ Olhu^ / iOkY f f.NvjMjAiUHCONVfAuAT lUuCNG? 

c » H .\w \i.r iuICuT»<.w."ijr»\Gujw* lUorli-tfiJ jH2 

cum.tu'iv OLA j -j/ r A ) ( 'j : ) , j I j 0 ) , umj c j O) , -^na 1 5 o> ,-.n js , ig cl <i5) ■; 

1 riwl (l3) IG) ,AL.2(1 p) tMuc2(l5) ,:|iA ,lAHoliCJ , 

2 Ai_r».(l^) jOHAiSuCly) *'jS£'^\(ly) «rjwZ(l5) »JAij(i:j) »I S'Jw X ( 24 J 

L L H L *’• / J L M >J 1 0 / G ^ N M A A f J * r f" ( a ) y c A i t*\ r) L f I L C. 0 t 

1 iN Ju\N, I louM ( 1 0 0 ) , I 22o'-U 1 j 0 ) , U ( 1 00 ) , M'J JN ( 1 0 0 ) , U3 ( 100 ) , UN ( t 0 0 ) 

1 JGr- , nN J 

<^'3) ,.<V£Nr, lAIR, JAIR,RlMGTH,^HOSR, 

1 A 5 P R A T f F C5R j J W A , T RC'J A r , J 2 W , K WA L'L' 

GjNr1uS/IOFLui>/I3rAGPK»I31tiij2»iP\lNT H»IPRlNrQjIPrAlYJ.i»IPR3M|<h£OPR 
cj'),ijN/rifids/T.NibMr,rHo.-AN,rt;v£.N,u£wri,aEuTi Ti.1t: . _ 

Jl.ic.NSI Jfi £07(12) 

c 

I.NTdGER LMY,uAYl,0AY2,rR, DATE, 0, NITS -- - ■ - 

c 

Gmi m £0T/‘”ll*2t'*13*9»“7 •5tl«lf3*3t'"i«.+ f*"D*2t'“2*4?7*5 $ 
llt?*4tl3*2ijl*o/ —————— — 

C 

a.i = uHG CTIic) ♦QhCONV £ T A M3 = T A 0 ( I T I i£) 

lA-lAiJ(l) S T(IA)=T4i3 
G-'SJ-A.'cnEaTSOURGc. 

lF(NJi._£.0) aETU.RH - - - - - - 

OJ 140 I=itNHAX 
ItO J(I)=cP(i)=0. 

aH2=Oi.'i: = C. 5-QrRAN = (jT 

I F ( uh . L . p . ) KETURN 

OiZ-uH _ _ _ _ 

C -- SJwA.-^ rii£ AHO HJUk ang.e 

SJ.NT IMt = Tj..i£X+F^OAT (ISO-X ( ITIiE) ) *-EOT (.iC) / 6 C.-3 lGNG/ 15 . - .5 

riR= 15. » ( 12. -SUNT liE ) - 

G J S n.R- C 3 S ( 0 £ GRA j * H R J S i I I'iH ■«.= 31 N ( jEGkA J*H R) 

G bUN MuTiTUGc 

S I N « L T = 2 C S L A T *■ GO S 0 E C * CO SH R+ S I NL A T »b' I NO ECT 

IFtjIi'iAwT.LL.Ci.) RETURN 

A-T = A£In (SxNAlT) 3 COSAlT=COS ( A..T ) . . . 

G SUN AZ.k.iJ7.A 

SxNMZl = :ObO£C*SINH.R/COSALT 5 A ZI .i = AS IN ( S I NA ZI ) 

IF ( AtS ( rir%) . UT .H.T.CRIT) A Z i M = ( Px -Aa S ( A 2 I M ) ) ^ A ZI i/ Ab S ( A ZI .i) 

G -- GU-LcCTCr INPUT 

IF (TILTS. L£. 1.0) GO TO 153 

C A _ L C 0 ^ L E C 7 ( G H t A Z I i , S r iN A L r Tt; US Arrv PP )“ '5~GO " 7 a ~l 7 0 : 

160 OEXr = u5P-^SI NAL7 6 PP = Q,A/Q£XT 

170 FOJh = 1.79*PP-0.55 

IF ( fUL.S . l7. G. ) FUON = 0. i I F ( FOQN . GT . 1 . ) FUO.N=l. 

UJ.N = FG JN^UONiAX 

UO £ r\= uO N ^ 3 1 NAL 7 3 SuiF = OH-0 jER - 

G 

c — sjla.r heat sourcc. nodes 

■ 00 oCQ J = l, NJS 

I=IiCx ( J ) 

WAZ.l = w A Z li i J ) S NG^ = .NG_Z(J) 

3lHriuT = 3lN (uEGRA J^TIlT ( J ) ) S GOo 7ILT = C0S (DEG.RAO^TI-T ( J) ) 

J — — AHta.L'^xOL^K AZxMUFH 

GAii“.*Zx''1 — oc.uR/.D^WaZI i COc>G.^i“bOb(GAt1) — — ■ - - - — - — - 

s — amwx-3olar. Angle of ingidpnse 

COSxNS^SCSALT-SIiNTI^T^SGiGAi ^3INALT *003 TILT 

•" • iF (uGSINC.lE.G.) GUSINS = , 0001 

A..Nu~.nLOb(oUSx>).^) ^ C)xNli'‘iL“^j.N(AxNL) 

3 — F.xESkEw RElaTIUin for R£-LECfION 

TiXAiJi-l. 3 AtnEF-C. 3 xF(NuL.EU.O) uO TO 230 
NGF = 2*,N3L-1 

IF ( hINS . EG. 0 , ) GO TO 200 

A R Er =M 3 X N ( S In I .NS / RI NO Ex ) S X = Mx.NG-ARtF S Y = AINS+AP£F 
R1 =(S:n(a)/ 3IN(Y ) ) *--^2 S R2= (T/..N( X) / TA.Ni Y) ) ♦’‘a 
~ ■■ Tl=il.-<1)/(1. »<1 -^,NgF) S T2=ll.-R2)/(l.«-R2"^NbF ) 

T.rA,y3 = 0 . E-' ( T 1+T2) £ GO 1 0 230 

20 0 %= ( (RxNOEX-1.) /(RI.NDEXH. > ) *^^2 S TRA.NS= ( 1 .-.R) / ( 1 . NGF) 

230 GONi x^.Gcl 

« — ASbOroTICN IN GuAZING 

GxAo3 = ExP(-EX’-TjlZ(J)/S 03 (A.RtF) ) 
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PASOLE: A GENERAL SIMULATION PROGRAM 
FOR PASSIVE SOL.VK ENERGY 

by 

Robert D. McFarland 



ABSTRACT 

The PASOLE computer program was developed to do simulations of 
passive solar heated buildings. Modeling is done using a general thermal 
network method that allows for heat sources and thermal storage. Sun posi- 
tion equations are used with a global-to-direct solar radiation correlation to 
develop solar heat sources from measured insolation data. Models of a par- 
ticular class of south-mass-wall passive buildings have been developed and 
are described in this report. 



INTRODUCTION 

PASOLE (PAsslve SOLar Energy) is a computer program that incorporates a general, thermal 
network solution of an interconnected system of uniform temperature regions with time-varying 
boundary conditions and heat sources. PASOLE was designed specifically for detailed analyses 
of passive solar heated structures and includes algorithms for calculating solar sources in a 
general way. Because of its generality, PASOLE can be used to simulate hybrid or active solar 
heating systems. The PASOLE format allows the user to describe the thermal network model by 
specifying nodes that represent finite regions, connections between the.se nodes, and parameters 
associated with the nodes and connections. These specifications are made in FORTRAN sub- 
routines in the program. P.ASOLE is not a standard "user oriented" program; that is, to use the 
program for other than the specific models for which it is programmed, the user must have some 
knowledge of heat transfer principles and FORTRAN programming. 

In describing the program, we will refer to a particular class of south-mass-wall (.SMW) 
passive solar heating sy.stems — the network models for which thi> version of P.-XSIRE is 
programmed. .Although ikASOLE was not written specifically for these particular models, they 
have been used for extensive parametric analyses and serve as an example of how simulation 
models can be created. For these models, hour-by-hour simulations for one full year, u.sing real 
weather data, take about 1 min of central processor (CP) time on the CDC 7000. 

Most of the variables in the program are given in .A[)penclixes A-C. The FORTRAN listing in 
Appendix D includes numerous comment.s. A sample output is given in Appendix E. 



♦I 



0i 



I 



Tlll-UMAL NCTWOHK MODKL OF IMIYSICAI. SYSTFM 



I'he sv>(em‘s })hysical and environmental structure is represented by a network of N nodes, 
eacli representin'' a region of uniform temperature. Fach node mas be connected to :\ny other 
node bv a thermal conductance K, suclt that the rate of heat flow from node i to ncKle j is linear. 



where j denotes each of the other nodes in the network system. Many of the Ktj’s u>ually are zero, 
but Ktj is always zero for j = i. Many nodes have no heat source term (Si = Rt = 0) and many 
have essentially no heat capacity (M, = 0). When the node is mas.sless its temperature is not 
directly dependent on its temperature history, so it requires a sliyhily cliffereni mathematical 
treatment, as shown below. 

By rearranging^ Fep (I) to obtain the time rate of node temperature change (dT/dt)i, and by 
scparal ing variables and integrating over the time interval from to t, we obtain the node tem- 
perature change over time increment At by 



The average values of the integrands are a.ssumed to be a linear function of their vabie.s at 
tintes t® and t. Letting I be a general integrand 



where l = fl + (l— QI®- Although f may be varied in the })rogram and may take any value from 
zero to unity, a value of 0.5 usually is used because it is the most accurate and results in stable 
solutions. 

Some nodes, such ns outside air toin|jerature, have ’’fixed" or known temperatures at a given 
lime, whereas other nodes have "variable" or unknown temperatures that must be determined. 
Let XV be the number of variahle-tcmpernture nodes, NF the number of fixed-terni)erature 
ntides, and N the tot.al nunihrr of nodc'^. fly inserting the integrand approximation into Fri. (21 
and collecting terms that contain the unknown T,, the variable and the constant terms for 
each of the XV vjiriahle-temperature nodes, wc obtain a set of XV linear equations with NV 
unknowns. 



(hi = Ivij (Tj — I'j). Fach node may have a heat capacitance M and a heat source rate K. 'The 
source rate F may be a linear function of the node temperature K = S + B-T. 



I'hc heat balance at any given node i at time t is 




+ H. (dT/dt). = S. + B.T. , 



( 1 ) 




dt , 



( 2 ) 



where the superscript zero deiK>tes conditions at the "old" time, t” = t - .it. 



t 




dt = I At , 



1 



1. NV, 



(3) 
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j=l 



where 



= -f K,j for j M. 
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bi . H, T°/it f f (S, 



NF 

+ I K., T 



j = l 



IJ 



fj) 



( f oo N o /o 

* jf, '=1.i Vj 



- T, 



;) 



For Mi = M® = 0, bj reduces to 



N 



f S. + Z K. . T. 

J=> ^ 



from Eq. (1), and f cancels out of the equation. 



SOLUTION 



Temperatures are determined using Eq. (3) and a standard linear-equation solving technique. 
In this case the Los Alamos Scientific Laborator>'’5 (LASL) lil^rary routine LSS was used to solve 
the variable temperatures, after the values of the a and b arrays had been determined. Because 
many of the variables used to determine the a and b arrays can be functions of temperature, thus 
making the problem nonlinear, an iteration loop is provided at each time step to obtain tem- 
perature dependencies where the temperature determined from the previous time step is not 
satisfactory. This problem arises particularly for natural convection (thermocirculation or ther- 
mosiphon) connections and thermal radiation connections. 

If the coefficients a,; in Fa\. (3) are constant — meiining the conductances are cemstant — the ^et 
of linear eciuations can be solved at the beginning of the program run. This results in a set of coef- 
ficients Fij, such that the node temperatures can then be determined l)y 



h 



NV 

2 



b.; i . 1. tIV, 



(4) 



M 




where the b arrays are as dezicrihecl in &|. {[\). 

Use of thi> option can reduce computation time considerably for nunlels with numerous nodes. 
In l^\SOLE, four sets of Fu’s are calculated — for daytime or ni;^httime o|ieration with a fixed or 
floating ’’control node” temperature. 

CONTROL TEMPERATURE 

A particular temperature, such as room temperature, is monitored by a control node to main- 
tain it within t)rescribed hounds. At j)re>eiU, only one control node is allowed. When this tem- 
perature is within the prescribed bounds, the cont rol node is a variable-temperature nude and no 
auxiliary heat source is required. When the control node temperature goes beyond a temperature 
bound during a time increment, the control temperature is fixed at the bound and an aa.xiliary 
heat source is calculated to maintain that temperature. 

H 

= 2 S-j ("i - Tj) - 
j=1 

where T| is the fixed bound temperature of the control node. In some models, the thermal con- 
ductance value of a connection to the control node can be varied within certain bounds (.such as 
for proportional control dampers). This variability reduces or eliminates the need for nn aux- 
iliary heat source to maintain a control node teinj^erature bound. 

When the operating mode of the control node changes during the time step, the time step i.s 
subdivided and two or more calculations are made. Running sums of positive (source) and 
negative (sink) auxiliary heat sources are made. The sink (cooling) calculation is divided f :.'"her 
into "ventilation” and “auxiliary” cooling, depending on whether the cooling occurs when ‘ ::e out- 
side air temperature is below or above the control node upper bound temperature. 

Running sums aUo are kept of the total heat into and out of individual nodes, the total node 
licat source and sink other than auxiliary, and total connection heat flow.s. For bookkeeping })ur- 
poses the auxiliary heating or cooling of the control node is not combined with any other sources 
for that node. 



THERMAL CONNECTIONS 

The most common thermal connection is the ”UA" type 

Sj ■ \ <« 

where Ux is the overall iieat transfer coefficient for connection k between nodes i and j. and is 
the heat tran.sfer area of connection k. Many niodel.s have an advection connection that results 
from I hermocirculatinn. The algebraic fc^r.m of this connection conductance depends on the 
assumed temperature distribution.s in ti:e leg.s of the circulation path and tlie relation between 
the temperat ure distributions and the temperature of the nodes that represent these legs. A com- 
mon example is a Tromhe wall collector that has an advection connection l^etween the room and 
the healed air s[^aco l)etween the wall and gla/.ing. If it is assumed that t he temperature cli>tribu- 
tion in the air space is linear, such that tht" air space node i> at the arithmetic average of the inlet 



and outlet temperatures, and that the room air is, fully mixed, the thermal conductance is ^iven 
by Kij = 2pCpV, where C,, is the air specific heal, p is the air density, and v is the volumetric air 
flow rate. Because the value of v depends on the a\erage air column lemperatures, the thermal 
conductance is hij<hly teiu[)erature dependent and iterations are required. 

In the example model for a Trombe wall, the therniocirculation flow is calculated by 

V = >/gH AT/T , (7) 

where 

Cd is the vent discharge coefficient— taken to be 0.8, 

Av is the area of one row of thermocirculation vents — top or bottom (assumed to be the same), 
g is the acceleration due to gravity, 

H is the air column height, 

AT is the difference between the average temperatures of the wall/glazing air space and the 
room behind the wall, and 

T is the average absolute temperature of the air space. 

Equation (7) was derived by considering a driving force caused by the density difference between 
the two air columns and a flow resistance from two rows of vents in series, assuming the vents are 
the dominant resistance to air flow. The equation must be modified if there is another significant 
flow resistance, such as that to Bow in the wall/glazing air space. 

For problems involving one-dimensional transient conduction, such as in thick masonry mass 
walls, an adequate model can be made using several internal nodes in series, in which each node 
has a part of the total mass associated with it, and massless nodes on each surface. The nodes are 
connected by appropriate conduction U.A values in the direction of heat flow. Other schemes in 
which each node has some mass associated with it may also give good results. 



INSOLATION 

Any node may have sever^al solar heat sources with various insolation areas, glazing tilts and 
azimuths, number and thickness of glazings, and solar absorptivites. In some cases, the ap- 
plicable external insolation is read directly from the weather data file; in other cases, the weather 
data insolation is for a horizontal surface and a correction to the insolator orientation is required. 
Also required are estimates of reflected insolation and glazing transmission. The equations used 
to determine the insolation on an internal solar collector surface are given in Appendix A (most 
of the equations are from Ref. 1). The solar sources are added to other heat sources that have 
been specified for that node. The solar source for node i is calculated by 

S. = QTRAN.AGLZ-ALFA , (S) 

where QTRAN is the solar flux transmitted through the glazing, AGLZ is die area of absorbing 
surface associated with node i, and ALFA is the .^olar absorhtivity of the surface. Solar radiatiiM'i 
not absorbed by the surface is assumed to be lost; that is, internal reflections are not accounted 

for. 

In the present P.ASOLE model, glazings of the primary solar source are represented by nodes. 
These nodes are given heat .«Joiirces equal to the solar radiation absorbed by the glazing. 'Phe 
model could be expanded to allow for heat source nodes in the glazings of all solar sources. 



PHASK-CliANGE MATERIALS 



Because temperature is the primary dependent variaole, it is not coin'enient to handle hear-of- 
fusion as an isoth.crmal enthalpy change. Problems involving i)hase-change materials haw hcren 
solved by representing the heat-of-fusion as an iiiLTea>e in heat capacity over a given tern* 
peralure range of about 10 to 20'^C. For these problems the heat capacitance becomes a 
temperature-dependent parameter. 



PROGRAM STRUCTURE 

PA30LE consists of a main program, seven subroutines, and one function (Fig. 1). All but 
one (COLLECT) of the subroutines are called directly by the main program. In addition to the 
PASOLE subroutines, several routines in the LASL computer librarv' also are called by the main 
program. 

Three time step loops are set up in the main prf)gram. The outer loop advances the month, the 
next loop advance.s the day of the month, and the inner loop advances the basic time 
step — usually 1 h. Another loop inside the inner time step loop is used for iterations, if recjuired. 
for temperature clej)endence of the coelficients. Coefficients a,j and b, of Eq. (3) are computed in- 
side the iteration loop. Equation (3) is solved using the L.ASL library routine LSS, or the tem- 
peratures are determined directly by solving Eq. (4). If any of the calculated node temperatures 
deviate from their previous iteration value (or from the previous time step, if the first iterati(»n) 
by more than the specified value of TOLT, another iteration is made using the calculated tem- 
peratures as the new values. The iteration j)rocess continues until convergence is reacheci or until 
ITM.A.X iterations have been made. If convergence is not reached the run is not stopped, hut the 
output counter IvERR is advanced by one. 

If the control n()de changes operating mode (from fixed to variable temperature, or vice versa) 
during a time step, the basic time step is divided into twoor more smaller time steps, but no new 
solar/weather data are obtained. 

At the end of each time step, heat flow sums are updated, \-ariahles are set for the next time 
step, and film fie variables are set if graphic output is to he obtained. .After all the time step 
loops are completed, the summary output is made. For parametric studies the entire pn^gram is 
put in a "problem” loop, in which specific parameters are varied. 



SUBROUTIMvS 

Tlie ilata required to run the problem and the parameters (Apj)ondix B) that describe the 
model structure are set in subroutine INDAT.A. We also could use INDATA to read data from 
cards or from another file. .All data written into INDAT.A are given in FORl'IUAN. Except for a 
block of default values at the beginning, IND.ATA probably vrould be changed completely for 
each ditferenc model. 

IVeliminary calculations are done in subroutine PRLME, using the data supi)lied in INDAT.A. 
Included in the>c calculations are the .setting up of arrays of fixed- and floating-temperature 
nodes and the determination of the coefficients Fij of Eq. (4), if apidicahle. 

Subroutine D.AYLY is called once each time through the day loop. Here, the daily weather 
data are read from a data file or are calculated. In addition, the solar declination and ex- 
Iratcrrestrial normal solar flux are computed for the current day. Weather data required are the 
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total S'^iar radiation Hux on a known tilt (QHD) and the ambient air temperature n'AD). Also 
useful is the wind velocity (VHLD). If there are rixed-temperature nodes other than the outside 
ambient, DAYLV is a convenient place to enter the ternf)erature data for these nodes. Ijecau.-r 
weather data formats are different, lhi,> subroutine probably will have to be adapted to the i;ar- 
licular data set being used. 

SUNSRC calculates the heat sources from solar radiation absorption. Most of the equations in 
Appendix A are solved by SUNSRC. If the solar radiation measurements read in DAYLY ar- iioi 
taken on the horizontal plane, the subroutine COLLECT is called. COLLECT’ solves the 
correlation backwards; that is, it obtains the equivalent total horizontal radiation from the ’iiven 
measured radiation, the till (TILTC), and azimuth (.AZIMC) of the measuring surface and the 
a.ssumed diffuse ground reflectance (RHOC) pertaining to the measurement. The computed 
total horizontal radiation is then used to proceed with the calculations. Different optical <y .sterns 
than those assumed here (for example, one with internal reflections) would require change^ in 
SUNSRC. 

Model structure parameters that should be furnished by INDATA, hut that are temperature 
and/or time dependent, are calculated in subroutine PROP, which is inside the temperature 
iteration loop. PROP is called just before the soluti^m of E({. (3) for node temperatures is ob- 
taiiied. This subroutine, like INDATA, probably will require changes for each different mo«.ieL 

Subroutine CONTROL determines whether a change has taken jdace in the operating mode of 
the control node. If a change has taken place the flags KICHNG and KIC are set, the fixed and 
variable node number arrays are adjusted, the time at which the node change occurs is deter- 
mined, and temperatures at this time are calculated by linear interpolation. A new time incre- 
ment is then determined for the remainder of the original time increment. 

Descriptions of the variables found in the common blocks are given in Appendix C. 

COMPUTER SYSTT.M REQUIREMENTS • 

Most of the LASL library routines, such as DATEl, CLOCKl, LSS, SPLOT, PLOT, and 
• DLCH, will have to be replaced by local equivalents. D.ATEl and CLOCKl, which merely return 
real date and time values for output, could be eliminated. SPLOT, PLOT, and DLCH are used 
to make SC-4020 CRT' files for graphic output. If no equivalent routines exist or if graphic output 
is not required, that .section of the program may also be eliminated. The routine LSS, called from 
the main program and from the subroutine PRIME, is used to solve the set of linear equations, 
Eq. (3), for the Tj's given the coefficients an and bi. Any computer facility of reasonable size 
should have software equivalent to LSS. 

PASOIT*! requires a field length of 57 000 octal words for compilation and 122 000 octal words 
at e.xecution including all library routines. No attempt has been made t(> minimize the core- 
memory. The problems shoum in Appendix E and in Figs. 2 and 3 require 4 to 5 min of CP time 
on a CDC 6G00, or about 1 min on a CT>C TOGO, for a yearly calculation. I’.sing the linearized 
method of Eq. (4), the executi«>n time is reduced to about 1.5 min on a CDC nnOO, 

The loaders used at LASL automatically clear the registers, a feature used in PASOLE 
])rogramming. Not all computer .system> are sot up this way, so it may be nect*s^a^y m-ikt; 
special provi.sions to clear the register.s. There are ot her differences between t he CDC FOHTFLAN 
and other computer .systems, such as multiple replacement statement.s and packed FOHTR-AN. 
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Water wall model with massive exterior wall (KWALL = 1, NSGW = 2). 

(9) 




MODELING RESTRICTIONS 

For simulations in the southern hemisphere, changes to the solar angle and time equations are 
required. There are no restrictions on initial temperature conditions, but a small heat balance 
error will be observed if the initial temperatures are not in equilibrium. The initial value of the 
ambient temperature(s) is set to the first %veather data values In SllNSRC. 

The number of nodes is now limited to 50 because of dimensions. If more detail is required, a 
different type of simulation program probably would he better suited. The computer rime re- 
quired is approximately proportional to the square of the number of nodes, A model with 1^0 
nodes takes about 40 min of CP ti me on the CDC GnOO for a yearly hour-by-tiour calculation. This 
time could be reduced considerably, however, if the linearized calculational procedure, Eep (4), 
could be used. 






RX AMPLE MODELS 



}3uilt rnLo subroutines INDAl’A and PROP of the j:)resent PASOLE version are two types of 
south-inass-wall (SMW) pas.'^ive solar heating syslenus: (1) Water wall, in which the thermal 
storage and ab^o^ber svall are assumed to be isothermal and, therefore, representable by a single 
node, and (2) Trombe wall, in which the transient conduction thrcnigh the SMW is modeled l)v 
several nodes in series assuming the heat tlow is one dimensional. Figures 2 and 2 show the ther- 
mal network representation of these two types of SMW.sy>iems. \\‘all type is .selected i^v setting 
the flag KWALL to 1 for a water wall and to 2 for a 'Prondje wall. I'hermocirculation of warm air 
from the wall/gla/.ing air .space into the interior space is controlled by the tlag KVIONT, which is 
set to 0 for no tliermocirculacion, 1 for unlimited therniocircuhUion, 2 fur backdrafc prevention, 
and 4 for thermostatic control of the air llow. d'he lumAier of gla/.ings nuw i)e varied by changing 
the N(iL value. The number of segments in the Trcunhe wall may be varied by changing the 
NSEG value. 

The basic models have heat capacitance associated only with tlie SMWs. Ho\vever, if the 
parameter NSGW is nonzero, an additic>nal heat fiow path is modeled through NSGW' segments 
of a massive external wall, as shown in Fig, 2. This one-diinen'>ional heat How path, like chat 
through the Trombe wall in Fig. 3, consi:>ts of several internal nodes that have the wall mass 
associated with them and mas.>le.ss nodes on the surfaces. 

Parameters peculiar to these models and used only in IND.ATA and PROP are explained in the 
Appendix D program listing. The INDATA listing in Appendix D shows the first block of data is 
the set of default values for various ’’Program Variables." Here, program variables are defined as 
those set in IND.ATA and PROP that are needed to run the program — largely those listed in Ap- 
pendix A. "Model Parameters" are detlned as those u.sed strictly in INDATA and PROP to set up 
the specific model, that is, to calculate values of the program variables. 

After default values are set. the specillc model programming begins with the setting of values 
for the mode! parameters. The.^e values are most often changed when making j)arametric studies 
with the model. More model parameter.s are set further down for each wall type. Progrvam con- 
nection parameters and node parameters then are evaluated using a methodology deveh)ped for 
these models. Next, constants u.^ed in PROF to evaluate temperature-dependent conductances 
are calcuUited. 

The optional nmssive exterior wall is modeled for NSGW > 0 by adding nodes and connections 
to those set above. Next is a section for linearizing tlic model coin[)letc*ly, if so desired, by 
calculating effective constant conductances for all connections. This linenri/.atiot'i must be done 
if the simplified calculation (KC.Al.C ~ 2) method is to be used, as exj)lainecl in Pf\IME. 

A provision for an outside insulation node on the S.M W is used mostly for iionglazed south wall 
calculations. Finally, an output li.^ting is made fur many of the model parameters and resulting 
program varialdes. 

In subroutine PROP the volumetric thermocirculntion flow rate is calculated for EVENT non- 
zero using the i)reviously calculated temperatures, after which conductances based oi^ this flow 
rate are comT)uted. Radiation and convection conductances between glazings, wall, air space, 
and outside air are calculated, again using the temperatures obtained from the t)revioirs itera- 
tion. 'ri'.e tirral >ection of FM^OP. \vhich calculates valuers ot [irogram \'ariables U, CON'D, and 
SCON, .should he retained regardless of the model chosen. Even when all coiuluctani'es are con- 
^tanl, these last calculations are made in PROF^. Only the night (UN) and clay (UD) values of 
conductance are set before these last calculations are made. 
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SAMPLE OUTPUT 



Appendix E shows output listings obtained from running the Appendix D program with the 
specific model parameters of Figs. 2 and 3 and using Los Alamos weather data beginning .Sep- 
tember 1, 1972. The first line of each listing ir\ .Appendix E gives real time and date, and the 
second line gives the model starting date. Parameter values given for "conductance connections" 
are connection number J; node numbers II and 12 (IlCON and I2CON) connected by connection 
J; and UD (day) and UN (night) values of conductance in Btu/h-®F-ft- consistent with .AGON, 
which is the connection heat transfer area per unit SMW glazing area. Positive heat flow in con- 
nection rl is from node II to node 12. Some of the conductances are shown as zero in the listings. 
The zero value usually means that the conductances are temperature-dependent and have not 
been calculated in INDATA. 

Ambient temperature node and solar heat source information is given next. Following that is a 
block of integers and a block of real variables, which are a combination of program variables and 
model parameters whose names correspond to those defined in Appendixes B, C, and D. Except 
for the first line in the listings, the information is generated in IND.ATA when the flag KHEDPR 
is nonzero. 

The summary table, always generated at the end of the problem in the main program, gives 
monthly totals. 

DEGDAY - heating degree days 

QOUT - heat loss back through SMW glazing 

QHEAT - auxiliary heat required by the control node 

QCOOL - total cooling required by the control node 

QSOIA.R - solar radiation absorbed in primary solar source 

QLO.AD - building heat load for other than SMW 

QACL - auxiliary cooling required 

QSINC - solar radiation incident on glazing of SMW 

PCTSOL - percentage of solar heating defined by PCTSOL = 100 (1-QHEAT/QLOAD). 

Headings given in the summ.ar^' table are MTIME: total number of basic time increments (hours 
in this case); NSTEP: total number of time steps calculated; NCALC: total number of tem- 
perature solutions (including iterations); and KERR: number of convergence failures. 

If IPRSM is nonzero, tables of the individual node and connection heat flow sums are printed. 
These sums are defined in ^Appendix C. Heat flows are given in Btii/it^ of SMW glazing. 

In addition to the printed output, graphic output can be obtained using LASL library routines. 
Graphic output is generated starting when DATE is equal to lOl through the day before 102. Ex- 
amples of plots for the model of Fig. 3 are shown in Appendix E (Figs. E-1 and E-2) for the same 
Los Alamos weather data as above, December 31, 1972, through January 6, 1973: (lOl = 123172. 
102 = 10773). Figure E-1 shows the time variation of the temi)erature.s of nodes 1, 6, 8, and 7. 
Figure E-2 shows the rate of solar radiation absorption in node I (the primary solar heat source), 
the rate of heat flow thrcuigh connection 8 (the heat load), the total heat flow into the 
through connections 7 and 9, and the total heal lo.-ss from the stnith side thrtmgh connc'crions 1 
and 11. 

Further results obtained using PASOLE may be found in Refs. 3 and 4. 
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APPENDIX A 

EQUATIONS FOR SOf.AR RADIATION 
1. Solar declination (deg)^ 

DEC = 23.279 - cos (30 - (MO ~ 1 + DAY/32 - 5.7)], 
where MO - month of year (MO - 6 for June, etc.) and 
DAY ~ day of month. 

'2. Extraterrestrial normal solar radiation' 

QSP = A1 - A2 • .sin [360 - (272.1 + IDAY)/365|, 
where Al = 426.93, A2 = 13.50 for QSP in Btu/b-ft^, 

A1 = 1346.1, A2 = 42.56 for QSP in W/inh and 
IDAY = d:iy of year. 

3. Equation of lime (EOT) is given in tabular form, one value per month (min) 
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4. Solar time (!i) 

SUNTIME = TIME + EOT/60 - DLONG/15 ~ 0.5. 
where TIME = local standard time, and 

DI.ONG = diiference between local longitude and standard time meridian (deg). 

(0.5 is subtracted when insolation data are averaged over the preceding hour.) 

5. Hour angle (deg) 

HR - 15 - (12 ~ SUNTIME) 

6. Critical hour angle (hour angle for which solar azimuth is 90^) 

COS(HRCRIT) = TAN (DEC)/TAN(LAT) 

with constraint that 0 < COS(HRCRIT) < 1. 

7. Sun altitude (deg) 

SIN(ALT) = COS(L.-\T)^COS(DEC)-COS(HR)TSIN(L\T)-SlN(DEC), 
where LAT = latitude (deg). 

8. Sun azimuth (deg) 

SIN(AZr) = COS(DEC)-SIN(HR)/COS(ALT) 

AZr will always be between -90'' and +90®. 

To find the tnie value of the azimuth. AZIM 

If |HR| ^ HRCRIT, AZIM = AZI\ and 

if |HR| > HRCRIT, AZIM = (ISO- jAZI'D - AZIV]AZri. 

9. Does correlation^ for direct normal and diffuse radiation measurements 

(a) PP = QH/[QSP^SIN(ALT)1, 

where Qhl is total horizontal radiation. 

(b) FQUN = l.TD^PP - 0.55 with constraint 
that 0 < FQDN < 1.0. 

(c) QDN = FQDN‘QDNMAX. 

where QDN = direct normal radiation and 
QDNMAX = 1000 W/m'^ = 317.2 

(d) QDIF = QH - QDN^SIN(ALT), 
where QDIF = sky diffuse radiation. 

10. Wall/solar azimuth (deg) 

GAM = AZIM - WAZI, 

where WAZI is the wall azimuth; 0 when facing south and positive when facing east of sovith. 

11. Wall/solar angle of incidence (deg) 

COS(INC) = COS(ALT)'SIN(TILT)*COS(GAM) + SIN(ALT)*COS(TILT). where TILT 
is wall tilt from horizontal — [jositive toward south, and 0 < INC < 90®. 

12. Radiation incident on wall from specular reflector. The equations are only for a horizontal 
reOcctor adjacent to a vertical wall. East-west dimensions of wall and reflector assumed 
equal. See also item 14. 

(a) RLKFF - C03(C;AM)/TAN(ALT). 

If RLEFF > KLNGTI I. RLKFF is set equal to RLNGTH, where RLNGTH is the ratio of 
the reflector length (N-S) to the wall height. 



{!)) DWOW - RLE[T-TAN(GAM)/ASPRAT, where ASFIUT is the ratio of the reHector 
width (EAV) to its length (N-S). 

(c) ASR = RLEFF-(l-D\V0\V/2) if DWOW < 1, 

ASR = 0.5^^SPR.■^T/TAN(GA^^) if DWOW > 1. 

(d) QINSR = QDN *ASR'RHOSR'*SlN(ALT), where QINSR = reflected radiation incident 
on wall, RHOSR = reflectivity of retlector. 

13. Shading from overhang^ 

AFAC r is the ratio of direct normal radiation incident on wall with overhang to that without 
overhang (Fig. A-1). 

AFAGT = 1 - [OHANG^TAN(BEFF) - 0SEPR]/[S1N(TILT) + COS(ldLT) 

-TAN (BEFF)] 

Constrained: 0 < AFACT ^ 1, where TAN(BEFF) = TAN(ALT)/COS(GAM). 



14. Total radiation incident on wall (collector) 

QINC = QINDN -r QINSR + QINDF + QINRF. 

(a) QINDN is incident direct radiation 
QfNUN = QDN-AFACT^COS(INC). 

(b) QINSR is incident retlected direct radiation (see item 12. d). 

(c) QINDF is incident sky diffuse radiation 
QINDF = QDIF -[1 -h COS (Tl[/r ))/2 
(assuming infinite horizon). 

(d) QINRP' is incident reflected diffuse — includes diffuse ground reflection w'ith infinite 
horizon and specular reilection of sky diffuse radiation. 

QINRF - QH-RHO^{0..a-(l - C05(TILT)1 - FCSR) + CJDIF-FCSR-RHOSR. 
where RHO = diifuse ground reflectivity, FCSR = view factor between wall and 
specular retlector. 

P'or this retlector model (see item 13) 

FCSR = 0.5 (RLXGTH + 1 -yjKLSCTlV + 1) appro.ximately. 

tAs-siimed to be infiiiite in K-W Hinirnsion (no end ettects). 




OHANG = B/A 
OSEPR = C/A 



Fig. A- 1. 

Overhang geometry. 



15. Transmission through j^lazin^s 

QTRAN = [(QIXDN -f QINSR) -TRANS^GL.M^S f- (QINDF + QINRK)- 
D'ni\NS-DGLABSl'M*RCOAT 

(a) TRANS = transmission resulting from surface reflection losses only 
TRANS = 0.5MTI •}- T2), 
where Tl = (1 - Rl)/[1 + Rl*(2"NGL - 1)], 

T2 = (I - R2)/(l + R2'(2 vNGL - 1)], 

R1 = (SINdNC - ARP:F)/SIN(1NC -f AREF)1^ 

R2 = (TANdNC - AREF)/TAN(1NC + AREF)1A 
SIN(ARKF) = SINdNO/RlNDRX. 

NGL = number of glazings in series, and 
RINDEX = refractive index of glazing surfaces. 

GLASS = transmission resulting from absorption in the glazings only^ 

GRABS = EXP [-EX*TGLZ-NGL/COS(AREF)K 
where EX = glazing material extinction coefficient., and 
TGLZ = thickness of one glazing layer. 

TRCO.AT = transmission of one glazing coating. 

TRANS and GRABS are for direct or beam radiation. An approximation to the 
transmission for diffuse radiation (DTR.ANS and DGLABS) is found by setting the 
angle of incidence (INC) to 60® in (a) and (b) above. 



(b) 



(c) 

(d) 



tThe variable CI.ABS in Appendix D is transmission for a sin{<!e glazing layer. 



z\PPENDIX B 

BASIC INPUT DATA , 

The basic input data must be entered into the program by the IND.ATA and PROP sub- 
routines. Some parameters have internally set ’'default" values (shown in brackets). Use of the 
units must be consistent; either ASHRAE units (Btu, ft, h, ®FJ or S.I. units (W, m, s, ®0) may be 
used with no significant changes. Frequently, many of the data are calculated from other input 
data in the subrcjutine INDATA, which is executed at the beginning of the problem. Some of t he 
data (such as U values) may vary during the problem. In this case the data are recalculated in 
the subroutine PROP, which is svithin the time and iteration loops. 

(A) Nodes (I) — 1 values not necessarily contiguous, but max I must be < NM.\X, which is 
currently 50. 

KF(I): 0 for no node 

1 for variable-temperature node [O] 

2 for fixed-temperature node 

T(I): Node temperature — initial values must be specified for all nodes. For 

fixed-teiiqierature node?, 'FiD imi?r. lie ?pecifitrd fnr all time step^. 

CPM(I): Mass-heat capacity [heat capacitance, M in Eq. (1)1 

S(I) Heat source = 5(1) -r SP iD^TG) [0] 

SP(I) See Eq. (1) 

(B) Connections (J) — J value.? .should be contigaious starting \vith 1. 

NCONN: Number of connections (Js) 



I)C.ON(fl)) connected by thermal connection J 

]-^CON(d)/ 

UD(J)1 Day and night specific conductance 

UN(J)/ values of connection 

ACON(J): Heat flow area of connection 

U(J) is set to UD(J), ifTMORN < TIM?: < TEVEN; UN(d) othenvise. 
Then the conductance is COND (Il,I2) = U(J) ‘•ACON(J), 
where lI = IlCON(J), 12 = I2CON(J). 



(C) Solar heat sources 

NJS: Number of solar heat sources 

ISOL(J): Node number of solar heat source J 

TILT(J): Collector tilt angle from horizontal for solar heat source J (deg) 

WAZIM(J): Collector wall azimuth of source J (deg) 

AL?'A(J): Receiver solar absorptivity for source J 

AGLZ(J): Glazing area for source J 

NGLZ(J): Number of glazings in series for source J 

OHANG(J): Overhang divided by glazing height (see Fig. A-1) 

OSEPR(J): Overhang vertical separation divided by glazing height (see Fig. A-1) 

TGLZ(J): Glazing thickness (feet or meter per layer) 

D.‘\BS(J): Glazing diffuse absorptivity per layer 



(D) Ambient temperature nodes 



NJA: 
lAMB(J): 

(E) Integers 
IC: 

IDPERYR: 

INDATE: 



ion 
102 1 
IPRSM: 
ITMAX: 
KAlkXCL: 



KCALC: 



KCONU: 



KCTY: 

KD.ATA: 



Number of ambient temperature nodes 
Node number of Jth ambient temperature node 



Node number of control node 
Number of days per year 

Initial date for which calculations start — of form MMDDYY, 
where MM is month, DD is day of month, YY is year 
Beginning and ending dates of hourly print and 
film plots (102 is day after last day) 

Nonzero for print of all heat flow sums 
Maximum number of iterations per time step 
= 0 for no auxiliary cooling (that is, no cooling 
when TA.MB > TCMAX) 

5^ 0 for auxiliary cooling 

1 for standard solution (Eq. (3)] 

2 for siTnplified solution [Eq. (4)] 
see comments in PRIME 

0 for variable conductances 

1 for constant conductances 
see comments in PRIM?^ 

Flag for weather data format (see subroutine D.AYLY) 

= 0 for daily data (ambient temperature, wind velocity, 
solar radiation) rend from tape (T.APhX) 

7^ 0 for data from other source 



[365] 



[ 0 ] 

[01 

(11 

m 



[1] 



[01 



[01 
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KDELT: = 0 for no subclrvrding of basic time increment 

0 standard subdividing [1] 

see comments in PRIME 

KFILM: = 0 for no film plots [0] 

5^ 0 film plots 

KHEDPR: = 0 for no printout of ITs’DATA information [1] 

0 print 

MOSHl ] Beginning and ending months of summer period 

MOSH2 ) when night insulation operation is reversed ^ [13, 0] 

NDAYl: Number of sequential daily calculations if < NDAY (MO) [32] 

NDAY(MO): Number of days in month MO [set in data 

statement — if IDPERYR == 366, NDAY(2) is reset to 291 
NHOUR: Number of time increments on each film plot frame [168] 

NMAX: Max number of nodes, that is, dimension [50] 

NMO: Number of successive monthly calculations 

NTIME: Number of time increments per day: At = t/NTIME, 

where At is time increment length and t is total number of time 
units per day (DELTl and TFERDAY) [24] 

PRINTD: = 0 no daily print [0] 

0 daily print 

PRINTH: = 0 no hourly print [0] 

7^ 0 hourly print 

PRINTM: = 0 no monthly print [0] 

7^ 0 monthly print 

UNITS: = 1: Btu, h, ft, "F (ASHRAE units) [I] 

= 2: J, s, m, °C (S.I. units) 



(F) Real variables 
ALAT: 
ASPRAT: 
AZIMC: 
DLONG: 

EX; 

FAC: 

FACI: 

FCSR: 

QDNMAX: 



QHCONV: 

QICUT: 

RDIFF(NGL): 



RESNI: 

RHO: 

RHOC: 



Latitude of locality (deg) 

Specular reflector aspect ratio (EAV to N-S) 
Pyranometer^^ “wall azimuth" (deg) 

Local longitude minus standard time meridan (deg) 
Extinction coefficient of glazing material (ft*0 
f in Eq. (2a) 

Same as FAC for heat (low integrals 
View factor betsveen collector and specular reflector 
Maximum direct normal insolation in direct/diffuse 
correlation for measured data (Btu/h-ft^ or W7nr)13l7.2 
for UNITS = 1, 1000 for UNITS = 21 
Conversion factor from weather data insolation 
Value of incident solar radiation below which night 
insulation is applied 

HeHectance of glazings for NOL glazings in series for 
diffiu^e radiation (0.16, 0.24, 0.29, 0.33, 0.35) 

Resistance of night insulation on collector '^F'/(Htu/hTtM 
or "CAW/m’) 

Diffuse reflectivity of ground or other external reflector 
Same for pyranometcr^^ 



[ 0 ] 

(G.Ol 

10.5) 

[0.51 

[01 



[-1 X 10^1 



[0.3] 

(0.3) 



is npplit*il when TI.MK > TIA'K.N and TIMK ^ TMOR.N h*r winder 
ttl’yranonurter for which we;»U>fr data utanlation is 



I 



RHOSR; 

RINDRX; 

RLNGTH: 

SBCON: 

TBASE: 

TCM.-VXD; 

TCMAXN: 

TCMIND: 

TCMINN'; 

TE\^N: 

TILTC: 

TMORX: 

TOLT; 

TPERDAY: 

TRCOAT; 

TZERO: 



Refleclivity of specular reflector 
Refractive index of plazinj^s 

Specular rellector length (N-S) divided hy glaring height 
Stofan-Boltzmann constant (1.7Io2 x 10'’ Btii/h-fl’-“F* 
or 5.67 X 10-* W/m’-“C‘) 

Base temperature for heating degree-days (65°F or 1S®C) 

Max allowable control ncKle temperature for day 

Max allowable control node temperature for night 

Min allowable control node temperature for day (°F or °C) 

Min allowable control node temperature for night ("F or °C) 

Time at which "night" starts 

Tilt of pyranometer^ (degree from horizontal) 

Time at svhich "daytime" start.s 

Tolerance on temperature iteration (°F or °C) 

Number of time units per day 
Transmission of coating on inner surface of 
inner collector glazing 

Difference between zero and absolute zero on 
temperature scales 



(Ul 

( 1.5261 

(01 



(17 h| 
[ 0 ] 
(7 hi 
(l.OI 

24 h/day or 86 400 s/dayl 
(1.0] 

(460“F or 273“C] 



tP\*ranoine'ler for which \seaiher data insolation is piven. 



APPENDIX C 

PARAMETERS IN CO.MMON 



All parameters in Appendix B are also in COMMON. Starred parameters must be generated 
in D.AYLY (weather data). Routines in brackets are where |)arameter is calculated. 



A(I,J) 

H(f) 

COND(I,.J) 

COSDEC 

COSL.-\T 

COSTC 

DATE 

D.^Y 



DAYl 

DAY2 

DD 

DEC 

DEGR.‘\D 

DEBT 

DELTl 



Temperature coefficients in Eep (3) 

Source terms in Eq. (3) 

'Conductance of connection between nodes I and J 
K„ in Eq. (1) 

COS(DEC) 

C03(AI^T) 

COS(TILTC) 

Integer date of same form as IND.^TE 
(see Appendix B) 

Index of day loojj in main program 
(integer day of month) 

First day in day loop (integer) 

Last day in day loop (integer) 

Heating degree days 

Solar declination (deg) 

r/180 — conversion from degrees to radians 

Current time increment 

Basic time increment 



(M.MN'l 

(MAIN) 



(PROP) 
(DAYLY) 
(PRIME] . 
(PRIME] 

(DAYLY) 

(M.MNl 

(MALN'l 

(MAIN'l 

(DAYLY) 

(DAYLY) 

(MAIN) 

(.MAIN) 

(PRIME) 



ERRT(I) 

FD{I,.J.K) 

FN(l.J.K) 

FRAC 



FSGC 
FSSC 
. HRCRIT 
lAIR^ 

IDAY 

INDAY 

ISOLX(IH)* 



rriME 

IXFILM 

JAIRt 

JCPt 

•JCWt 

JWA^ 



KCOOL 



KIC 



KICHNG 

KND 

KSHUT 



KVENT^ 



Deviation of node temperature between 
successive iterations 

Coefficients Fij in Eq. (4) for "day" and "ni^^ht" 

K = 1 for IC a variable-temperature node; K = 2 for IC 
a fixed-ternperature node 

Fraction of DEL'F since beginning of time increment 

or lust control mode change 

View factor between pyranometer and ground 

View factor between pyranometer and sky 

Hour angle at which azimuth is 90"^ 

(For SMW models) node number of air in 
glazing/wall space 
Day of year 

Initial day of month (from INDATE) 

Weather data time offset 

(see equation for SUNTIME in Appendix D, 

subroutine SUNSRC) 

Index of time increment loop 

Index of film plot arrays 

Number of connection between LMR and IC 

Connection number of nonmass-associated load 

Connection number of mass-associated load 

Connection number between mass-wall surface node 

and lAIR. This is the first in a sequence of 

connections from the wall through the glazings to 

which glazing conductance calculations are keyed 

= 0 for no night insulation 

= 1 for winter operation 

= —1 for summer operation 

(see MOSHl, MOSH2 in Appendix B) 

Control node mode indicator 

-1 for T(IC) at TCMIN limit. QCIN > 0 

0 for TCMIN < T{IC) < TCMAX, QCIN = 0 
+ 1 for T(IC) at TCMAX limit, QCIN < 0 

+ 2 for T(IC) fixed for all time (TCMIND = TCMAXD, 
TCMINN = TCMAXN) 

= 0 for no mode change 
= for mode change 
= 1 for night 

= 2 for day (TMORN < TIMEX < TEVEN) 

= KCOOL for night 
= —KCOOL for (lay 

Night insulation is used if KSHLH' = + 1 
= 0 for no thvrmocirculat ion 

1 for unlimited ihermocirculation 

2 for thermocirculation with backdraft dainj^ers 

4 for thermocirculation with thermostatic control 



I MAIN I 
[PRIME! 



[CONTROLl 

[PRIMEI 

[PRIME] 

[DAYLVj 

[INDATAl 
[PRIME, MAIN] 
[PRIME! 
[DAYLY] 



[MAIN] 

[MAIN] 

[INDATA] 

[INDATAl 

[INDATA] 



[INDATAl 



[MAIN] 



[CONTROL] 



[PRIME] 

[CONTROL! 



(MAINI 

[MAIN] 



[INDATAl 



tFor SMW mndeU. 



KWALL* 


= 1 for water wall model 
= 2 fur Trombe wall model 


[IN DATA] 


LF([F) 


Array of fixed-temperature nodes 


(PRIME! 


LT(IT) 


Array of all nodes 


[PRIME! 


LV(IV) 


Array of variable-temperature nodes 


[PRIME! 


MO 


Current month number (MO = I for January, etc.) 
The month loop index is MONTH, which always 






goes from 1 to NMO 


[MAIN] 


MOl 


Initial value of MO, from INDATE 


[PRIME) 


MTIME 


Number of basic time increments from beginning 






of problem 


[MAIN] 


KERR 


Number of nodes not converging 


[MAIN! 


NF 


Number of fixed-temperature nodes 




NPR1.NPR2 


First and last problem numbers 


[MAIN] 


KPROB 


Index of problem Icn^jp 


[MAIN] 


NT 


Number of nodes, total 




NV 


Number of variable temperature nodes 




PI 


7T 


[MAIN] 


QACL 


Integral auxiliary' cooling of control node 


[MAIN] 


QAHT 


Integral auxiliary* heating of control node 


[MAIN] 


QC12(J) 


Net heat rate through connection J 

[positive from node IlCON(J) to node I2C0N(J)] 


[MAIN] 


QCIN 


Heat source or sink applied to control node to 






maintain temperature limits 


[MAIN] 


QCINO 


QCIN from prevnuus time step 


[MAIN] 


QCINT 


Integral of QCIN over one time step 


[MAIN] 


QCOK(I) 


Net heat rate conducted into node I 


[MAIN] 


QFACO(I) 


Time rate of temperature change of node I, 






previous time step 


[MAIN] 


QHDdH)* 


Weather data insolation (pyranometer data) 


[DAYLY] 


QHZ 


QHD(IH) converted to correct units 


[SUNSRC] 


QINC 


Solar heat flux incident on collector 
(primary solar heat source) 


[SUNSRC] 


QSP 


Normal, extraterrestrial solar radiation 


[DAYLY] 


QSRC(I) 


Net heat source rate in node I 


[MAIN] 


QTRAN 


Solar heat flux transmitted through collector 






glazing(s) 


[SUNSRC] 


QVCL 


Integral vent cooling of control node 


[MAIN] 


QVHT 


Integral vent heating of control node 


[MAIN] 


RCOK(I)^ 


Array of constants generated in INDATA for use 






in PROP, or just for general use 


[INDATA] 


REMAIN 


Remaining fraction of basic time increment 


[CONTROL] 


imndH)' 


A fourth weather data parameter (relative 






humidity in some cases) 


[DAYLY] 


SCON(I) 


Sum of conductances (COND) of connections to 






node I 


[PRIME, PROP] 


SINDEC 


SIN(DEC) 


[DAYLY] 


SINLAT 


SIN(ALA'D 


[PRIME] 



tFor SM\S^ nuuleln. 



SINTC 


RIN'(TILTC) 


[PRIME] 


SQCONl(I), 


InlegraU of positive and negative values 




SQCON2(l) 


of QCON(I) over time 


[MAIN'l 


SQC12K1), 


Integrals of positive and negative values 




SQC 122(1) 


of (jCr2(I) over time. 


[MAINI 


SQSKCl(I). 


Integrals of positive and negative values 




SQSKC2(I) 


of (JSRC(I) over time. 


[MAIN] 


TAD(IH)- 


Weather data — anibient temperature 


[DAYLVl 


TAMB 


Ambient temperature 


[SUNSRCI 


TBAR 


Daily average ambient temperature 


[DAYLY] 


TCMAX 


Ma.ximum allo'^vable control node temperature 


[MAIN'I 


TCMIN 


Minimum allowable control node temperature 


[MAIN'l 


TCOOL 


Threshold ambient temperature above which any cooling 
of control node cannot be done by venting. Now set 






G(jual to TCM.AX 


[MAIN) 


miE 


I'iine of day in consistent units 


[MAINI 


TIMEX 


Time of day in hours 


[MAINl 


TMAX 


Daily ma.ximum ambient temperature 


[DAYLY 1 


TMIN 


Daily minimum ambient temperature 


[DAYLY] 


TO(I) 


Node temperature from previous time step 


[MAIN] 


U(J) 


Conductance per unit area of connection J 


[PROP] 


VELD(IH)* 


Weather data — wind velocity 


[DAYLY] 


VOLF^ 


Therrnocirculation volumetric flow rate 


[PROPI 


YR 


Two-digit year number 


[ DAYLY 1 


tFor SMW models. 


• 
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APPENDIX E 
SAMPLE OUTPUT 
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- ph: 


i I VE 
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9 u 1 73 
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IE 




uii 


UM 


MCDf^ 
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er 
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Zf 
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1 . niirntE-*- II H 
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4 
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1 . 0 ij 1 1 iiE + 11 11 


5 


1 


4 


1.1 . 




u. 


1 . 0 1.1 1 1 i .iE 1.1 1.1 


E. 


4 


cr 


u. 




0. 


1 - 11 U 1.1 1.1 E 11 11 


•y 

f 


c 


E* 


1.1 . 




0. 


1 . uui.ii'E-i-i'iM 


3 


A 


3 


0 . 




1*1. 


1 • 1 j 1 j 0 11 E + I’.i 0 


9 


E* 


( 


1 . 


I'.i 1*1 1*1 E > 1 1 1.1 


1 . 5 11 0 11 E '1' U 1*1 


, 1*1 11 II r»f^ + I'l fl 


1 0 


“7 

f 


3 


•“< 


3 U 0 - I E ^ 0 1’l 


3 . c 1.1 11 1 'E + 11 


Ui'.iiJiiE +• Ui*i 


n 


Cj 


9 


1 . 


rf. M I.IME l.in 


1 . E 1.1 1) i>£ + ' 1 H 


S’. 0 u 1.1 1 'E> 1.1 »!i 
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1 
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1 0 




1 . 


ZIZIE-Ol 


1 . ZlZlE-m 


S - 0 11 1'l OE + 0 0 



ULOhD= 1 . i.njO 



hMETEMT TEMFERmTUPE UCjbE‘1 
I = 

SDLhP HEhT SnUPCEi: 

J I SUL nGLZ TILT MhZ 



1 1 


3 3 I'.I . 


1.1 r.i 


NMmX 


IC 


MTinE 


5i;i 


C. 


Z4 


Ki.^hLL 


} iIhLC 


^ I'ELT 
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1 


1 


JCP 
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3 


13 


cr 



: hLFh 
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3 1 j 
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y L TV 
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T HIjVilL h VEMT 
1 0 
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0 . 



OHlDMV 


hLhT 


DLUnG 


PiniiEX 
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3 . 0 1'.i i.'i 0 E - l.l 1 
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1 . 0 Ij 0 i.’E > l.l l.l 


TiSnUiM 


T c n I nil 
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5 0 0 0 E + 0 1 


Z , 5 tj l.l 0 E 1 1 1 
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3 . 1.1 i'.i 0 0 E ~ 0 1 


T-. iji'iij I'F-O 1 


TCLm.-Ti 


THCnn 


THICK. 


VDLSP 


ULDHTi 


PESni 


T* 0 0 1.1 E + i.'i 1 


0 . 


ij . 


0. 


1 • 0 (i 0 l.l E + l.l 0 


0. 


PLMGTH 


PHOIP 


H SPPhT 


FCZP 


CPUP 


PSni.ii 
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ij . 


(1 . 


4 . 5 0 0 OE 0 1 


0. 


PMDU^ 


iji.inuT 


THCM 


THi^ l.l 
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Hl.l 
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"3 . l.l 0 IJ ijE ♦ ij 0 






4 ' 



-f. 

J 



..j* 





■i 






MT me 



inFP= l jTP 1 KEPP = 



iia I'd 


K Tim* 


ijOUT 


OHEPT 


oc nni 


ornLHp 


1 " LClHl' 


C'hCl 


1 I uc 


pct:dl 




E' Ov 


*r* .'• *Z‘ o 


0 


-)37^ 




L* C. J tZ z* 


1 03 3 7 


— !■ ^■' 


::4. i.i 1 3 : 


1 0 


1 0 


S I I 


1 1 (' 1 4 


533^. 


-x34 


r3x-3 3 


1 5^.35 


(1 


3-^3 


r:* . 1 0 


1 1 


1 0 :-c 


15417 


1 4 r> X. 




1 


34 4. IS 


35 4.5 


II 


XX 713 


4 1 . 4-:* 


IE* 


1 144 


1754? 


1 4.3 1 H 




1 


33F43 


3xr«7l 


0 


4 I 07 7 


43. 35 


1 


1 .^0^ 


1X5::P 


157XX 




1 


.:’.34.S4 


33534 


1:1 


44 . 44 ^ 


4»:.. 5 3 


■Zf 

1 — 




14x57 


1 xx-^3 




1 


35344 


E* 1 . 1 7 4. 


0 


3''. 7 3 4* ;. 


43. 


3 


0 


14i;r^i:i 


1 4.S 1 1 




1 


c* 1 1 5 ?. 


3 j.x; 7 X 


0 


Xo 1 :i:4 


5 1' . 35 


4 


75r* 


1 1 :f: 1 X 


1 1447 




1 


1 •? 1 4 1 


I'iMS;:: 


\\ 


Xl“.’4. 1 iz 


3; 7 . 


5 


4* 'd- 


755 0 


4’c’ o J 




1 


147 14 


1 114 1 ri 




35543 


»7;x. 4..:; 


r» 


\a< 


55 X 0 


■-. r 
•J. J* 


-. c* "-*■ 

C. • 1. 




1 44.74 


0 5 ' * 'T* 


“130X 


C- 

C. • d* • 


34 . 57 


? 


E’o 


^ 


0 


- 3 1 3 


1 


1 :.m4 1 


3331 


- 1 X j: 1 


334.43 


1 0 0 . 0 0 


o 


1 6 


5V79 


0 


~4-^:-:l 


13517 


4.343 


“ i!L 'T‘ 4 7 


■-• s ^ s 


1 0 0 . 0 


SUI):S 


7:-: 50 


137E4X 


’? 7 4 3 


“1373: 




3551 13 


30434 1 


— ool 1 


3i*?75 1 5 


53. 44 . 


TDThLS 




UVHT 


OVCL 






mhHT 


•?mCL 




nn 






0. 




4 li'jE + 


Ox *?. 




'4 37:E-»- 04 


“6 • 1 03E 1 


;i3 7.3; 


;5 05E 03% 





1 


S 14 SPC 1 






.303 PC3 


SOlDHI 


OOODME 


1 


3. 551 1 E-^05 




1 ) , 




0. “3. 


544 .OE + O 5 


3 


0. 




0. 




9 . '•:'*4 04E**- 05 “4.. 


51 03E+03 


4 


0. 




0. 




1 . 97X3E-09 “E. 


5 0 07 E- 1 II 


sr 


1 • 94. 1 4E-»“04 




II . 




7 • 9 9 7* 4 E “ 1 1 1 “ 1 . 


9^:. 79E4 U4 


t. 


3. 443ZE + "4 




0 . 




0. 4. i::-.3E-^ 01 -3. 


44 . 1 nE-i-04 


7 


0. 




0. 




3*‘. 9 1 Oi'E — 04 0. 




c* 


0. 




0. 




9. 573’. 1 E + 0 -x. 


5 055E-I- 0 7 




0 . 




0. 




9 • iT.’ ^ c* 3 E i.i “ 9 . 


35-5R^07 


1 0 


0 . 




0. 




1.5::4.xE-10 “3. 


737 7E-I-0 0 




0. 




u. 




1 • 3 7 3.7 E 04 “9 . 


7 4 3; 4 E + I.I 4 


j 


I 1 13 




30 


l: 1 3 1 


300 1 33 




1 


1 5 


c- . 


S ‘I* 


3E+ 04 


- 3: , 5 . 3 ; 0 \']^ m3 




3 


1 3 


1 . 


1 73 


OE 05 


-1 .x455F^>i« 7 








1 . 


44 0 


3E-I-05 


“1 . 4.573E + 07 




4 


4 3 


0. 






0 . 




er 


1 4 


cr 


1 9 II 


5E-^04 


“ 1 . 95 74.E-I- 03 




r> 


4 5 


c 


1 9 '.I 


5E + 04 


-1 . 95 7EE>i'3 




”7 

/ 


5 4. 


1 , 


^ r 


9 E + 1 . 1 5 


“ '< . 9 I.I 7 4 E -I' 0 1 1 




3 . 


4* 


1 . 


c*. 1 4 


::E+05 


11 . 




3 


C. t" 


K* . 


141 


7E ^ 1.14 


■“1 . •’ C. 4 E I.I .*• 




1 (1 


7 X' 


4. . 


141 


9E-I- 04 


- 1 . 734.4E + II 7 




1 1 


0 9 


cr 


971 


7E-^04 


-9. 7E75E+M1 




13 


9 1 0 


, 


0 7 0 


E "1“ I.I 4 


- 1 . 73 4 xE”!- I.I 7 




1 3 


1 0 3 


4. . 


09 1 


OE 4 04 


“ 1 . Z'34xE + 0 7 





Er^D DF FILE 



16 



k 



PhIOLE 


- ph: 


: I V’ E 


lULHP EriEFCV lZ.nC.ZZ 


73 ' 0 :Z ' 


TPDMFE 


MmLL 


- r-EP i.itUT uPLL 


hPFm - IfU'HTE 


~ Z"17Z 


CDIU'UC 


iHr^ilE 


>: cfiMF.i: T 1 ntii 






J 


1 1 


IZ 


UD 


ur< 


HL 0r« 


1 


1 


1 0 


u. 


0. 


1 . 0 0 I I 0^4^ 0 0 


z 


1 


c 


5 . j: :-E i*' M 


S . Z Z Z ZE -I- 0 1*1 


1 . 1.1 1 ' 1 1 hE 4- n I I 


'S 


z 




Z . C^-o7E > U 0 


c.’ . ^ ^ t' E ”*■ 1 1 1 1 


1 . 0 I I I I I I E -4 0 ij 


4 


:2: 


4 


Z. CCE*7E + ”H 


Z . E-'r-^-7E-^ 1111 


1 , 1 j 1*1 r» 1 1 E 4 1 1 0 


5 


4 


.1 


Z . CCC'i'F + " 0 


Z . CC»r*7E-** 11 11 


1 . U 1.1 1 1 1 >£ 4 n 0 


c 


cr 


♦rv 


. Z Z :• :‘E 0 1'« 


. Z Z Z ZE -♦* 0 11 


1 . 1 1 0 1 III E 4 II 0 


7 


c 


7 


1 . S 11 1 1 mE ^ 0 0 


1 . Z 0 0 1 'E + 11 0 


1 • C' 0 11 OE 11 11 




e 


•t; 


1 , 0 u 0 r»E 4* 0 1*1 


\ . 0 0 1.' ij E 1.1 1 1 


1.11 0 II 1*1 E 4 1*1 ij 


♦4 




•7 


0 . 


0. 


1 • 1.1 11 1 * OE 4 1*1 ij 


1 r.i 


1 




0 . 


ij . 


1 . ij 1 1 1 j 11 r. 4 n 1*1 


1 1 




1 11 


1*1 . 


Cl . 


1 • 0 ij ij hE 4 I’l I'l 


IZ 


1 1.*1 


1 1 


0 . 


0. 


1 . Cl 0 1*1 1 j E 4 n 0 


IZ 


1 1 


X 


0. 


u. 


1 . Cl 1.1 1 j Cl E 4 ij Cl 



M % 



ULDHri = 



I • U M 



HMBIEMT TEMPEPmThPE r^□DEZ: 
I = o 



OLhP HEmT 


fauFi: 


C " 


J IZDL 


nGLZ 


TILT 1.! 


1 1 


C. 


*Z 1 j . ij 1 j 


riMHX 


1C 


MTIME 


50 


7 


Z4 


KMhLL 


y cmlc 


PDELT 




1 


1 


JCP 


jr.u\ 


JMH 


:2< 


Cl 


1 Ij 



Hi hLFh hi?LZ nhHM'? 
0 0 I . 0 'j ij 1 . i!i 0 ij . 0 " M ' I 



riZEC 


riMD 


MDriYl 


4 


iz 


ZZ 


Kcnnu 


F CTY 


\\ l |ZI*1 


ij 


1 


Cl 


IHIP 


..IH I P 





T»?LZ I'hFv 
U lu4 .u7I*:l 

UfUTS: KrihTPf 

I 1.1 

K hUXUL KVEriT 

1 Z 



OH c or IV 

Z. Z 1 ZOE ^ ui* 
TZEPD 

4 , E* u u I'lE uZ 
T CM inn 
C , so 0 OE -1“ 1 1 1 
TCMMXri 
7 • 5 0 0 11 E 0 1 
PLHCTH 

0 . 



HLhT 

3 • S C' 0 E 1 1 1 
I? DM MHZ 

Z. 1 rzoE-^oz 
TC MIMD 
4. . S n 0 0 E ij 1 

THcnn 

1 , 0 I I I I I I E 1 1 0 

PHD.: F 

0. 



DLDriC 

1 . Z 0 0 OE ^ 0 11 
EHn 

?. . 0 0 1*1 1 'E - 0 1 

TILTC 

0 , 

THIC K 

I . s I I M rtp -1“ 0 M 
mIPPmT 

0 . 



PIMDEX 
1 . 5ZCi^E^-oi.i 

e:<: 

E* . 0 0 1 1 r.ic + 1.1 0 
hZIMC 
i‘i . 

vnLi p 

. U ij 0 OF + 0 1 
FC iP 

M . 



niEPP 

0 . 0 0 0 1 1 

ITMhX 
2: 0 

KFILM 

0 



TMDPr< 

7 . 0 0 1*1 "E 0 u 
FhC 

s . 0 r» 11 i.'E ~ 0 i 
PHno 

Z . 0 0 0 1 j E — 1 j 1 
ULDhD 
1 . 1.1 0 0 OE IJ it 
CPMP 

4 . 5 U 1 j I j E I j 1 



TEVEM 
1 . 7 1.1 0 1*' E 1* 1 
T2LT 

1 . 0 1.1 1 1 0 E + 1.1 0 
Fhi»: I 

7- , 0 0 n i;E - 0 1 
PE :ra 

0 . 

P i MM I 

0 . 



47 





MTTME= 3 


74 0 


MiTEP= 3 


:• 3 


MLHl «: = ix 


4 73 


h FPP = 


M 


MG 


DEG DmY 


onuT 


OHEmT 


Cm’.CGL 


0 : HLhP 


OLUmD 


Ohl L 


0 I* I ML 


FO T : GL 




3 MY 


:i: 1 S 1 


4 x3 


-4314 


33::- ?7 


34 1 4 


-1 1 7 


:>.M1 3 


34 . 


1 0 


jMG 


1 \) 1 :-:y 


733 m 


-3344 


3 Ti II 4 7 


lSii3- 


0 


3‘^ 3ii ? 


51.71 


1 1 


1 mS*3 




15341 


- ::33 


'T- 7 1 


35-l^ 


0 


<< < 1 X 


1*4 




1 144 


1 74 MS 


1 7353 


“373 


33771 


33 0 OS 


0 


4 1 073 


1 0 


1 


13 03 


IS 4 4 0 


1 7 4.33 


-1 04 4 


x34 1 0 


? Or- Mr. 


(1 


4 »:• 4 4 4 


4 3.3 0 


C 


’r*-;’3 


14xsri 


15137 


-431 


c 5 1 3 1 


343 1 4 


0 


3734 ? 


33. 114 


3 


V3 (• 


14 333 


1741? 


— 3 »r 3* 


31 033 


4 4 r. 


0 


3 ?1 34 


33. 71 


4 


75Y 


1 3 1 OS 


13514 


-1 MM 


13131 


1 3 3 4 0 


Ij 


3 0 1 3 


35. Ml 


5 


3o4 


7X4 s 


4::.'3‘4 


-333 


1 4 735 


1 M:r;x3 


- 1 03 


35543 


55. 1 M 


6 


133 


4314 


373 


-SMS? 


147U- 


4 057 


-3 XX M 


C- -» ”7 

C. r ^ 1 


3 ? . ? ? 


7 


33 


4 353 


45 


— 4 3 .**4 


1 3 043 


43r>3 


- ? ]• 4 7 


33443 


33 . 43 


s 


13 


4SIX 


43 


— 7 4ci5 


14413 


5 03 1 


-5513 


33555 


33. 04 


IMS 


73*5 u 


1 334 03 


1 1 0343 


-33540 


355144 


3054 ?.S 


-13335 


33751 5 


44 , 34 



tqthls 




QVHT 




3Vl 


L i->HHT 


|•■‘HLL DD 




0. 






-1 . 4575 


E 0 4 1,1 r.i E: 4 E + ‘’.‘5 - 


1.33 •*_: S E + 1 j 4 7 . 3 s M S E + 0 3 


I 


•Cj 


QSPC 1 




S'OSPl 


3 SOlDMI 


30C GM3 


1 


3.55 


17E-^MS 




0. 


O.34S3E-I-M0 - 3 . 


553x3-^05 


0 

4_ 


0. 






0. 


y , 343::*E^^04 -3. 


3 1 X M F + 0 4 


3 


0. 






0. 


3.7151E+04 -3. 


7 i.’ixl E-»' 04 


4 


0. 






0. 


1 . 3433E-^04 - 1 . 


3537 F + 04 


cr 

.^1 


0. 






u. 


1 . 1333E-I-M4 -1 . 


13S3E+M4 


6 


0. 






0. 


3 • 33 s ?E — 03 - 1 . 


SSxxE-Or:. 


3 


0. 






0. 


3. 31053+05 -7. 


F M 3 


3 


0. 






0, 


3. 1 177E + 03 - 1 . 


3 0 »■ rr- E + M c’ 


1 0 


1 . 34 


14E+04 




l.l • 


3.373:?:E-0l -1 . 


34 }:4E + 04 


1 1 


3 • 44 1 3 E 0 4 






3 . 33 33 E + M 1 — 3 . 


4Xrr.5E+04 


7 


0. 








3 . 3 7 1 •: c + M 4 — 1 . 


1 i>/E;E+ 05 


J 


I 1 


13 




si:‘C 1 3 1 


SQC 1 33 




1 


1 


1 0 


1 • 


01 03E-»-05 


11 , 




3 


1 


c. 


1 . 


3343E+H5 


. 1 37E>E + M4 




3 


iZl 




3» . 


3333E + 04 


-3. :3734E + i»4 




4 


A 


4 


4. 


55 i.i ::E+ 04 


-S . .?7 03E + M '< 




5 


4 


c 


4. 


1 443E + 04 


-1 . ::4SiE-^M.?: 




4 


cr 




4. 


3531 E> 04 


.44 73 3+0 




7 


4 


”7 


4. 


xS3lP. + n4 


-3. 4473E + M ?: 




c» 


7 


C* 


CL • 


M5*r'4L+ MS 


-3. r^roE + o:-: 




3 


•4 


7 




34 M7E ♦- 04 


- L. 344 7E+ M M 




1 0 


1 


•4 


1 . 


1 ?xSE+05 


M. 




1 1 


•4 


1 0 


• 


f’x OxE-*- 04 


— f7. , 47 ‘z* kl E + M *• 




13 


1 0 


1 1 


1 . 


531 I:E + ‘i5 


-X. 3415E+M 1 




1 3 


1 1 


3 


1 . 


?4::;4E + 05 


-1 .444 1E-KM1 





EHli Dr FILE 



J 






I 



LAST OAr r 10t73 











1 

1 
























r* 


:,0.. 




r.. s! 


i-=*AC, 


} 1 














— f 
1 






1 


...V 


, 

j:. <. • 


'.J 


1 


r.i:) 










1 

1 




! 






1 


i 

1 






! 


i 




K 


1 

1 




1 




i 

1 


1 



1 




[i 

111 








i 


1 


1 






i i 

1 !! 




ill 




{i\ 








1 


Ji 


1 


1 


' ' 1 




1 iii 








1 


1 






1 


!' 


\ 


! A 


\ 1' 


\ 


h 


j 


- 


\ 








1 

1 




\ 


1 1! 


-V l! 


t 

'vA. 


Ui\ 






1 




L-. 


! !_' 1 






\ '»■ ,7,; \ 












^ 


ixT. 




1 

! 


V 


i 


1 


1 !l 






1 1 ‘ 


1 




1 

! J 














r^’ 




i{- M’.IO 




[COM r 


:n) 






1 1 






r 


p NOC;. 


3 


rau.N 








1 ! 


1 

1 






i 




A. 


/=k 


1 


L 




'A 


1: Hki\ ! fi 






! 

I 


1 




' 


j 


r^'‘ 


i' ^ 


li 


Vi| 




! 1 




M 


1 


i 


I , 

1 
















1 

1 








j 


L.: 




1 



* » i 4 I *“ 'i 

DAYS 



Fig. E-1. 

Node temperatures for December 31, 1972, through. January 6, 1973 (see Fig. 3). 
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